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Abstract 

Helicopter  main  rotor  smoothing  is  a  maintenance  procedure  that  is  routinely  performed 
to  minimize  destructive  airframe  vibrations  induced  by  non-uniform  mass  and/or 
aerodynamic  distributions  in  the  main  rotor  system.  This  important  task  is  both  time 
consuming  and  expensive,  so  improvements  to  the  process  have  been  long  sought. 
Traditionally,  vibrations  have  been  minimized  by  calculating  adjustments  based  on  an 
assumed  linear  relationship  between  adjustments  and  vibration  response.  In  recent  years, 
artificial  neural  networks  have  been  designed  to  recognize  non-linear  mappings  between 
adjustments  and  vibration  response.  This  research  was  conducted  in  order  observe  the 
mathematical  character  of  the  adjustment  mapping  of  the  Vibration  Management 
Enhancement  Program’s  PC-Ground  Base  System  (PC-GBS).  Flight  data  from  the  UH- 
60,  AH-64A,  and  AH-64D  were  utilized  during  the  course  of  this  study.  What  has  been 
determined  is  that,  in  a  majority  of  situations,  the  neural  networks  of  the  PC-GBS 
produce  adjustments  that  can  be  reproduced  by  a  linear  algorithm,  thus  implying  that  the 
character  of  the  mapping  is  in  fact  linear. 
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A  COMPARISON  OF  MAIN  ROTOR  SMOOTHING  ADJUSTMENTS  USING 


LINEAR  AND  NEURAL  NETWORK  ALGORITHMS 


I.  Introduction 


1.1  General 

Since  the  early  days  of  rotary  wing  aviation,  helicopters  have  been  well  known  for 
their  unique  ability  to  take  off  and  land  vertically  as  well  as  for  their  tendency  to  vibrate 
while  doing  so.  Along  with  providing  lift,  helicopter  main  rotors  are  the  source  of  some 
of  the  most  destructive  vibrations  known  to  the  aircraft  industry.  These  one-per-rev 
vibrations  are  brought  about  primarily  by  non-uniform  rotor  mass  distribution  and  non- 
uniform  aerodynamic  properties  among  the  rotor  blades. 

In  order  to  alleviate  these  destructive  vibrations,  the  helicopter  main  rotor 
smoothing  procedure  is  periodically  performed.  Helicopter  main  rotor  smoothing  (MRS) 
is  a  process  in  which  vibration  magnitude  and  phase  are  recorded  both  on  the  ground  and 
in  flight.  Then,  based  on  a  predetermined  relationship  between  vibrations  and  corrective 
adjustments,  changes  are  determined  for  blade  weights,  pitch  links,  and  trim  tabs.  The 
process  usually  requires  several  flights  to  reduce  the  vibration  to  acceptable  levels  and  is 
therefore  time  consuming  and  expensive. 
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Efforts  have  long  been  made  to  improve  the  main  rotor  smoothing  process  so  as  to 
decrease  the  number  of  required  flights  and  thus,  save  time  and  money.  Over  the  years, 
data  acquisition  equipment  has  become  more  accurate,  dependable,  and  user  friendly. 

In  recent  years,  a  substantial  effort  has  been  made  to  improve  the  computer  software  used 
to  convert  raw  vibration  data  into  corrective  main  rotor  adjustments  [1,2,3].  These 
improvements  are  centered  on  the  abandonment  of  the  linear  assumption  relating  main 
rotor  adjustments  to  vibration  response. 

Due  to  the  complexity  of  the  dynamics  involved  with  a  main  rotor  system, 
simplifications  have  to  be  made  about  the  relationship  between  vibrations  and  corrective 
adjustments.  One  such  simplification  is  to  assume  a  linear  relationship  exists  between 
main  rotor  adjustments  and  vibration  changes.  Most  MRS  algorithms  utilize  the  linear 
assumption  to  calculate  main  rotor  adjustments.  It  is  typical  that  several  flights  are 
required  in  order  to  smooth  vibrations  to  within  acceptable  levels.  The  blame  for  this 
poor  performance  has  recently  been  directed  toward  the  linear  assumption. 

To  this  end,  a  new  paradigm  has  been  introduced  to  the  practice  of  main  rotor 
smoothing.  Artificial  neural  networks  have  been  trained  to  recognize  the  relationship 
between  main  rotor  adjustments  and  their  resulting  change  in  vibrations.  This  new 
algorithm  makes  no  assumptions  about  the  mathematical  character  of  the  solution  space, 
so  it  is  therefore  believed  that  a  higher  degree  of  adjustment  accuracy  can  be  realized. 
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1 .2  Problem  Statement 


Until  recently,  the  mathematical  character  of  the  vibration/adjustment  solution 
space  had  been  assumed  to  be  linear;  however,  this  has  not  been  verified  by  a  rigorous 
study.  This  thesis  shall  examine  the  validity  of  the  linear  assumption  and  shall  detennine 
whether  higher  order  algorithms  are  necessary  to  accurately  calculate  main  rotor 
adjustments. 

1.3  Objectives 

The  objective  of  this  research  was  to  characterize  the  vibration/adjustment 
solution  space  as  it  is  understood  by  a  trained  neural  network,  and  to  search  for 
differences  in  adjustment  solutions  produced  by  a  linear  algorithm. 

1.4  Research  Methodology 

The  objectives  of  this  research  were  met  using  two  general  techniques.  The  first 
technique  was  essentially  a  detailed  comparison  of  the  calculated  adjustments  from  a 
neural  network  algorithm  to  those  of  a  traditional,  linear  algorithm.  If  the  two  algorithms 
behaved  the  same,  then  the  conclusion  would  be  that  the  neural  network  behaved  in  a 
linear  fashion.  If  they  produced  different  results,  however,  the  conclusion  would  not  be 
as  clear  cut. 

The  second  technique  involved  using  the  neural  network  to  create  linear 
coefficient  matrices  at  multiple  locations  in  the  vibration  and  adjustment  spaces.  These 
matrices  then  became  the  core  of  new  linear  algorithms  that  were  used  to  calculate 
adjustments  for  further  comparisons. 
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1.5  Chapter  Summary 

Main  Rotor  Smoothing  is  an  expensive  procedure  that  must  be  performed 
periodically  on  all  helicopters.  Due  to  the  high  cost  of  fuel,  maintenance  time,  and 
aircraft  unavailability,  improvements  are  sought  to  reduce  the  number  of  flights  required 
to  complete  a  MRS  procedure.  At  the  forefront  of  this  effort  is  the  creation  of  robust 
algorithms  that  calculate  more  accurate  adjustment  sets  than  today’s  linear  algorithms. 
This  study  will  quantify  the  major  differences,  if  any,  between  a  modern,  non-parametric 
algorithm  and  traditional,  linear  algorithms. 
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II.  Background 


2.1  Introduction 

Vibration  is  an  inherent  artifact  of  all  machinery.  Typically,  in  rotating 
machinery  vibrations  are  often  caused  by  a  machine’s  rotating  components  having 
centers  of  mass  that  do  not  he  at  the  center  of  rotation.  For  rotating,  aerodynamic 
components,  vibrations  can  also  be  caused  by  non-unifonn  aerodynamic  properties  such 
as  blade  pitch  or  camber.  In  helicopters,  vibrations  can  be  quite  powerful  and  lead  to 
shortened  component  lives  as  well  as  crew  fatigue.  With  non-aerodynamic  components 
vibrations  can  be  minimized  using  precise  manufacturing  techniques  and  mass  balancing 
procedures.  The  rotating,  aerodynamic  components  of  the  main  rotor  and  tail  rotor 
systems  require  that  specialized  vibration  reduction  procedures  be  performed.  This 
research  is  focused  on  the  largest  contributor  to  aircraft  vibrations:  the  main  rotor  system. 

2.2  General  Classification  of  Vibration 

Vibrations  are  usually  defined  by  the  oscillatory  motion  of  an  object  as  given  by 
equation  (1)  [4]. 

x(t)  =  A-  Cos{cot  - (j))  (1) 

In  this  equation,  A  is  the  vibration  amplitude,  or  maximum  displacement  of  the  object 
from  its  mean  position,  co  is  the  frequency  of  the  oscillation  and  4>  is  the  phase  shift  of  the 
steady  state  response  due  to  damping.  In  most  MRS  applications,  the  vibration 
displacement  x  is  measured  with  one  or  more  accelerometers  mounted  in  specific 
locations  on  the  aircraft. 
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The  frequency  of  the  vibration  that  MRS  is  concerned  with  is  the  same  as  that  of  main 
rotor  system  rotation.  The  phase  lag  is  a  product  of  the  structural  makeup  of  the  aircraft 
and  is  dependent  on  the  location  from  which  vibrations  are  measured.  Therefore, 
accelerometers  must  be  mounted  in  the  same  location  each  time  vibration  data  is 
recorded. 

2.3  Classification  of  main  rotor  vibrations 

The  main  rotor  system  produces  vibrations  that  are  quantified  in  a  number  of 
different  ways.  Typically,  vibration  is  described  as  being  a  vertical  vibration  or  a  lateral 
vibration  with  respect  to  the  helicopter  frame  of  reference.  The  amplitude  of  the 
vibration  is  typically  measured  as  velocity,  in  units  of  inches-per-second,  or  IPS.  While 
equation  (1)  requires  amplitude  in  terms  of  distance,  here  we  are  using  amplitude  to  mean 
the  total  distance  traveled  by  the  accelerometer  in  one  second.  The  tenns  amplitude  and 
magnitude  are  often  used  interchangeably  when  referring  to  the  strength  of  the  vibration. 
The  frequency  of  the  vibration  is  typically  a  multiple  of  the  main  rotor’s  period.  A  1/rev 
(pronounced  1-per-rev)  vibration  produces  one  vibration  cycle  for  each  revolution  of  the 
main  rotor  system.  A  1/rev  vibration  is  typically  brought  about  by  mass  and/or 
aerodynamic  asymmetries  in  the  main  rotor  system.  An  nl rev  vibration  produces  n  cycles 
per  revolution,  where  n  is  the  number  of  blades  on  the  main  rotor  system.  While  both 
types  of  vibrations  are  of  great  importance  to  helicopter  users  and  designers,  only  the 
1/rev  vibrations  shall  be  considered  in  this  study. 
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In  order  to  further  categorize  a  vibration,  a  quantity  known  as  phase  angle  is 
assigned  to  the  location  of  peak  amplitude.  Phase  angle  is  measured  from  an  aircraft 
specific  origin,  such  as  the  nose  or  tail  of  the  aircraft,  with  increasing  phase  opposite  the 
direction  of  main  rotor  rotation.  Test  equipment  is  used  to  record  vertical  and  lateral 
vibration  amplitude  as  well  as  the  main  rotor  angular  position.  The  vibration’s  peak 
amplitude  is  matched  with  angular  position  of  the  main  rotor  such  that  a  vibration  map 
(see  Figure  1)  can  be  plotted  with  respect  to  the  stationary  main  rotor.  This  is  particularly 
useful  in  determining  the  corrective  action  to  take  in  response  to  vibrations. 


figure  1.  Polar  Chart  of  AH-64D  vibrations  from  PC-GBS.  The  arrows  indicate 
predicted  vibrations  following  corrective  adjustments. 
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2.4  Flight  Conditions  for  Vibration  Measurement 


Now  that  we  can  quantify  a  helicopter  vibration  in  terms  of  direction,  magnitude, 
and  phase,  we  must  assign  a  flight  condition  at  which  the  vibration  should  be  measured. 
Typically,  a  helicopter’s  vibration  characteristics  are  measured  at  several  specific  flight 
states  ranging  from  on  the  ground  with  flat  blade  pitch  at  100%  RPM  (FPG-100)  to 
maximum  level  speed  (or  close  to  it).  Table  1  shows  the  flight  states  at  which  vibration 
measurements  are  recorded  for  two  types  of  helicopter.  The  quantity  of  flight  states  is 
carefully  chosen  so  as  to  give  a  good  representation  of  an  aircraft’s  vibration  signature 
while  minimizing  the  duration  of  the  test  flight. 

Table  1 .  Vibration  Measurement  Flight  Conditions 


UH-60 

AH-64 

FPG-100 

FPG-100 

Hover 

Hover 

80  Knots 

60  Knots 

120  Knots 

80  Knots 

145  Knots 

100  Knots 

120  Knots 

140  Knots 

2.5  Types  of  Main  Rotor  Adjustments 

With  a  clear  understanding  of  what  defines  a  main  rotor  vibration,  it  is  now 
important  to  determine  what  form  the  corrective  actions  should  take.  In  determining  this 
we  examine  a  few  of  the  phenomena  that  could  be  responsible  for  main  rotor  vibrations. 
Like  any  rotating  component,  the  main  rotor  system  could  suffer  from  a  mass  imbalance. 
There  is  also  the  likelihood  that  the  rotor  blades  are  not  aerodynamically  identical 
resulting  in  non-uniform  lift  production  on  the  rotor  disk. 
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It  is  possible  that  blade  stiffness  varies  from  blade  to  blade  allowing  some  blades  to  bend 
and  twist  more  than  others  at  equivalent  aerodynamic  loads.  While  this  is  not  a  complete 
list  of  causes  for  vibrations,  it  lends  insight  into  the  methods  we  may  employ  to  help 
eliminate  them.  These  methods  shall  include  modifying  the  aerodynamic  qualities  of  the 
blades  as  well  as  modifying  their  mass. 

2.5. 1  Adjustments  for  Mass  Imbalance. 

To  correct  for  mass  imbalances,  weights  can  be  added  to  or  subtracted  from 
individual  rotor  blades.  Rotor  blades  may  be  balanced  prior  to  installation  using  devices 
such  as  the  Universal  Static  Balance  Fixture  shown  in  Figure  2.  This  piece  of  equipment 
allows  one  to  ensure  consistency  with  the  spanwise  center  of  mass  of  each  rotor  blade 
prior  to  installation.  This  balance  technique  is  an  important  step  in  the  overall  MRS 
process  as  it  greatly  reduces  the  time  spent  balancing  the  mass  of  the  main  rotor  system 
after  blade  installation. 


Figure  2.  Universal  Static  Balance  Fixture.  This  piece  of  equipment  is  utilized  to 
statically  balance  the  spanwise  moment  prior  to  installation  [5]. 
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The  only  weight  adjustment  method  being  considered  in  this  study  is  the  addition  of 
blade  weights  in  response  to  flight  test  data.  Figure  3  shows  the  rotor  hub  of  an  AH-64 
Apache  with  blade  weights  installed.  What  follows  is  a  brief  explanation  of  why  mass 
imbalances  can  occur  in  a  main  rotor  system. 


Figure  3.  AH-64  Apache  with  blade  weight  installed  [5] 


2.5.2  Origin  of  a  Mass  Imbalance. 

Mass  imbalances  can  originate  at  the  manufacturing  level,  evolve  over  time,  or 
result  from  specific  events.  At  the  manufacturing  level,  main  rotor  system  components 
are  fabricated  to  meet  specific  engineering  tolerances.  The  size  of  these  tolerances  must 
be  small  enough  to  satisfy  the  demands  of  the  engineer,  yet  large  enough  for  affordable 
mass  production. 
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If  large  enough,  which  they  usually  are,  these  tolerances  can  lead  to  mass  imbalances  in 
the  main  rotor  system,  thus  requiring  weight  adjustments. 

Another  cause  for  mass  imbalance  is  that  over  time,  blade  erosion  occurs  due  to 
countless  impacts  of  the  rotor  blades  with  sand,  dirt,  and  other  debris.  Blade  erosion 
causes  the  rotor  system  to  slowly  lose  weight.  If  one  of  the  eroded  blades  is  replaced 
with  a  new  blade,  a  mass  imbalance  is  likely  to  be  introduced.  Since  most  helicopter 
users  do  not  change  all  of  their  blades  at  the  same  time,  mass  imbalances  are  almost  sure 
to  be  introduced  every  time  a  blade  is  changed. 

One  final  avenue  for  blade  mass  imbalance  is  from  a  specific  event  such  as  a 
blade  repair.  A  blade  repair  usually  consists  of  applying  a  small  amount  of  patching 
material  to  the  blade  to  mend  a  hole  or  crack.  The  addition  of  this  material  to  the  rotor 
blade  not  only  increases  the  blade  mass  but  also  alters  its  center  of  moment,  thus  leading 
to  an  increased  number  of  MRS  iterations. 

2.5.3  Aerodynamic  Imbalances . 

Aerodynamic  imbalances  may  be  created  at  the  same  time  as  mass  imbalances. 
When  rotor  blades  are  fabricated,  slight  imperfections  in  blade  shape  cannot  be  avoided. 
Even  if  they  were  avoided,  the  overall  shape  of  the  blade  would  be  changed  by  blade 
erosion  during  normal  flight  operations.  These  aerodynamic  imbalances  can  be  dealt 
with  in  a  number  of  ways.  One  way  is  to  change  the  angle  of  attack  of  the  blade  by 
adjusting  the  length  of  the  pitch  control  rod  (PCR).  Another  is  to  change  the  camber  of 
the  blade  by  bending  trim  tabs  located  on  the  blade  trailing  edge  (see  Figure  4).  While 
PCR  adjustments  change  the  AOA  of  the  entire  blade,  trim  tab  adjustments  only  affects 
the  camber  of  the  blade  over  the  length  of  the  trim  tab. 
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Figure  4.  A  PCR  on  a  UH-60  (left)  and  a  Trim  Tab  on  an  OH-58  [5] 

2.6  Main  Rotor  Adjustment  Fidelity 


Now  that  we  have  discussed  adjustment  types,  it  is  time  to  understand  the 
limitations  in  fidelity  of  the  adjustment  sets.  Every  adjustment  type  has  a  minimum 
allowable  magnitude  that  is  set  by  the  manufacturer.  For  instance,  the  blade  weights  may 
come  in  one  ounce  bars  or  the  PCR  may  have  a  locking  turnbuckle  that  allows  for  no  less 
than  V*  turn  of  adjustment.  An  adjustment  set  with  a  higher  degree  of  precision  than  this 
basic  unit  will  be  rounded  off  to  the  nearest  basic  unit. 

There  is  also  the  fundamental  principle  that  a  calculated  adjustment  move  may 
have  a  phase  angle  anywhere  from  0-360°  but  can  only  be  applied  on  actual  blade 
locations  (0°,  90  °,  180  °,  270  °  in  4  bladed  rotor  systems).  When  an  adjustment  is 
calculated  at  a  location  between  two  blades,  the  adjustment  must  be  divided  between  the 
two  blades  to  produce  an  equivalent  result.  Most  of  the  time  a  perfect  division  of  the 
solution  is  impossible  and  more  round  off  error  is  incurred. 
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2.1  History  of  Track  and  Balance  [6,71 

Over  the  decades  since  the  introduction  of  the  helicopter,  vibration  reduction 
procedures  have  changed  greatly.  Prior  to  the  incorporation  of  onboard  vibration 
measuring  equipment,  rotor  tracking  was  the  limit  of  our  capabilities  for  vibration 
reduction.  Essentially,  rotor  tracking  is  an  effort  to  make  all  rotor  blades  fly  the  same 
path,  thus  forming  a  flat  rotor  disk  at  100%  RPM.  The  principle  is  that  the  rotor  blades 
are  aerodynamically  similar  if  they  fly  the  same  path.  By  incorporating  this  tracking 
procedure  with  a  static  blade  balancing  as  discussed  earlier  in  section  2.5.1,  some 
vibration  reduction  is  achieved.  Unfortunately,  early  blade  tracking  techniques  could 
only  be  accomplished  using  ground  equipment,  thus  preventing  measurements  of 
aerodynamic  differences  at  various  flight  speeds. 

Later,  onboard  vibration  recording  equipment  was  developed  and  used  in 
conjunction  with  onboard  optical  tracking  equipment  for  in  flight  tracking  and  balancing 
of  the  rotor  system.  Highly  skilled  technicians  analyzed  the  magnitude  and  phase  of 
vibrations  as  well  as  track  split  data  in  order  to  determine  appropriate  adjustment  sets  for 
blade  pitch,  trim  tabs,  and  blade  weights.  Again,  this  technique  required  a  great  deal  of 
time  and  skill  to  accomplish. 

Finally,  with  the  advent  of  high  speed  digital  computers,  diagnostic  equipment 
was  created  that  could  measure  vibrations  and  blade  track  as  well  as  calculate  adjustment 
sets.  This  new  equipment  not  only  sped  up  the  rotor  smoothing  process  but  also  reduced 
the  requisite  skill  level  of  the  operator.  For  the  past  20  years  very  little  has  changed  in 
the  overall  concept  of  a  main  rotor  smoothing. 
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Vibrations  and  rotor  track  are  still  measured  in  much  the  same  way,  albeit  more 
accurately,  and  adjustments  sets  still  consist  of  blade  weight,  pitch  link,  and  trim  tab 
corrections. 

2.8  Linear  Algorithms 

One  area  that  has  seen  great  efforts  for  improvement  in  recent  years  is  the 
computer  software  used  to  calculate  adjustment  sets.  There  have  been  efforts  to  improve 
this  software  not  only  through  a  more  sophisticated  user  interface,  but  also  through  the 
algorithm  that  converts  raw  vibration  and  track  data  into  useful  adjustment  sets.  Two 
types  of  algorithms  are  being  considered  in  this  study.  The  first  type  is  an  algorithm  in 
wide  use  today  that  is  built  on  the  same  linear  assumption  used  since  computer  automated 
track  and  balance  began.  The  second  is  a  state-of-the-art  neural  network  algorithm  that  is 
being  used  by  a  select  few  units  in  the  US  Army. 

2.8.1  The  Aviation  Vibration  Analyzer. 

The  Aviation  Vibration  Analyzer  (AVA)  is  a  helicopter  vibration  measurement 
system  that  utilizes  a  linear  algorithm  in  its  main  rotor  smoothing  software.  That  is,  we 
assume  that  a  linear  relationship  exists  between  an  adjustment  and  the  resulting  change  in 
vibration.  While  this  greatly  simplifies  our  calculations,  there  is  the  potential  that  this 
assumption  may  not  be  accurate  and  could  lead  to  poor  adjustment  sets  in  environments 
where  the  rotor  system  behaves  in  a  non-linear  fashion. 
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2.8.2  The  Linear  Algorithm  Defined. 

By  making  the  linear  assumption,  we  are  able  to  employ  a  relatively  simple 
algebraic  model  to  describe  the  cause  and  effect  of  adjustment  sets  on  helicopter 
vibrations. 

[c]nxm  ■  [ML,  =  [^Vib Lxi  (2) 

In  Equation  (2),  [C]  is  a  sensitivity  coefficient  matrix  that  is  specific  to  a  type,  model,  and 
series  of  helicopter,  [Adj]  is  an  adjustment  set,  and  [AVib]  is  the  change  in  vibrations 
brought  about  by  the  adjustment  set.  N  stands  for  the  number  of  specific  flight  regimes 
over  which  data  is  collected  (see  Table  1),  and  M  stands  for  the  number  of  adjustment 
types  (i.e.,  blade  weight,  pitch  link,  and  trim  tab).  The  concept  of  a  linear  algorithm  can 
be  better  understood  if  we  first  know  how  the  sensitivity  coefficient  matrix  is  determined. 
In  doing  so,  we  will  look  at  the  case  of  the  UH-60. 

2.8.3  Determining  the  Linear  Sensitivity  Coefficient  Matrix. 

A  set  of  sensitivity  coefficients  for  a  helicopter  is  determined  through  a  series  of 
test  flights.  First,  a  baseline  vibration  set  is  recorded  at  the  N  flight  conditions  before  any 
adjustments  have  been  made.  Next,  a  single  adjustment  move  is  applied  to  the  main  rotor 
system  (e.g.  adding  a  5  oz  weight  to  the  yellow  blade).  The  helicopter  is  then  flown 
again  at  the  N  flight  conditions  and  a  new  vibration  set  is  recorded.  After  landing,  the 
adjustment  is  removed  and  a  different  adjustment  type  is  applied.  This  process  is 
repeated  until  vibration  responses  for  all  adjustment  types  have  been  recorded.  When 
complete,  the  C  matrix  is  populated  in  accordance  with  equation  (3). 


2-11 
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Typically,  several  coefficient  matrices  are  combined  together  in  a  root-mean- 
square  sense  in  order  to  detennine  a  mean  coefficient  matrix.  Due  to  several  factors, 
such  as  electronic  noise,  weather,  and  differences  between  like  airframes,  the  average 
coefficient  matrix  is  only  about  20%  accurate,  and  in  many  cases,  worse  [3].  For  this 
reason,  it  is  very  likely  that  several  test  flights  will  be  required  to  adjust  the  main  rotor 
system  to  an  acceptable  vibration  level. 

The  resulting  mean  coefficient  matrix  used  in  MRS  procedures  is  obtained  by 
rearranging  equation  1  and  solving  for  the  adjustment  set.  This  is  not  easily  done 
without  the  aid  of  a  computer,  due  to  the  fact  that  the  linear  system  typically  contains 
more  equations  than  unknowns.  In  the  case  of  the  AH-64D,  there  are  five  unknowns  to 
solve  for  seven  flight  regimes.  By  incorporating  a  solution  algorithm  that  minimizes  the 
sum  of  the  squares  of  predicted  vibration  magnitudes,  the  best  overall  adjustment  set  can 
be  calculated. 

2.9  Application  of  Adjustment  Sets 

Usually  a  main  rotor  smoothing  algorithm  is  written  to  perform  a  majority  of  its 
calculations  in  the  complex  plane.  In  this  fashion,  a  set  of  vibrations  enters  the  algorithm 
and  is  first  converted  to  a  complex  number  using  Equation  (4). 

Vib  ix  =  Mag  x  [Cos (phase)  +  i  x  Sin( phase)]  (4) 
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This  complex  vibration  is  then  converted  into  a  complex  adjustment  set  as  previously 
described.  The  complex  adjustment  set,  referred  to  as  a  reduced  adjustment  set,  must  be 
implemented  as  real  adjustments  on  individual  blades.  This  is  then  called  the  detailed 
adjustment  set.  There  is  typically  more  than  one  detailed  adjustment  set  that  corresponds 
to  a  single  reduced  set.  For  instance,  on  four-bladed  helicopters,  a  positive  adjustment 
applied  to  one  blade  has  an  equivalent  effect  on  vibration  as  a  negative  adjustment 
applied  to  the  opposite  blade  (180°  phase  shift).  The  difference  between  applying  one 
versus  the  other  lies  in  the  effect  the  adjustment  has  on  the  track  split. 

Since  a  small  track  split  is  desired,  an  additional  algorithm  is  typically 
incorporated  in  MRS  software  to  convert  reduced  adjustments  to  detailed  adjustments 
while  minimizing  track  split.  Such  an  algorithm  would  incorporate  track  data  collected 
with  equipment  such  as  the  Universal  Tracking  Device  (UDT)  shown  in  Fig.  5.  The 
UTD  is  an  electro-optical  device  that  is  usually  mounted  on  the  exterior  of  the  aircraft 
during  MRS  operations  and  is  able  to  “see”  the  flight  path  of  individual  blades  as  they 
enter  its  field  of  view.  In  recent  years  it  has  been  shown  that  track  split  is  not  a 
contributor  to  main  rotor  vibrations  (post-ground  balance)  and  is  optimized  mainly  for 
aesthetics. 
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2.10  The  Linear  Assumption  Comes  into  Question 

About  a  decade  ago,  Taitel  et  al[3]  questioned  the  viability  of  the  linear 
assumption  in  main  rotor  smoothing  algorithms.  It  was  postulated  that  the  accuracy  of 
the  adjustment  set  might  be  compromised  under  certain  conditions  if  higher  order  main 
rotor  interactions  were  neglected.  This  could  then  result  in  MRS  iterations  that  require 
multiple  flights.  Since  multiple  flights  are  typically  required  to  perform  a  MRS,  it  was 
believed  that  the  performance  of  MRS  algorithms  could  be  improved  by  including  the 
higher  order  interactions.  To  this  end,  neural  network  architectures  were  applied  to  the 
main  rotor  smoothing  problem.  The  premise  behind  this  approach  was  to  train  neural 
networks  using  vibration  data  from  multiple  test  flights  such  that  the  true  mathematical 
relationship  between  1/rev  vibrations  and  main  rotor  adjustments  could  be  learned. 
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2.1 1  The  Artificial  Neural  Network  Algorithm 


In  the  late  1990’s,  Wroblewski,  et  al  [2]  designed  and  implemented  a  neural 
network  based  software  system  for  use  on  the  US  Army’s  helicopter  fleet.  The  neural 
networks  were  trained  with  data  from  approximately  thirty  test  flights  per  type  of 
helicopter.  One  of  the  drawbacks  of  neural  networks  is  that  they  are  unable  to  extrapolate 
the  character  of  the  solution  space  for  vibration  regimes  in  which  they  haven’t  been 
trained.  Since  thirty  flights  are  insufficient  to  fully  train  a  neural  network,  modeled  data 
(linear  and  other)  was  also  incorporated  into  the  training  process  in  order  to  fill  gaps. 

The  intent  was  that  as  more  flights  became  available  the  neural  network  training  could  be 
updated,  thus  allowing  the  vibration/adjustment  mappings  to  mature  over  time. 

A  neural  network  is  essentially  a  non-linear  algorithm  that  is  able  to  make 
predictions  based  on  a  history  of  observations.  Figure  6  is  a  schematic  of  a  neural 
network  undergoing  training.  What  follows  is  a  simplified  explanation  of  how  the 
network  is  trained  to  calculate  MRS  adjustments. 


Figure  6.  Schematic  of  Neural  Network  Undergoing  Training  [8] 
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The  adjustment  portion  of  the  neural  network  is  trained  by  first  recording  a 
baseline  vibration  set.  An  adjustment  is  then  applied  to  the  main  rotor  system  and  a 
second  vibration  set  is  recorded.  The  difference  vector  between  the  second  vibration  set 
and  the  baseline  vibration  set  is  passed  through  the  neural  network  as  the  “Input” 
vibration  vector  in  Fig.  6.  A  set  number  of  weighted  functions  operate  on  the  vibration 
vector  to  create  an  adjustment  vector,  or  “Output”.  The  target  value  that  the  output 
adjustment  set  is  compared  to  is  the  negative  of  the  adjustment  set  that  was  applied  prior 
to  recording  the  second  vibration  vector.  Based  on  the  comparison,  the  weights  within 
the  neural  network  are  adjusted  until  the  output  matches  the  target.  This  process  is 
repeated  with  a  variety  of  adjustment  sets  until  the  network  is  fully  trained. 

The  appealing  aspect  of  using  a  neural  network  for  MRS  procedures  is  that  neural 
networks  can  be  designed  to  recognize  higher  order  relationships  between  adjustments 
and  their  vibration  response.  For  this  reasons,  a  properly  trained  neural  network  is 
believed  to  be  capable  of  describing  the  relationship  between  vibrations  and  adjustments 
more  accurately  than  a  linear  algorithm. 

2.11.1  The  PC-Ground  Based  System. 

The  neural  network  MRS  algorithm  that  was  considered  in  this  study  was  that  of 
the  Vibration  Management  Enhancement  Program’s  PC-Ground  Based  System  (PC- 
GBS).  This  algorithm  was  designed  to  work  for  the  same  flight  regimes  as  AVA.  In  fact, 
vibration  data  collected  with  AVA  equipment  is  able  to  be  analyzed  by  the  PC-GBS 
software.  After  vibration  data  has  been  entered,  the  default  “Best  Overall  Solution”  is 
calculated  to  reduce  vibration  magnitudes  to  as  close  to  zero  as  possible  in  the  fewest 
number  of  moves. 
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A  secondary  consideration  for  the  PC-GBS  is  to  apply  the  moves  such  that  the  main  rotor 
track  split  is  also  reduced. 

2.11.2  PC-GBS  Solution  Options. 

The  user  may  choose  other  solution  options  such  as  “Resolve  to  Vibration 
Limits”,  “Limit  Solution  to  x  Number  of  Moves”,  or  create  a  “Manual  Solution”.  If  the 
user  chooses  the  option  of  “Resolve  to  Vibration  Limits”,  an  adjustment  solution  is 
calculated  that  will  reduce  vibration  magnitudes  to  within  minimum  acceptable  limits 
with  the  fewest  number  of  moves.  If  the  user  chooses  “Limit  Solution  to  x  Number  of 
Moves”,  the  number  of  adjustment  moves  will  not  exceed  (but  can  be  less  than)  the  x 
number  of  moves.  Finally,  the  experienced  user  may  create  a  solution  using  a  graphical 
user  interface  by  selecting  the  “Manual  Solution”  option.  In  all  cases,  PC-GBS  predicts 
the  vibration  levels  and  track  split  following  application  of  the  adjustment  set. 

2.11.3  The  Neural  Networks  of  PC-GBS. 

The  architecture  used  in  the  PC-GBS  algorithm  consists  of  four  neural  networks. 
These  are  the  Adjustment  Evaluation  Networks  (AEN),  the  Vibration  Prediction  Network 
(VPN),  the  Track  Optimization  Network  (TON),  and  the  Solution  Optimization  Expert 
(SOE)  [1],  A  vibration  set  passes  through  these  networks,  entering  as  a  complex 
vibration  vector  and  exiting  as  a  set  of  specific  adjustments  with  a  corresponding 
vibration  and  track  prediction.  Each  network  performs  a  specific  function  that  shall  be 
described  next. 

The  AEN  is  a  set  of  networks  that  provide  the  mapping  of  the  complex  vibration 
vector  to  the  corresponding  complex  adjustment  space. 
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Several  candidate  adjustment  sets  are  calculated  and  then  passed  into  the  VPN  where 
they  are  mapped  back  to  the  vibration  space  in  order  to  predict  a  vibration  response. 

Next,  the  candidate  adjustment  sets  are  passed  into  the  TON  where  they  are  converted 
from  reduced  adjustment  sets  into  detailed  adjustments  for  specific  blades.  This  net 
incorporates  track  data  from  the  UTD  when  converting  the  complex  adjustments  into  real 
adjustments.  Finally,  the  SOE  reviews  the  candidate  adjustment  sets  and  chooses  the  best 
one,  based  on  the  predicted  vibration  levels,  the  number  of  actual  corrective  moves,  and 
the  predicted  track  split.  This  final  selection  network  is  trained  through  an  analysis  of 
real-life  examples  where  the  optimal  solution  was  chosen  by  human  experts  [1]. 

2.11.4  PC-GBS  Performance  From  an  Earlier  Study. 

Wroblewski  et  al  [1]  discuss  the  initial  results  of  the  PC-GBS  algorithm  on  the 
AH-64  Apache  helicopter.  They  noted  that  2-4  flights  were  typically  required  to  smooth 
the  rotor.  It  was  not  feasible  to  run  parallel  flight  tests  with  AVA  adjustments,  so  it  was 
hard  to  say  which  algorithm  was  superior;  however,  the  adjustment  solutions  were 
compared  between  the  two  algorithms.  The  neural  network  consistently  offered  solutions 
with  fewer  numbers  of  moves  and  lower  predicted  vibrations  than  those  of  the  linear 
algorithm.  This  alone  represents  a  noteworthy  improvement  of  the  neural  network  over 
the  linear  algorithm. 
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2.12  Chapter  Summary 

The  problem  of  main  rotor  smoothing  has  long  been  an  expensive  challenge  in  the 
helicopter  industry.  Over  the  years,  many  improvements  have  been  made  to  reduce  the 
number  of  flights  required  to  smooth  1/rev  vibrations  to  within  acceptable  levels.  The 
most  recent  improvement,  the  artificial  neural  network,  has  shown  potential  in  reducing 
the  number  of  iterations  (flights)  involved  in  a  typical  smoothing  process.  While  the  core 
of  the  software  is  protected  by  proprietary  rights,  its  analytical  capabilities  can  be 
employed  to  characterize  the  main  rotor  solution  space  as  it  is  known  to  the  trained  neural 
networks. 
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III.  Methodology 


3.1  Introduction 

The  primary  goal  of  this  study  was  to  gain  insight  into  a  main  rotor  system’s 
mathematical  nature  by  extracting  vibration/adjustment  relationships  from  trained  neural 
networks.  By  comparing  a  large  population  of  flight  data  over  a  broad  spectrum  of 
conditions,  characteristics  of  the  main  rotor  solution  space,  as  learned  by  the  neural 
network,  may  be  identified. 

3.2  Origin  of  Research  Data 

Comparisons  were  made  for  three  different  aircraft;  the  AH-64A  Apache,  AH- 
64D  Longbow,  and  UH-60  Blackhawk.  Each  of  these  helicopters  has  a  distinct  set  of 
trained  neural  networks  contained  in  the  PC-GBS  software,  as  well  as  its  own  set  of 
linear  sensitivity  coefficients.  For  the  purposes  of  this  study,  the  coefficients  from  the 
Aviation  Vibration  Analyzer  system  were  used  (Appendix  A).  Flight  test  data  from 
multiple  flights  of  each  type  of  aircraft  was  compiled  in  several  Microsoft  Excel 
databases  (Appendices  B,  C,  D)  and  was  evaluated  using  Matlab.  The  flight  data  utilized 
in  this  research  was  downloaded  from  an  online  VMEP  database  compiled  by  Intelligent 
Automation  Corporation  (IAC)  and  was  evaluated  using  PC-GBS  version  3.0  Build  439 
Service  Pack  2.  For  each  flight,  an  Excel  database  entry  was  made  consisting  of 
vibration  vectors  (magnitude  and  phase),  predicted  vibration  vectors,  and  detailed 
adjustment  values.  The  data  were  recorded  directly  from  the  Vibration  tab  (Fig.  7)  and 
the  Rotor  Smoothing  Solution  tab  (Fig.  8)  of  the  PC-GBS. 
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Vibration  Adjustment  —  UH60  |  79  23280  |  05/14/2005  10:39:37  |  FLI... 


Vibration  Values  |  Rotor  Smoothing  Solution  ]  Vibration  Plot  |  Track  Trend 
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Quality 

State 
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FpgIOO 

A-B 

0.03 

0.05 
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0.10 
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Print... 


Close 


Figure  7.  Vibration  Values  Tab  of  PC-GBS.  Predicted  vibration  values  in  this  figure 
correspond  to  adjustment  values  of  Fig.  8. 


Figure  8.  Rotor  Smoothing  Solution  Tab  of  PC-GBS 
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3.3  Data  Acquisition  Techniques 

Before  proceeding  further,  it  is  important  to  understand  the  method  used  to  extract 
data  from  the  PC-GBS.  There  are  four  neural  networks  contained  within  the  PC-GBS. 
These  are  the  AEN,  the  VPN,  the  TON,  and  the  SOE.  Of  these  four  networks,  only  the 
AEN  and  the  VPN  have  learned  the  relationship  between  main  rotor  adjustments  and 
vibrations.  The  effects  of  the  other  two  networks  must  therefore  be  removed  from  the 
data  in  order  to  obtain  a  clear  picture  of  the  main  rotor  solution  space. 

3.3.1  Removing  the  Effects  of  the  Solution  Optimization  Expert. 

The  purpose  of  the  SOE  is  to  choose  from  an  array  of  possible  adjustment  sets  in 
accordance  with  training  it  received  from  human  experts.  The  intent  of  this  training  was 
to  teach  the  SOE  to  offer  a  solution  with  the  minimum  acceptable  number  of  moves. 
Figures  7  and  8  show  the  adjustment  set  and  predicted  vibration  vector  when  the 
maximum  number  of  moves  have  been  requested.  Note  that  the  adjustment  set  consists  of 
six  adjustment  moves.  Had  the  “Best  Overall  Solution”  been  used,  the  adjustment  set 
would  have  been  a  single  move  of  +  1  notch  on  the  yellow  pitch  link.  The  selection  of 
this  single  adjustment  move  was  determined  by  the  SOE.  This  is  obviously  a  superior 
adjustment  set  from  the  mechanic’s  point  of  view,  because  it  reduces  the  time  to 
implement  adjustments  as  well  as  reduces  the  possibility  of  making  a  mistake.  The 
tradeoff  for  this  simpler  solution  is  slightly  higher  (but  still  favorable)  predicted  vibration 
magnitudes  and  track  split.  While  the  SOE  is  an  incredibly  important  contributor  to  the 
PC-GBS ’s  performance,  for  the  purposes  of  this  study  it  only  serves  to  mask  our 
understanding  of  the  AEN. 
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The  SOE  was  effectively  turned  off  by  calculating  adjustment  sets  using  the 
option  of  Limit  Solution  to  X  Moves  where  X  was  set  to  the  maximum  number;  as  shown 
in  Fig. 8.  For  the  AH-64  and  UH-60,  Jwas  set  to  20  and  12  respectively.  This  number 
was  determined  for  the  UH-60  by  multiplying  the  three  adjustment  types  over  the  four 
blades,  and  for  the  AH-64  by  multiplying  the  five  adjustment  types  over  the  four  blades. 
By  allowing  the  maximum  number  of  moves,  PC-GBS  produces  the  most  highly  detailed 
adjustment  set  attainable  from  the  AEN.  Since  maximizing  the  number  of  moves  is  the 
antithesis  of  the  SOE,  it  is  believed  that  this  procedure  effectively  disables  it. 

3.3.2  Removing  the  Effects  of  the  TON. 

Unlike  the  AEN,  the  TON  is  concerned  primarily  with  minimizing  the  main  rotor 
track  split.  When  a  reduced  adjustment  set  is  calculated  by  the  AEN,  the  adjustment 
magnitudes  are  all  positive  numbers  with  associated  phase  angles.  One  could  apply  this 
reduced  adjustment  set  to  the  main  rotor  system  as  given,  thus  resulting  in  a  beneficial 
change  in  the  main  rotor  vibration  vector.  One  could  just  as  well  apply  the  negative  of 
the  adjustment  magnitude  at  a  180  degree  phase  shift  and  achieve  the  exact  same  change 
in  the  vibration  vector.  The  difference  between  applying  the  positive  versus  the  negative 
adjustment  magnitude  lies  in  the  effect  it  has  on  the  main  rotor  track.  One  of  the 
purposes  of  the  Track  Optimization  Network  is  to  detennine  which  adjustment  moves  to 
keep  as  positive  magnitudes  and  which  to  make  negative  magnitudes  in  order  to 
minimize  the  main  rotor  track  split.  The  TON  makes  this  determination  as  it  converts  a 
reduced  adjustment  set  to  a  detailed  adjustment  set  on  actual  rotor  blades. 
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The  method  that  was  used  in  this  study  to  convert  the  detailed  PC-GBS 
adjustment  sets  into  reduced  adjustment  sets  automatically  converted  all  magnitudes  back 
to  positive  values.  The  next  section  provides  a  thorough  explanation  of  this  conversion 
process. 

3.3.3  Converting  Detailed  Adjustment  Sets  into  Reduced  Adjustment  Vectors. 

As  discussed  earlier,  vibration  vectors  were  converted  from  real  vectors  to 
complex  vectors  in  accordance  with  equation  (4),  repeated  here  for  convenience. 


Vibcmptx  ~  Mctg  x  [Cos (phase)  +  i  x  Sin(  phase)] 


(4) 


The  detailed  adjustment  sets  were  converted  to  complex  numbers  in  much  the  same  way. 
In  the  case  of  the  UH-60,  equation  (5)  was  used. 


Adj cmpix  ~  AdjMag  x  \Cos(  Blade  Phase-  90)  -  i  x  Sin(BladePhase-  90)]  (5) 


Here,  the  value  for  BladePhase  can  only  take  on  values  of  0°,  90  °,  180  °,  or  270  °,  as 
these  are  the  respective  phase  angles  associated  with  the  black,  yellow,  blue,  and  red 
blades.  The  value  for  AdjMag  is  the  magnitude  of  the  adjustment  move  on  the  blade 
being  considered.  Figure  9  shows  the  top  down  view  of  a  UH-60  Blackhawk  with 
appropriate  phase  values  assigned  to  each  blade. 
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270° 
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Figure  9.  UH-60  Blackhawk  Blade  Map. 
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For  the  UH-60,  the  Matlab  code  in  Appendix  E  converts  the  magnitude  of  the 
detailed  adjustment  on  the  given  blade  to  a  positive  or  negative,  real  or  imaginary  number 
in  accordance  with  Figure  9  and  Equation  5.  The  complex  numbers  for  each  adjustment 
type  (weight,  pitch  link,  trim  tab)  were  then  summed,  thus  producing  a  single  complex 
number  for  each  adjustment  type.  This  complex  adjustment  vector  was  then  converted 
back  into  real  values  of  magnitude  and  phase.  Table  2  contains  sample  calculations  for 
two  different  sets  of  pitch  link  adjustments.  These  two  sets  of  adjustments  represent  the 
same  reduced  adjustment  vector,  thus  emphasizing  the  fact  that  the  effects  of  the  TON 
can  be  removed  from  the  analysis.  A  similar  mapping  can  be  demonstrated  for  the  AH- 
64 A/D  adjustments. 


Table  2.  Sample  Calculation.  Conversion  of  UH-60  Pitch 
Link  adjustments  into  a  Reduced  Adjustment  Vector. 


Adjustment  Set  1 

Blade  Color 

Yellow 

Blue 

Red 

Black 

Pitch  Move 

2 

1 

1 

2 

Complex  Value 

2 

-i 

-1 

2i 

Complex  Vector 

] 

1  +i 

Real  Vector 

Mag:  1.4142 

Phase:  45  degrees 

Adjustment  Set  2 

Blade  Color 

Yellow 

Blue 

Red 

Black 

Pitch  Move 

-1 

-2 

-2 

-1 

Complex  Value 

-1 

2i 

2 

-i 

Complex  Vector 

1 

1  +i 

Real  Vector 

Mag:  1.4142 

Phase:  45  degrees 

3.4  Overview  of  Analysis  Methods 

With  the  effects  of  the  TON  and  SOE  effectively  removed  from  the  data,  the  full 
capabilities  of  the  AEN  and  VPN  could  be  analyzed.  Two  methods  were  used  to  extract 
information  from  the  Adjustment  Evaluation  Networks  and  the  Vibration  Prediction 
Network. 


3-6 


The  first  method  used  a  graphical  and  statistical  comparison  of  the  reduced  adjustment 
sets  of  PC-GBS  to  the  reduced  adjustment  sets  of  a  linear  algorithm  based  on  the  AVA 
sensitivity  coefficients.  The  second  method  also  relied  heavily  on  graphical  and 
statistical  comparisons;  however,  the  comparisons  were  made  with  multiple  sets  of  ad 
hoc  linear  sensitivity  coefficients  derived  from  the  PC-GBS.  Both  methods  were  applied 
to  databases  representative  of  the  full  spectrum  of  vibration  regimes. 


3.5  Explanation  of  Vibration  Ranges 

The  comparison  of  the  neural  network  algorithm  and  the  linear  algorithm  was 
carried  out  for  multiple  flights  in  each  of  the  four  vibration  categories,  Good,  Above 
Goal,  Caution,  and  Exceed.  A  vibration  category  is  assigned  to  each  flight  based  on  the 
highest  vibration  magnitude  encountered  during  the  flight.  Table  3  shows  the  vibration 
magnitude  ranges  associated  with  each  category  by  type  of  aircraft.  By  analyzing  the  full 
spectrum  of  vibration  categories,  non-linear  effects  associated  with  vibration  magnitude 
should  be  revealed. 


Table  3.  Vibration  Magnitude  Ranges  by  Category 


Category 

UH-60 

AH -64  A/D 

Good  Vert  (Lat) 

0.0-0.25  (0.0 -0.2) 

0.0 -0.3  (0.0 -0.2) 

Above  Vert  (Lat) 

0.25-0.5(0.2-0.5) 

0.3 -0.5  (0.2 -0.5) 

Caution  Vert  (Lat) 

0.5 -0.8  (0.5 -0.8) 

0.5 -0.8  (0.5 -0.8) 

Exceed  Vert  (Lat) 

0.8+  (0.8+) 

0.8+  (0.8+) 

3.6  Calculating  the  AVA  Adjustment  Set 

The  Matlab  codes  of  Appendices  E,  F,  and  G  were  used  to  calculate  reduced  AVA 
adjustments  for  the  UH-60,  AH-64A,  and  AH-64D  respectively,  using  equation  (6). 

[^L,=[cl„4a^L,  (6) 
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In  this  equation,  the  adjustment  set  produces  the  vibration  response,  AVib,  when  applied 
to  the  main  rotor  system.  Since  we  wish  to  zero  out  our  measured  vibrations,  the  AVib  in 
equation  (6)  is  set  to  the  negative  value  of  our  measured  vibration  vector. 

The  C  matrix  is  comprised  of  the  AVA  sensitivity  coefficients  from  Appendix  A. 
Since  Equation  (6)  is  overdetermined  (more  equations  than  unknown),  the  Matlab  Left 
Divide  operation,  \,  essentially  solves  Equation  (7)  where  s  is  a  vector  that,  when  added 
to  AVib,  satisfies  the  equality. 

[^^L,=[CL„H[AF»]m+[£]„,)  (7) 

The  value  for  s  is  calculated  such  that  the  sum  of  the  squares  of  its  magnitudes  is 
minimized.  This  value  for  s  may  be  viewed  as  a  predicted  vibration  vector  following 
application  of  the  Adj  vector;  however,  it  corresponds  to  a  reduced  adjustment  set  that  has 
not  been  distributed  to  actual  rotor  blades.  Due  to  round  off  error  associated  with 
mechanical  limitations,  i.e.  weight  adjustments  in  whole  ounce  increments,  the  predicted 
vibration  vector  will  change  after  a  detailed  adjustment  set  has  been  determined. 

3.7  The  Ad  Hoc  Sensitivity  Coefficient  Method 

The  second  method  used  in  this  study  offered  a  good  deal  more  insight  into  the 
neural  network’s  mathematical  nature  than  the  first  method.  In  the  second  method,  an  ad 
hoc  sensitivity  coefficient  matrix  was  created  for  each  flight  by  using  the  PC-GBS  as  an 
adjustment/vibration  simulator.  This  method  not  only  allowed  for  a  straightforward 
comparison  of  the  ad  hoc  coefficients  to  those  of  AVA,  but  also  to  ad  hoc  coefficients  of 
other  flights.  In  this  manner,  specific  flight  regimes  could  be  studied  for  signs  of  non¬ 
linearity. 
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3.7.1  Constructing  the  Ad-Hoc  Sensitivity  Coefficient  Matrix. 

The  procedures  for  determining  a  set  of  linear  sensitivity  coefficients  were 
discussed  in  detail  in  section  2.8.3.  These  same  procedures  were  followed  in  order  to 
determine  the  ad  hoc  sensitivity  coefficients  from  PC-GBS.  Equation  (3),  repeated  here 
for  convenience,  shows  the  mathematical  model  used  to  determining  each  coefficient. 


c  = 

n ,  m 


fKfierj,  ~  ViKfore, 

Adj,„ 


(3) 


To  determine  the  ad  hoc  coefficients,  the  PC-GBS  was  used  to  predict  the  vibration 
vector,  Vib after  n,  following  the  application  of  a  single  adjustment,  Adjm,  to  a  flight  with 
measured  vibration,  Vibteforen •  This  process  was  repeated  for  each  adjustment  type  in 
order  to  populate  the  ad  hoc  sensitivity  matrix  for  each  flight. 

3.8  Database  Entries 

Table  4  shows  an  example  of  a  database  entry  for  a  single  UH-60  flight.  The 
AVA  comparison  method  uses  the  measured  vibrations,  the  detailed  adjustment  set,  and 
their  corresponding  predicted  vibrations.  The  ad  hoc  coefficient  method  uses  additional 
vibration  predictions  based  on  application  of  single  adjustment  moves  to  single  blades 
using  the  Manual  Solution  option.  In  the  case  of  Table  4,  the  smallest  allowable 
adjustment  moves  were  used  with  PC-GBS  to  produce  the  ad  hoc  predictions. 
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For  each  flight,  a  second  database  was  created  using  the  largest  allowable 
adjustment  moves.  Table  5  indicates  the  smallest  and  largest  allowable  moves  for  each 
aircraft  considered  in  this  study.  The  reasoning  behind  comparing  ad  hoc  coefficient 
matrices  based  on  small  versus  large  adjustment  moves  was  to  discover  non-linearity  as  a 
function  of  adjustment  magnitude. 


Table  4.  Typical  database  entry  for  the  UFI-60  Blackhawk.  Magnitude  is  in  IPS.  Weight 
Wt)  is  in  ounces,  Pitch  Links  (P/L)  is  in  Flats,  and  Tab  is  in  Mils.  Data  entered  by  hand. 


Measured  Vib  Detailed  Adj  Ad  Hoc  Prediction  Ad  Hoc  Prediction  Ad  Hoc  Prediction 

Batch 

1 

Before 

Mag 

Before 

Phase 

Pred 

Mag 

Pred 

Phase 

Yellow 

Wt 

After 

Mag 

After 

Phase 

Yellow 

P/L 

After 

Mag 

After 

Phase 

Yellow 

Tab 

After 

Mag 

After 

Phase 

FpgIOO  A-B 

0.032707 

247 

0.0197 

85 

5 

0.069268 

307 

1 

0.051184 

152 

2 

0.041347 

198 

Hover  A-B 

0.083903 

82 

0.01359 

322 

0.065411 

85 

0.1288 

108 

0.096912 

91 

80  Kt  A-B 

0.045817 

167 

0.047175 

220 

0.048443 

188 

0.067723 

147 

0.051565 

155 

120Kt  A-B 

0.121328 

70 

0.072691 

59 

0.103722 

71 

0.149222 

85 

0.132964 

76 

145Kt  A-B 

0.064676 

79 

0.02042 

143 

0.047172 

84 

0.100951 

102 

0.075388 

91 

Hover  A+B 

0.164067 

274 

0.156546 

268 

0.164286 

274 

0.177728 

274 

0.172191 

272 

80Kt  A+B 

0.117083 

323 

0.074428 

285 

0.116287 

323 

0.091119 

303 

0.100755 

307 

120Kt  A+B 

0.149976 

349 

0.008687 

244 

0.147891 

349 

0.080251 

342 

0.096581 

336 

145Kt  A+B 

0.181986 

1 

0.087178 

108 

0.179131 

1 

0.097419 

360 

0.097789 

354 

Yellow 

Blue 

Red 

Black 

Wt 

P/L 

Tab 

Wt 

P/L 

Tab 

Wt 

P/L 

Tab 

Wt 

P/L 
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Detailed 

7 

0 

0 

-6 

-2 

0 

0 

1 

-6 

0 

0 

-7 

T able  5.  Adjustment  Move  Sizes  used  to  detennine  the  ad  hoc  coefficients . 


UH-60 

UH-60 

AH-64 

AH-64 

Weight 

5  oz 

80  oz 

113  grams 

1017  grams 

Pitch  Link 

1  Notch 

30  Notches 

0.5  flats 

12  flats 

Trim  Tab 

2  Mils 

20  Mils 

0.5  degrees 

5  degrees 

3.9  Chapter  Summary 

The  purpose  of  this  study  was  to  search  for  non-linear  relationships,  as  learned  by 
a  trained  neural  network,  between  main  rotor  adjustments  and  their  resulting  change  in 
vibrations.  This  comparison  was  made,  in  part,  with  the  use  of  a  linear  algorithm  based 
on  the  AVA  sensitivity  coefficients  of  each  aircraft. 
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Comparisons  were  made  over  a  broad  spectrum  of  measured  vibration  amplitudes  in  an 
effort  to  discover  whether  the  shape  of  the  adjustment  mapping  changes  due  to  the 
roughness  or  smoothness  of  the  test  flight.  The  ad  hoc  coefficient  matrices  were 
constructed  using  both  small  and  large  adjustment  moves  in  an  effort  to  discover  non¬ 
linear  effects  based  on  adjustment  amplitude.  By  using  these  two  approaches,  the  two 
avenues  for  non-linear  effects  have  been  rigorously  explored  and  any  non-linear  effects 
should  reveal  themselves  in  the  graphical  and  statistical  analysis. 
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IV.  Analysis  and  Results 


4.1  Introduction 

The  Blackhawk,  Apache,  and  Longbow  were  studied  over  a  broad  spectrum  of 
vibration  and  adjustment  magnitudes.  A  solid  picture  has  been  developed  as  to  how  the 
PC-GBS  converts  measured  vibration  vectors  into  adjustment  vectors.  Using  graphical 
and  statistical  analysis,  the  objectives  of  this  thesis  have  been  met. 

4.2  Comparison  of  PC-GBS  Adjustments  to  AVA  Adjustments 

A  linear  algorithm  was  created  for  each  aircraft  based  on  the  sensitivity 
coefficients  of  the  US  Anny’s  Aviation  Vibration  Analyzer.  Figure  10  is  a  bar  chart  of 
the  20  UH-60  flights  included  in  this  study.  The  height  of  each  bar  represents  the 
magnitude  of  the  difference  between  reduced  adjustment  vectors  as  determined  by  the 
AVA  algorithm  and  the  PC-GBS  algorithm.  Differences  that  do  not  exceed  the  basic 
adjustment  unit  are  considered  to  represent  identical  adjustment  calculations. 
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Pitch  Link 


Trim  Tab 


Figure  10.  UH-60  Bar  Chart  Comparison  of  AVA  to  PC-GBS.  This  figure  s' 
difference  between  PC-GBS  adjustments  and  AVA  adjustments. 


lows  the 


Only  the  Pitch  Link  and  Trim  Tab  adjustment  differences  are  depicted  in  Figure  10  as 
AVA  does  not  produce  weight  adjustments  based  on  flight  data  for  the  UH-60.  It  is 
obvious  that  the  amount  of  difference  between  the  two  adjustment  sets  is  significant. 
These  differences  were  calculated  by  measuring  the  length  of  the  vector  separating  an 
individual  AVA  solution  from  its  counterpart  PC-GBS  solution.  These  quantities  are 
graphically  depicted  in  Figure  11. 
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UH-60  Pitch  Link  Adjustment  (notches) 


Figure  11.  Graphical  Depiction  of  Vector  Difference. 
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Figure  12  shows  the  magnitude  differences  between  the  AH-64A  adjustment  sets. 
Again,  for  the  twenty  flights  considered,  the  differences  between  the  neural  network 
solution  and  the  AVA  solution  are  often  large.  This  does  not  imply  that  a  non-linear 
mapping  exists  in  the  PC-GBS.  It  only  means  that  AVA  and  PC-GBS  produce  very 
different  adjustments  for  the  same  vibration  vector. 


Weight  Pitch  Tab  8-10  Tab  6-10  Tab  4-10 
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Figure  12.  AH-64A  Bar  Chart  Comparison  of  AVA  to  PC-GBS.  This  figure  shows  the 
difference  between  PC-GBS  adjustments  and  AVA  adjustments. 
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Figure  13  shows  the  magnitude  differences  between  the  AF1-64D  adjustment  sets. 
Unlike  the  previous  comparisons,  here  the  adjustment  sets  are  identical.  This  indicates 
that  the  PC-GBS  has  learned  a  linear  mapping  for  the  AH-64D  that  is  a  near  perfect 
match  to  the  AVA  coefficients.  This  simple  analysis  has  helped  to  show  the  similarity  of 
the  two  algorithms  for  one  aircraft  in  this  study;  however,  a  different  method  was  used  to 
better  understand  the  other  two. 


Figure  13.  AFI-64D  Bar  Chart  Comparison  of  AVA  to  PC-GBS.  This  figure  shows  the 
difference  between  PC-GBS  adjustments  and  AVA  adjustments. 

4.3  Ad  Floe  Coefficient  Method 

In  order  to  draw  conclusions  about  the  characteristics  of  the  adjustment  space  for 
the  UH-60  and  the  AH-64A,  sets  of  ad  hoc  sensitivity  coefficients  were  developed  for 
every  flight.  These  ad  hoc  coefficients  were  analogous  to  the  AVA  coefficients  of 
Appendix  A.  These  new  coefficients  were  then  used  to  calculate  adjustment  sets  parallel 
to  those  of  the  PC-GBS  for  the  20  flights  analyzed. 
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4.3.1  Small  Moves  Ad  Hoc  Coefficient  Analysis. 

Figure  14  is  the  polar  plot  of  the  UFI-60  weight  adjustment  sets  as  determined  by 
the  ad  hoc  coefficients  and  the  PC-GBS  for  all  20  flights.  The  ad  hoc  coefficients  used  in 
this  set  of  analyses  were  developed  using  the  smallest  adjustment  moves  allowable. 


Weight 

90 


50 


Figure  14.  UFI-60  Small  Moves  Ad  Hoc  Weight  Adjustments.  This  chart  represents  20 
flights  as  calculated  by  the  PC-GBS  (+)  and  the  small  moves  ad  hoc  sensitivity 
coefficients  (O). 


For  the  majority  of  flights,  the  +  and  O  markings  are  on  top  of  one  another.  This 
indicates  that  the  adjustment  sets  as  calculated  by  the  ad  hoc  coefficients  and  the  AEN 
are  nearly  identical.  The  largest  magnitude  difference  between  the  PC-GBS  weight 
adjustments  and  those  of  the  ad  hoc  method  is  2.178  oz. 
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Similar  polar  plots  could  be  generated  for  the  AH-64A  and  AH-64D;  however,  a  more 
concise  technique  for  observing  similarities  in  parallel  adjustment  sets  is  to  tabulate  the 
largest  magnitude  of  dissimilarity  for  each  adjustment  type  (weight,  pitch  link,  tab). 

For  every  flight,  a  difference  vector  separating  parallel  adjustments  was 

determined,  as  was  shown  in  Fig  11.  Table  6  contains  the  magnitude  of  the  largest 

difference  vector  per  adjustment  type.  As  a  reference,  the  table  also  contains  the 

minimum  adjustment  unit  that  is  mechanically  allowable  on  the  UH-60.  Tables  7  and  8 

contain  similar  information  for  the  AH-64A  and  AH-64D. 

Table  6.  UH-60  Small  Moves  Ad  Hoc  Comparison.  This  table  shows  the  largest 
difference  between  small  moves  ad  hoc  adjustment  sets  and  PC-GBS  adjustment  sets. 


UH-60 

iwarni 

Pitch  Link 

Tab 

Largest  Adjustment  Difference 

2.178 

1.253 

4.699 

Basic  Adjustment  Unit 

1  oz 

1  Notch 

1  Mil 

Table  7.  AH-64A  Small  Moves  Ad  Hoc  Comparison.  This  table  shows  the  largest 


difference  between  small  moves  ad  hoc  adjustment  sets  and  PC-GBS  adjustment  sets. 


AH-64A 

Weight 

Pitch  Link 

Tab  8-10 

Tab  6-10 

Tab  4-10 

Largest  Adj  Difference 

43.78 

0.417 

0.960 

0.838 

0.471 

Basic  Adjustment  Unit 

0.25  Flats 

BifcliHIl 

ESI 

0.5  deg 

Table  8.  AH-64D  Small  Moves  Ad  Hoc  Comparison.  This  table  shows  the  largest 


difference  between  small  moves  ad  hoc  adjustment  sets  and  PC-GBS  adjustment  sets. 


AH-64D 

Pitch  Link 

Tab  8-10 

Tab  6-10 

Tab  4-10 

Largest  Adj  Difference 

96.37 

0.492 

2.063 

1.248 

0.437 

Basic  Adjustment  Unit 

52  grams 

0.25  Flats 

0.5  deg 

0.5  deg 

0.5  deg 

For  all  three  aircraft,  it  is  seen  that  the  largest  adjustment  difference  is  usually  greater 
than  the  basic  adjustment  unit.  This  indicates  that  the  PC-GBS  may,  under  certain 
circumstances,  produce  slightly  different  adjustment  sets  than  the  ad  hoc  linear 
coefficient  algorithm. 
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Since  the  AH-64D  has  already  been  shown  to  behave  like  AVA,  yet  still  has  magnitude 
differences  exceeding  the  minimum  adjustment  unit,  the  values  noted  in  these  tables 
should  not  be  interpreted  as  proof  of  non-linearity  in  the  AEN  or  VPN  of  PC-GBS. 

4.3.2  Method  of  Root  Mean  Squares  of  the  Small  Moves  Ad  Hoc  Coefficients . 

During  the  course  of  this  analysis,  it  was  observed  that  the  ad  hoc  coefficient 

matrices  were  virtually  identical  from  flight  to  flight.  For  each  type  aircraft,  the  ad  hoc 

coefficient  matrices  were  summed  into  a  single  root-mean-square  (RMS)  coefficient 

matrix.  All  flights  for  each  aircraft  were  then  re-evaluated  using  the  respective  single 

RMS  matrix.  Tables  9,  10,  and  1 1  are  analogous  to  Tables  6-8  in  that  they  show  the 

maximum  difference  encountered  between  the  PC-GBS  solution  and  the  ad  hoc  RMS 

solution.  Also  contained  in  these  tables  are  the  largest  differences  between  individual  ad 

hoc  coefficient  adjustments  and  PC-GBS  adjustments  from  tables  6-8. 

Table  9.  UH-60  Comparison  for  RMS  of  Small  Moves  Ad  Hoc.  This  Table  shows  the 
largest  difference  between  RMS  of  small  moves  ad  hoc  and  PC-GBS  adjustment  sets. 


UH-60  Small  Moves 

WMStnzm 

Pitch  Link 

Tab 

Largest  Individual  Difference 

2.178 

1.253 

4.699 

Largest  RMS  Difference 

1.402 

0.475 

1.166 

Basic  Adjustment  Unit 

1  oz 

1  Notch 

1  Mil 

Table  10.  AH-64A  Comparison  for  RMS  of  Small  Moves  Ad  Hoc.  This  Table  shows  the 
argest  difference  between  RMS  of  small  moves  ad  hoc  and  PC-GBS  adjustment  sets. 


AH-64A  Small  Moves 

BWJigffiaa 

Pitch  Link 

Tab  8-10 

Tab  6-10 

Tab  4-10 

Largest  Individual  Difference 

43.78 

0.417 

0.960 

0.838 

0.471 

Largest  RMS  Difference 

33.901 

0.410 

0.426 

0.363 

0.337 

Basic  Adjustment  Unit 

52  grams 

0.25  Flats 

0.5  deg 

0.5  deg 

0.5  deg 

4-8 


Table  11.  AH-64D  Comparison  for  RMS  of  Small  Moves  Ad  Hoc.  This  Table  shows  the 


argest  difference  between  RMS  of  small  moves  ad  hoc  and  PC-GBS  adjustment  sets. 


AH-64D  Small  Moves 

Pitch  Link 

Tab  8-10 

Tab  6-10 

Tab  4-10 

Largest  Individual  Difference 

96.37 

0.492 

2.063 

1.248 

0.437 

Largest  RMS  Difference 

14.698 

0.066 

0.219 

0.125 

0.043 

Basic  Adjustment  Unit 

52  grams 

0.25  Flats 

0.5  deg 

0.5  deg 

0.5  deg 

The  most  important  thing  to  leam  from  studying  these  tables  is  that,  for  each  type 
of  helicopter,  a  single  set  of  linear  sensitivity  coefficients  is  able  to  produce  reduced 
adjustment  sets  that  are  virtually  identical  to  those  of  the  AEN  of  PC-GBS.  This  is  proof 
that  the  AEN  calculates  reduced  adjustment  sets  using  a  linear  mapping.  Table  12  offers 
the  standard  deviation  of  the  difference  between  PC-GBS  adjustments  and  those 
produced  with  the  RMS  of  the  large-move  ad  hoc  method.  The  values  in  this  table  are 
less  than  the  basic  adjustment  unit,  thus  proving  that  a  majority  of  ad  hoc  adjustments  are 
identical  to  those  of  the  PC-GBS.  Figure  15  is  a  similar  polar  plot  to  figure  14  except 
that  the  single  RMS  matrix  of  the  small  moves  ad  hoc  coefficients  was  used  to  calculate 
parallel  adjustments  to  those  of  PC-GBS. 

Table  12.  Standard  Deviation  of  Differences  of  RMS  of  Small  Moves  Ad  Hoc.  This 
table  shows  standard  deviation  of  the  difference  in  adjustments  as  determined  by  the  PC- 


GBS  and  the  RMS  of  the  small  move  ad  1 

loc  coefficients. 

Weight 

Pitch  Link 

Trim  Tab 

UH-60 

0.342  oz 

0.1  Notch 

0.36  Mil 

Weight 

Pitch  Link 

Tab  8-10 

Tab  6-10 

Tab  4-10 

AH-64A 

7.86  g 

0.10  Flats 

0.12  deg 

0.10  deg 

0.08  deg 

AH-64D 

3.43  g 

0.01  Flats 

0.04  deg 

0.03  deg 

0.01  deg 
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Figure  15.  Polar  Chart  of  UFI-60  Weight  Adjustments.  This  chart  represents  20  flights 
as  calculated  by  the  PC-GBS  (+)  and  the  single  RMS  of  the  small  moves  ad  hoc 
sensitivity  coefficients  (O). 

The  preceding  analysis  included  an  equal  distribution  of  flights  from  each  of  the 
four  vibration  categories.  The  results  of  that  analysis  showed  a  high  degree  of  linearity 
exists  in  the  AEN  and  VPN  networks  regardless  of  measured  vibration  magnitude.  The 
next  section  will  address  the  effects  of  high  magnitude  adjustments  on  the  VPN. 

4.3.3  Large  Moves  Ad  Hoc  Coefficient  Analysis. 

In  order  to  determine  whether  large  adjustment  magnitudes  produce  non-linear 
predictions  from  the  VPN,  a  second  set  of  ad  hoc  coefficients  was  created  for  each  flight 
using  the  maximum  allowable  adjustment  magnitude  of  Table  5. 
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It  was  believed  that  these  new  coefficients  would  produce  different  adjustment  sets  than 

the  previous  small  move  coefficients  if  any  non-linear  mappings  had  been  learned  by  the 

VPN.  As  tables  13-15  show,  this  is  not  the  case.  Here  are  tabulated  the  largest 

magnitude  of  the  difference  vector  between  adjustments  as  detennined  by  PC-GBS  and 

the  large-move-ad-hoc-coefficients.  Also  contained  in  these  tables  are  the  largest 

differences  between  individual  ad  hoc  coefficient  adjustments  and  PC-GBS  adjustments. 

Table  13.  UH-60  Large  Moves  Ad  Hoc  Comparison.  This  Table  shows  the  largest 
difference  between  large  moves  ad  hoc  and  PC-GBS  adjustment  sets. _ 


UH-60  Large  Moves  Ad  Hoc 

mwistsizm 

Pitch  Link 

Tab 

Largest  Individual  Difference 

1.069 

0.392 

0.394 

Largest  RMS  Difference 

0.779 

0.393 

0.523 

Basic  Adjustment  Unit 

1  oz 

1  Notch 

1  Mil 

Table  14.  AH-64A  Large  Moves  Ad  Hoc  Comparison.  This  Table  shows  the  largest 
difference  between  large  moves  ad  hoc  and  PC-GBS  adjustment  sets. _ 


AH-64A 

Weight 

Pitch  Link 

Tab  8-10 

Tab  6-10 

Tab  4-10 

Largest  Individual  Difference 

37.938 

0.672 

0.596 

0.333 

0.366 

Largest  RMS  Difference 

39.219 

0.647 

0.486 

0.360 

0.318 

Basic  Adjustment  Unit 

52  grams 

0.25  Flats 

0.5  deg 

0.5  deg 

0.5  deg 

Table  15.  AH-64D  Large  Moves  Ad  Hoc  Comparison.  This  Table  shows  the  largest 
difference  between  large  moves  ad  hoc  and  PC-GBS  adjustment  sets. _ 


AH-64D 

Weight 

Pitch  Link 

Tab  8-10 

Tab  6-10 

Tab  4-10 

Largest  Individual  Difference 

16.369 

0.076 

0.259 

0.106 

0.063 

Largest  RMS  Difference 

9.786 

0.041 

0.151 

0.055 

0.026 

Basic  Adjustment  Unit 

52  grams 

0.25  Flats 

0.5  deg 

0.5  deg 

0.5  deg 

These  tables  show  that  for  the  UH-60  and  the  AH-64D,  the  respective  RMS  values  of  the 
large  moves  coefficients  produced  adjustment  sets  that  were  virtually  identical  to  those  of 
the  PC-GBS.  For  the  AH-64A,  a  comparison  of  Tables  14  and  7  reveal  that  the  RMS  of 
the  small  moves  ad  hoc  coefficients  actually  provides  a  better  match  to  adjustments  of  the 
PC-GBS. 
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Figures  16,  17,  and  18  illustrate  the  magnitude  difference  between  adjustments  as 
calculated  by  the  PC-GBS  and  the  single  RMS  matrix  of  the  appropriate  ad  hoc 
coefficients  for  the  three  helicopters. 


Weight  Pitch  Link  Trim  Tab 


Flight  Number 

Figure  16.  UFI-60  Bar  Chart  Comparison  of  RMS  of  Large 
chart  shows  the  magnitude  difference  between  the  PC-GBS 
on  the  RMS  of  large  move  ad  hoc  coefficients.  The  red  bar 
adjustment  unit  magnitude. 
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Figure  17.  AFI-64A  Bar  Chart  Comparison  of  RMS  of  Small  Moves  to  PC-GBS.  This 
chart  shows  the  magnitude  difference  between  the  PC-GBS  adjustments  and  those  based 
on  the  RMS  of  small  move  ad  hoc  coefficients.  The  red  bar  indicates  the  basic 
adjustment  unit  magnitude. 


Figure  18.  AFI-64D  Bar  Chart  Comparison  of  RMS  of  Large  Moves  to  PC-GBS.  This 
chart  shows  the  magnitude  difference  between  the  PC-GBS  adjustments  and  those  based 
on  the  RMS  of  large  move  ad  hoc  coefficients.  The  red  bar  indicates  the  basic 
adjustment  unit  magnitude. 
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Figures  16-18  show  that  a  single  set  of  linear  coefficients  will  almost  always 
provide  reduced  adjustment  sets  that  are  identical  to  those  calculated  by  the  AEN.  Table 
16  contains  the  standard  deviation  of  the  difference  between  PC-GBS  adjustments  and 
those  produced  by  the  RMS  of  the  appropriate  ad  hoc  coefficients. 


Table  16.  Standard  Deviation  Between  Ad  Hoc  RMS  and  PC-GBS  Adjustments.  This 
table  shows  standard  deviation  of  the  difference  in  adjustments  as  determined  by  the  PC- 
GBS  and  the  RMS  of  the  indicated  ad  hoc  coefficients. 


Weight 

Pitch  Link 

Trim  Tab 

UH-60  Large  Moves  RMS 

0.187  oz 

0.1  Notch 

0.11  Mil 

Weight 

Pitch  Link 

Tab  8-10 

Tab  6-10 

Tab  4-10 

AH-64A  Small  Moves  RMS 

7.86  g 

0.10  Flats 

0.12  deg 

0.10  deg 

0.08  deg 

AH-64D  Large  Moves  RMS 

2.63  g 

0.01  Flats 

0.04  deg 

0.02  deg 

0.01  deg 

These  standard  deviation  values  are  all  less  than  the  size  of  the  basic  adjustment  unit. 
This  means  that  a  vast  majority  of  the  reduced  adjustments,  as  calculated  by  this  RMS 
method,  will  produce  identical  detailed  adjustments  to  those  of  the  PC-GBS. 

4.3.4  Ad  Hoc  Sensitivity  Coefficients  as  Determined  By  PC-GBS. 

Tables  17,  18,  and  19  contain  the  RMS  values  of  the  ad  hoc  sensitivity 
coefficients  for  the  UH-60,  AH-64A,  and  AH-64D  as  determined  in  this  study.  The 
coefficients  for  the  UH-60  and  the  AH-64D  were  detennined  by  taking  the  RMS  of  the 
large  move  ad  hoc  sensitivity  coefficients.  The  coefficients  for  the  AH-64A  where 
determined  by  taking  the  RMS  of  the  small  moves  ad  hoc  sensitivity  coefficients.  These 
coefficients  may  be  encoded  into  a  linear  algorithm  such  as  AVA  in  order  to  produce 
reduced  adjustment  sets  that  are  virtually  identical  to  those  determined  by  the  PC-GBS. 
Of  course  some  form  of  post  processing  will  be  required  in  order  to  convert  the  reduced 
adjustment  sets  to  detailed  adjustment  sets  with  the  same  performance  capabilities  as 
those  of  the  PC-GBS. 
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Table  17.  UH-60  Ad  Hoc  Sensitivity  Coefficients.  These  coefficients  were  detennined 


by  taking  the  RMS  of  the  large  moves  ad 


hoc  coefficients. 


Weight 

Pitch  Link 

Tab 

UH-60 

Mag 

Phase 

Mag 

Phase 

Mag 

Phase 

FPGIOO(A-B) 

0.012437 

335.1 

0.061253 

119.9 

0.014479 

146.7 

Hover(A-B) 

0.003777 

258.1 

0.062257 

140.9 

0.008360 

128.8 

80Kts(A-B) 

0.003672 

261.0 

0.027984 

114.8 

0.005148 

99.7 

mm 

0.000018 

309.2 

0.043495 

187.9 

0.016227 

197.2 

8fet»iagfei:W 

0.003595 

250.2 

0.042692 

127.8 

0.007671 

116.6 

0.000057 

273.8 

0.068672 

176.9 

0.028640 

191.4 

145Kts(A-B) 

0.003700 

249.2 

0.046724 

132.0 

0.008472 

139.5 

145Kts(A+B) 

0.000041 

311.3 

0.081355 

182.9 

0.040795 

190.3 

0.000006 

15.9 

0.013415 

278.9 

0.004451 

243.3 

Table  18.  AH-64A  Ad  Hoc  Sensitivity  Coefficients.  These  coefficients  were 


determined  by  taking  the  RMS  o 


the  small  moves  ad  hoc  coefficients. 


AH-64A 

Pitch  Link 

Tab  8-10 

Tab  6-10 

Tab  4-10 

Magnitude 

i  ps/gram 

ips/flat 

ips/deg 

ips/deg 

ips/deg 

FPGIOO(Lat) 

0.000593 

0.041779 

0.000165 

0.000139 

0.000135 

Hover(Lat) 

0.000515 

0.185371 

0.028842 

0.064399 

0.104124 

0.000358 

0.030565 

0.148199 

0.250506 

0.295098 

HSB 

0.000012 

0.079075 

0.334511 

0.501796 

0.751610 

lOOKt(Vert) 

0.000009 

0.129714 

0.357461 

0.546680 

0.797509 

120Kt(Vert) 

0.000011 

0.159954 

0.440764 

0.660609 

0.969146 

140Kt(Vert) 

0.000018 

0.235513 

0.562146 

0.824348 

1.323498 

Phase 

deg 

deg 

deg 

deg 

deg 

FPGIOO(Lat) 

169.0 

20.7 

338.3 

344.3 

0.8 

Hover(Lat) 

166.4 

54.7 

71.1 

61.6 

51.3 

224.9 

220.0 

247.9 

246.0 

259.4 

44.0 

268.5 

251.5 

264.2 

262.3 

msmbesi 

101.2 

267.4 

261.6 

256.1 

260.3 

120Kt(Vert) 

42.5 

248.6 

252.5 

256.0 

260.2 

140Kt(Vert) 

22.9 

246.1 

243.6 

249.8 

238.2 

4-15 


Table  19.  AH-64D  Ad  Hoc  Sensitivity  Coefficients.  These  coefficients  were 


determined  by  taking  the  RMS  o 


the  large  moves  ad  hoc  coefficients. 


AH-64D 

Weight 

Pitch  Link 

Tab  8-10 

Tab  6-10 

Tab  4-10 

Magnitude 

ips/gram 

ips/flat 

ips/deg 

ips/deg 

ips/deg 

FPGIOI(Lat) 

0.000491 

0.044240 

0.000000 

0.000000 

0.000000 

Hover(Lat) 

0.000456 

0.155225 

0.000000 

0.000000 

0.000000 

60Kt(Vert) 

0.000468 

0.038996 

0.159772 

0.336765 

0.647302 

0.000452 

0.063262 

0.170617 

0.287236 

0.735197 

mmssm 

0.000480 

0.113166 

0.187719 

0.314955 

0.655976 

120Kt(Vert) 

0.000451 

0.180973 

0.214903 

0.372574 

0.695178 

140Kt(Vert) 

0.000449 

0.242108 

0.305237 

0.444347 

0.903616 

Phase 

deg 

deg 

deg 

deg 

deg 

FPGIOI(Lat) 

163.0 

15.5 

165.9 

165.9 

165.9 

Hover(Lat) 

171.1 

57.6 

238.3 

238.9 

238.9 

60Kt(Vert) 

212.0 

286.3 

263.6 

270.5 

256.0 

204.7 

273.2 

261.6 

261.3 

259.2 

216.0 

262.1 

258.7 

268.3 

256.7 

219.7 

256.1 

255.8 

258.7 

255.0 

235.7 

247.3 

250.3 

260.4 

250.6 

4.5  Chapter  Summary 

The  Adjustment  Evaluation  Networks  and  the  Vibration  Prediction  Network  of 
the  PC-GBS  have  been  isolated  and  studied  in  detail.  Graphical  and  statistical  analyses 
have  shown  that  in  all  vibration  categories,  the  AEN  calculates  adjustment  sets  using  a 
linear  mapping.  The  prediction  capabilities  of  the  VPN  have  also  been  shown  to  behave 
in  a  purely  linear  fashion  with  respect  to  adjustment  magnitude.  This  overall  lack  of  non¬ 
linear  behavior  indicates  that  the  neural  networks  of  the  PC-GBS  have  learned  that  a 
linear  relationship  exists  between  adjustments  and  vibration  response,  thus  validating  the 
linear  assumption. 
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V.  Conclusions  and  Recommendations 


5.1  Chapter  Overview 

Due  to  the  difficulty  in  making  MRS  adjustments  on  many  helicopters,  it  was 
strongly  believed  that  the  linear  assumption  was  flawed.  For  this  reason,  neural  networks 
were  used  to  determine  a  non-linear  mapping  of  the  solution  space  of  main  rotor 
vibrations  in  order  to  produce  more  accurate  adjustments.  This  study  has  shown  that  the 
mapping  of  these  networks  is  essentially  linear.  Therefore,  the  performance  gains  of  the 
PC-GBS  are  not  due  to  unique,  non-linear  mappings,  but  rather  to  accurate  linear 
mappings  and  improved  post-processing  of  the  reduced  adjustment  sets. 

5.2  Conclusions  of  Research 

The  goal  of  this  study  was  to  characterize  the  vibration/adjustment  mapping  as  it 
is  known  to  a  trained  neural  network.  This  goal  was  achieved  through  analysis  of  the 
AEN  and  VPN  in  the  PC-GBS.  By  studying  multiple  flights  over  a  broad  range  of 
vibration  and  adjustment  magnitudes,  it  was  determined  that  the  PC-GBS  mappings  can 
almost  always  be  described  accurately  with  linear  sensitivity  coefficients. 

In  a  study  conducted  by  Wroblewski  et  al  [1],  the  PC-GBS  algorithm 
outperformed  the  AVA  algorithm  by  consistently  producing  adjustment  sets  for  the  AH- 
64  with  fewer  moves  and  lower  predicted  vibrations.  It  is  now  apparent  that  the  success 
that  the  PC-GBS  has  enjoyed  is  not  due  to  the  inclusion  of  higher-order  interactions  in 
the  vibration  response,  but  rather  to  improved  accuracy  in  a  traditional  linear  mapping, 
the  effects  of  the  SOE,  and  other  improvements. 
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The  SOE  has  made  large  improvements  in  decreasing  the  chances  of  human  error  in  main 
rotor  smoothing  iterations  by  selecting  adjustment  sets  with  minimal  numbers  of  moves. 
While  there  is  a  slight  tradeoff  in  predicted  vibration  levels,  this  is  more  than  justified 
when  these  few  adjustment  moves  are  applied  correctly  the  first  time.  Future  rotor 
smoothing  algorithms  must  also  account  for  the  human  element  in  order  to  surpass  the 
performance  of  the  PC-GBS.  Until  the  day  that  adjustment  mechanisms  are  engineered 
to  be  error  proof,  mistakes  will  continue  to  reduce  main  rotor  smoothing  performance. 

5.3  Significance  of  Research 

Helicopter  main  rotor  smoothing  is  an  extremely  expensive  and  time  consuming 
task  that  must  be  periodically  performed  on  all  helicopters.  Small  improvements  to  the 
process  have  the  potential  to  save  millions  of  dollars  annually.  In  recent  years,  efforts 
have  been  made  to  design  non-linear  algorithms  for  the  task  of  main  rotor  smoothing. 

This  research  has  shown  that  a  set  of  non-linear  neural  networks  have  essentially  learned 
that  the  vibration  response  is  linear.  Assuming  that  the  neural  networks  in  the  PC-GBS 
are  fully  trained,  the  significance  of  this  research  is  that  the  linear  assumption  is 
completely  adequate  for  calculating  adjustments.  Software  designers  can  safely  continue 
to  use  the  linear  assumption  and  seek  improvements  to  MRS  perfonnance  elsewhere. 
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This  research  has  also  shown  that  a  fair  amount  of  disparity  exists  between  the 
adjustments  offered  by  the  AVA  algorithm  and  those  of  the  PC-GBS,  which  is  now 
known  to  be  linear.  Due  to  the  improved  performance  that  has  been  noted  by 
Wroblewski  et  al,  it  is  believed  that  the  linear  mapping  of  the  PC-GBS  is  superior  to  that 
of  AVA  for  the  AH-64A  and  the  UH-60.  The  mapping  for  the  AH-64D  is,  of  course, 
identical. 

The  US  Army  currently  performs  MRS  procedures  on  these  three  aircraft  with  the 
Aviation  Vibration  Analyzer  and,  to  a  lesser  extent,  the  PC-GBS.  Until  such  time  as  the 
PC-GBS  is  used  on  all  Anny  helicopters,  the  sensitivity  coefficients  that  have  been 
determined  in  this  research  are  offered  as  a  free  upgrade  to  the  AVA  systems.  An 
improvement  in  perfonnance  should  be  expected  from  incorporating  these  updated 
coefficients. 

5,4  Recommendations  for  Action 

One  other  area  where  the  PC-GBS  has  shown  great  potential  in  improving  MRS 
performance  is  the  incorporation  of  the  Solution  Optimization  Expert.  The  SOE 
improves  the  MRS  process  by  selecting  adjustment  sets  with  minimal  numbers  of  moves. 
This  not  only  shortens  maintenance  time,  but  also  minimizes  the  chances  for  error. 

It  is  well  known  that  human  error  is  a  common  occurrence  during  adjustment 
application.  For  this  reason,  new  adjustment  mechanisms  must  be  devised  that  reduce  or 
eliminate  the  chances  for  mistakes.  Something  as  simple  as  an  engraved  numbering 
scheme  on  the  pitch  links,  much  like  on  a  micrometer,  could  virtually  eliminate  the 
chances  of  applying  pitch  adjustments  in  the  wrong  direction  or  on  the  wrong  blade. 
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This  research  has  shown  that  the  linear  assumption  is  still  valid;  however,  linear 
algorithms  must  still  be  replaced.  By  their  very  nature,  linear  algorithms  require  that 
specific  flight  profiles  be  flown  while  vibration  and  track  data  are  recorded.  Pilot  error 
can  lead  to  poor  data  acquisition  and  therefore,  poor  adjustment  sets.  A  rotor  smoothing 
process  that  utilizes  continuous  vibration  measurements  could  potentially  alleviate  the 
requirement  of  strict  adherence  to  flight  regimes.  This  type  of  system  has  been  examined 
by  Branhof  et  al  [9]  and  has  been  detennined  to  provide  adjustments  that  improve 
vibration  levels.  Further  efforts  should  be  made  to  develop  algorithms  that  do  not  require 
that  specific  flight  profiles  be  flown  while  recording  vibration  data. 

5.5  Summary 

There  are  many  improvements  that  still  need  to  be  made  to  the  main  rotor 
smoothing  process.  The  ultimate  goal  is  to  one  day  have  a  system  that  can  produce 
adjustment  sets  from  a  single  set  of  flight  data  that  will  reduce  any  vibrations  to  within 
acceptable  levels.  This  thesis  has  reported  on  one  of  the  latest  efforts  in  MRS  algorithm 
modernization.  While  the  use  of  neural  networks  in  MRS  applications  has  proven  to  be 
beneficial,  the  underlying  method  that  the  PC-GBS  neural  networks  uses  for  mapping  the 
solution  space  is  not  unique. 


5-4 


Appendix  A:  AVA  Sensitivity  Coefficients 


UH60  Coefficients 

adjustment  =  Hub  Weight,  Pitch  Link,  Trim  Tabs 
unit  =  Oz,  Notches,  Mils 

Hover  A  -  B 

Coeff  =  {0.0  IPS/oz,  0  deg},  (0.046  IPS/notch,  147.3  deg},  {0.0162  IPS/mil,  146.7  deg} 
80Kts  A  +  B 

Coeff  =  {0.0  IPS/oz,  0  deg},  {0.0410  IPS/notch,  196.2  deg},  {0.0296  IPS/mil,  196.2  deg} 
80Kts  A  -  B 

Coeff  =  {0.0  IPS/oz,  0  deg},  {0.0289  IPS/notch,  126.7  deg},  {0.0105  IPS/mil,  122.4  deg} 
1 20Kts  A  +  B 

Coeff  =  {0.0  IPS/oz,  0  deg},  {0.0516  IPS/notch,  189.6  deg},  {0.0413  IPS/mil,  191.7  deg} 
1 20Kts  A  -  B 

Coeff  =  {0.0  IPS/oz,  0  deg},  {0.0369  IPS/notch,  138.0  deg},  {0.0113  IPS/mil,  136.7  deg} 
145Kts  A  +  B 

Coeff  =  {0.0  IPS/oz,  0  deg},  {0.066  IPS/notch,  192.7  deg},  {0.053  IPS/mil,  192.0  deg} 
145Kts  A  -  B 

Coeff  =  {0.0  IPS/oz,  0  deg},  {0.046  IPS/notch,  147.0  deg},  {0.018  IPS/mil,  139.0  deg} 


AH64A  Coefficients 

adjustment  =  Hub  Weight,  Pitch  Link,  Tab  8-10,  Tab  6-10,  Tab  4-10 
units  =  Grams,  Flats,  Degrees,  Degrees,  Degrees 

FPG100  LAT 

Coeff  =  {.00033  IPS/gram,  165°},  {.044  I  PS/flat, 21°},  {0  IPS/0,  0°},  {0,0},  {0,0} 

Hover  LAT 

Coeff  =  {.00037  IPS/gram,  169°},  {.143  IPS/flat,  51°},  {0  IPS  1°,  0°},  {0,0},  {0,0} 

60KVERT 

Coeff  =  {.000331  IPS/gram,  231°},  {.054  IPS/flat,  274°},  {.141  IPS/0,  255°},  {.255  IPS/0,  246°},  {.286  IPS/0,  273°} 
80KVERT 

Coeff  =  {.000269  IPS/gram,  234°},  {.062  IPS/flat,  283°},  {.227  IPS/0,  267°},  {.329  IPS/0,  260°},  {.363  IPS/0,  270°} 
100KVERT 

Coeff  =  {.000390  IPS/gram,  232°},  {.106  IPS/flat,  265°},  {.264  IPS/0,  269°},  {.434  IPS/0,  260°},  {.485  IPS/0,  273°} 
120K  VERT 

Coeff  =  {.000369  IPS/gram,  242°},  {.156  IPS/flat,  240°},  {.405  IPS/0,  262°},  {.613  IPS/0,  253°},  {.630  IPS/0,  266°} 
140KVERT 

Coeff  =  {.000287  IPS/gram,  250°},  {.224  IPS/flat,  239°},  {.436  IPS/0,  249°},  {.664  IPS/0,  245°},  {.689  IPS/0,  267°} 
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AH64D  Coefficients 

adjustment  =  Hub  Weight,  Pitch  Link,  Tab  8-10,  Tab  6-10,  Tab  4-10 
units  =  Grams,  Flats,  Degrees,  Degrees,  Degrees 

FPG100  LAT 

Coeff  =  {.0004936  IPS/gram,  163.0°},  {.04448  flat/0,  15.5°},  {0  IPS/0,  0°},  {0,  0°},  {0,  0°} 

Hover  LAT 

Coeff  =  {.0004580  IPS/gram,  171.2°},  {.1560  flat/0,  57.5°},  {0  IPS/0,  0°},  {0,  0},  {0,  0°} 

60KTA  VERT 

Coeff  =  {.0004696  IPS/gram,  21 1 .9  °},  {.03924  flat/0,  286.6°},  {.1605  IPS  1°,  263.5°},  {.3385  IPS/0,  270.6°},  {.6507  IPS/0, 
256.1°} 

80KTA  VERT 

Coeff  =  {.0004546  IPS/gram,  204.7°},  {.06359  flat/0,  273.2°},  {.1714  IPS/0,  261.6°},  {.2886  IPS/0,  261.3°},  {.7391  IPS/0, 
259.2°} 

100KTA  VERT 

Coeff  =  {.0004825  IPS/gram,  215.9°},  {.1137  flat/0,  262.0°},  {.1888  IPS/0,  258.8°},  {.3168  IPS/0,  268.4°},  {.6596  IPS/0, 
256.7°} 

120KTA  VERT 

Coeff  =  {.0004547  IPS/gram,  219.6°},  {.1819  flat/0,  256.1°},  {.2161  IPS/0,  255.8°}  {.3747  IPS/0,  258.6°},  {.6992  IPS/0, 
255.0°} 

140KTA  VERT 

Coeff  =  {.0004519  IPS/gram,  235.6°},  {.2434  flat/0,  247.2°},  {.3071  IPS/0,  250.3°},  {.4467  IPS/0,  260.5°},  {.9085  IPS/0, 
250.5°} 
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Appendix  B:  UH-60  Database 


Note:  All  data  in  this  appendix  was  manually  typed  into  Microsoft  Excel  spreadsheets. 
The  data  was  originally  presented  with  a  graphic  user  interface  in  the  PC-GBS  program. 


UH-60  Flight  Log 

Flight  # 

Type 

BUNO 

Date 

Time 

Vib  Category 

1 

UH-60 

79-23280 

5/14/2005 

103937 

good 

2 

UH-60 

81-23550 

11/27/2003 

104807 

good 

3 

UH-60 

80-23470 

2/15/2005 

161139 

good 

4 

UH-60 

84-23981 

8/20/2004 

115128 

good 

5 

UH-60 

95-26658 

2/28/2004 

100113 

good 

6 

UH-60 

80-23470 

10/26/2003 

112939 

above 

7 

UH-60 

80-23470 

2/13/2005 

130352 

above 

8 

UH-60 

79-23280 

4/2/2003 

181046 

above 

9 

UH-60 

79-23280 

8/22/2003 

102945 

above 

10 

UH-60 

79-23280 

10/9/2003 

100510 

above 

11 

UH-60 

79-23280 

10/21/2003 

103923 

above 

12 

UH-60 

84-23981 

2/28/2003 

181616 

above 

13 

UH-60 

84-23981 

1/5/2004 

44045 

above 

14 

UH-60 

83-23921 

12/2/2005 

162722 

caution 

15 

UH-60 

80-23470 

2/13/2005 

145852 

caution 

16 

UH-60 

80-23470 

2/13/2005 

121244 

caution 

17 

UH-60 

95-26659 

3/17/2004 

102442 

caution 

18 

UH-60 

81-23550 

8/4/2003 

114453 

caution 

19 

UH-60 

84-23981 

11/23/2003 

81228 

exceed 

20 

UH-60 

91-26330 

9/16/2003 

132022 

exceed 

VMEP  Solution 

Flight  # 

Yellow 

Blue 

Red 

Black 

Wt 

P/L 

Tab 

Wt 

P/L 

Tab 

Wt 

P/L 

Tab 

Wt 

P/L 

Tab 

1 

7 

0 

0 

-6 

-2 

0 

0 

1 

-6 

0 

0 

-7 

2 

0 

0 

2 

0 

0 

-7 

0 

1 

0 

0 

-2 

0 

3 

-7 

-1 

0 

-5 

0 

0 

0 

0 

0 

0 

0 

0 

4 

0 

-1 

3 

0 

0 

4 

0 

0 

0 

0 

1 

0 

5 

0 

0 

0 

10 

0 

0 

5 

-2 

0 

0 

0 

5 

6 

16 

0 

0 

-11 

-3 

0 

0 

0 

0 

0 

0 

0 

7 

-11 

0 

0 

-16 

0 

-4 

0 

-2 

0 

0 

-3 

0 

8 

20 

0 

6 

16 

0 

3 

0 

-2 

0 

0 

0 

0 

9 

11 

0 

-2 

0 

2 

0 

0 

-1 

0 

-13 

0 

10 

10 

0 

0 

0 

0 

0 

-3 

0 

-2 

2 

0 

-3 

0 

11 

0 

0 

0 

0 

0 

-4 

6 

-2 

9 

-6 

-3 

0 

12 

-14 

0 

-8 

9 

2 

3 

0 

-1 

0 

0 

0 

0 

13 

0 

4 

0 

-6 

0 

0 

-44 

0 

6 

0 

0 

0 

14 

9 

0 

0 

7 

0 

0 

0 

0 

-16 

0 

0 

-2 

15 

-12 

0 

0 

-14 

0 

0 

0 

1 

0 

0 

-2 

-11 

16 

-16 

0 

0 

-26 

0 

0 

0 

-4 

9 

0 

-5 

20 

17 

-11 

-2 

7 

-15 

0 

0 

0 

0 

0 

0 

-4 

20 

18 

-8 

0 

0 

-9 

0 

0 

0 

0 

10 

0 

-4 

19 

19 

-28 

0 

0 

-16 

0 

-16 

0 

0 

20 

0 

0 

0 

20 

0 

0 

0 

-8 

0 

15 

0 

0 

20 

0 

0 

0 

B-l 


Flight# 


1 


FpgIOO  A-B 


Flover  A-B 


80Kt  A-B 


120Kt  A-B 


145Kt  A-B 


Flover  A+B 


80Kt  A+B 


120Kt  A+B 


Before 

Mag 


0.032707 


0.083903 


0.045817 


0.121328 


0.064676 


0.164067 


0.117083 


0.149976 


Full  Soln 


Pred 

Mag 


0.0197 


0.01359 


0.047175 


0.072691 


0.02042 


0.156546 


0.074428 


0.008687 


Pred  Yellow 
Phase  Wt 


145Kt  A+B  0.181986 


After 

Mag 


0.069268 


0.065411 


0.048443 


0.103722 


0.047172 


0.164286 


0.116287 


0.147891 


After  Yellow 
Phase  P/L 


After 

Mag 


0.051184 


0.1288 


0.067723 


0.149222 


0.100951 


0.177728 


0.091119 


0.080251 


After  Yellow 
Phase  Tab 


After 

Mag 


0.041347 


0.096912 


0.051565 


0.132964 


0.075388 


0.172191 


0.100755 


0.096581 


Flight# 


2 


FpgIOO  A-B 


Flover  A-B 


80Kt  A-B 


120Kt  A-B 


145Kt  A-B 


Flover  A+B 


80Kt  A+B 


120Kt  A+B 


145Kt  A+B 


Vib 

Mag 


0.123726 


0.027149 


0.004418 


0.029976 


0.04869 


0.163819 


0.125195 


0.183337 


0.176624 


Pred  Yellow 
Phase  Wt 


0.047083 


0.052165 


0.02832 


0.02966 


0.067846 


0.155879 


0.132697 


0.125869 


0.089133 


After 

Mag 


0.081083 


0.041329 


0.013692 


0.032046 


0.066435 


0.163762 


0.124163 


0.183109 


0.178737 


After  Yellow 
Phase  P/L 


After 

Mag 


0.18681 


0.081229 


0.031764 


0.073063 


0.042052 


0.162791 


0.09655 


0.197626 


0.243284 


After  Yellow 
Phase  Tab 


After 

Mag 


0.152567 


0.037668 


0.014868 


0.045286 


0.042265 


0.157802 


0.09819 


0.179883 


0.238748 


Flight# 


3 


FpgIOO  A-B 


Flover  A-B 


80Kt  A-B 


120Kt  A-B 


145Kt  A-B 


Flover  A+B 


80Kt  A+B 


120Kt  A+B 


145Kt  A+B 


Vib 

Mag 


0.03638 


0.084397 


0.141022 


0.188868 


0.104637 


0.081093 


0.078236 


0.067702 


0.139082 


Pred  Yellow 
Phase  Wt 


0.023388 


0.127113 


0.121409 


0.171848 


0.093494 


0.087003 


0.044569 


0.078086 


0.117409 


After 

Mag 


0.096599 


0.066325 


0.156057 


0.205598 


0.090338 


0.08104 


0.079445 


0.068724 


0.140772 


After  Yellow 
Phase  P/L 


After 

Mag 


0.036389 


0.112495 


0.134744 


0.188196 


0.147993 


0.076929 


0.120591 


0.114155 


0.198099 


After  Yellow 
Phase  Tab 


After 

Mag 


0.005196 


0.091682 


0.135442 


0.184501 


0.119648 


0.083628 


0.111787 


0.113317 


0.192785 


FpgIOO  A-B 


Flover  A-B 


80Kt  A-B 


120Kt  A-B 


145Kt  A-B 


Flover  A+B 


80Kt  A+B 


120Kt  A+B 


145Kt  A+B 


Vib 

Mag 


0.065383 


0.030869 


0.042057 


0.088015 


0.024951 


0.098902 


0.046546 


0.152539 


0.087242 


Pred  Yellow 
Phase  Wt 


0.035052 


0.081646 


0.025012 


0.069855 


0.028618 


0.103693 


0.027918 


0.101525 


0.014692 


After 

Mag 


0.117679 


0.047527 


0.024201 


0.087869 


0.020008 


0.098717 


0.046047 


0.15071 


0.085445 


After  Yellow 
Phase  P/L 


After 

Mag 


0.078001 


0.041575 


0.068372 


0.129982 


0.04203 


0.088973 


0.046151 


0.093166 


0.078107 


After  Yellow 
Phase  Tab 


After 

Mag 


0.050358 


0.015436 


0.052784 


0.102513 


0.020133 


0.089629 


0.047083 


0.112091 


0.087801 


Flight# 


5 


FpgIOO  A-B 


Flover  A-B 


80Kt  A-B 


120Kt  A-B 


145Kt  A-B 


Flover  A+B 


80Kt  A+B 


120Kt  A+B 


145Kt  A+B 


Vib 

Mag 


0.164886 


0.110756 


0.138877 


0.165934 


0.15523 


0.048921 


0.195477 


0.073415 


0.249527 


Yellow  After  After  Yellow  After  After  Yellow  After  After 


0.052168 


0.0843 


0.016932 


0.126838 


IHMEQEai 


0.062149 


Wt  |  Mag 


0.217559 


0.120421 


0.145008 


0.175768 


0.158113 


0.048873 


0.195211 


0.071677 


0.247478 


Mag  Phase  Tab  Mag  I  Phase 


0.151814 


0.047315 


0.117165 


0.122132 


0.110025 


0.045825 


0.195851 


0.053381 


0.193362 


0.142606 


0.09254 


0.132097 


0.149506 


0.136924 


0.051937 


0.189549 


0.036529 


0.1844 


B-2 


Flight# 


6 


FpgIOO  A-B 


Hover  A-B 


80Kt  A-B 


120Kt  A-B 


145Kt  A-B 


Hover  A+B 


80Kt  A+B 


120Kt  A+B 


145Kt  A+B 


Vib 

Mag 


0.071631 


0.235183 


0.121482 


0.163799 


0.283022 


0.125287 


0.40637 


0.27158 


0.085521 


Pred  Yellow 
Phase  Wt 


0.039252 


0.073489 


0.027723 


0.013991 


0.116078 


0.088157 


0.284325 


0.062613 


0.167117 


After 

Mag 


0.020056 


0.218362 


0.105876 


0.146516 


0.266013 


0.12524 


0.406641 


0.271377 


0.084935 


After  Yellow 
Phase  P/L 


After 

Mag 


0.123017 


0.240397 


0.136762 


0.181446 


0.288889 


0.121108 


0.422615 


0.282599 


0.102317 


After  Yellow 
Phase  Tab 


After 

Mag 


0.101456 


0.238165 


0.128917 


0.171929 


0.28187 


0.127871 


0.413406 


0.265621 


0.09478 


FpgIOO  A-B 


Hover  A-B 


80Kt  A-B 


120Kt  A-B 


145Kt  A-B 


Hover  A+B 


80Kt  A+B 


120Kt  A+B 


Vib 

Mag 


0.343202 


0.091605 


0.224173 


0.39421 


0.110259 


0.091364 


0.080271 


0.265887 


Pred  Yellow 
Phase  Wt 


145Kt  A+B  0.139524 


0.046506 


0.161325 


0.078613 


0.219135 


0.110074 


0.092323 


0.084629 


0.15114 


0.088816 


After 

Mag 


0.367909 


0.108966 


0.238623 


0.411574 


0.124712 


0.091223 


0.080538 


0.265406 


0.136674 


After  Yellow 
Phase  P/L 


After 

Mag 


0.394559 


0.143072 


0.233878 


0.409215 


0.106653 


0.101648 


0.08742 


0.296279 


0.218119 


After  Yellow 
Phase  Tab 


After 

Mag 


0.32991 


0.076931 


0.218869 


0.38694 


0.119414 


0.089653 


0.094859 


0.265814 


0.05463 


Flight# 


8 


FpgIOO  A-B 


Hover  A-B 


80Kt  A-B 


120Kt  A-B 


145Kt  A-B 


Hover  A+B 


80Kt  A+B 


120Kt  A+B 


Vib 

Mag 


0.090777 


0.110572 


0.212216 


0.32375 


0.018327 


0.111057 


0.361304 


0.461297 


Pred  Yellow 
Phase  Wt 


145Kt  A+B  0.333519 


0.023738 


0.134221 


0.13177 


0.176085 


0.16691 


0.138663 


0.189272 


0.166561 


0.161466 


After 

Mag 


0.05125 


0.101542 


0.219209 


0.328891 


0.024277 


0.111222 


0.360281 


0.45943 


0.330687 


After  Yellow 
Phase  P/L 


After 

Mag 


0.153956 


0.188618 


0.280973 


0.0336 


0.122806 


0.323352 


0.396855 


0.249109 


After  Yellow 
Phase  Tab 


After 

Mag 


0.06588 


0.116635 


0.21861 


0.336991 


0.033726 


0.108312 


0.386371 


0.507014 


0.413915 


Flight# 


9 


FpgIOO  A-B 


Hover  A-B 


80Kt  A-B 


120Kt  A-B 


145Kt  A-B 


Hover  A+B 


80Kt  A+B 


120Kt  A+B 


145Kt  A+B 


Vib 

Mag 


0.247293 


0.029988 


0.154695 


0.075135 


0.073574 


0.172685 


0.118899 


0.078203 


0.311925 


Pred  Yellow 
Phase  Wt 


0.053574 


0.114688 


0.10163 


0.035242 


0.017912 


0.209098 


0.180323 


0.072629 


0.084439 


After 

Mag 


0.301827 


0.048492 


0.149204 


0.06992 


0.075849 


0.172932 


0.118026 


0.079619 


0.31296 


After  Yellow 
Phase  P/L 


After 

Mag 


0.310507 


0.087807 


0.1328 


0.033244 


0.040392 


0.185022 


0.092378 


0.137788 


0.349783 


After  Yellow 
Phase  Tab 


After 

Mag 


0.275717 


0.046248 


0.147883 


0.06002 


0.056884 


0.177226 


0.102162 


0.120858 


0.34161 


Flight# 


10 


FpgIOO  A-B 


Hover  A-B 


80Kt  A-B 


120Kt  A-B 


145Kt  A-B 


Hover  A+B 


80Kt  A+B 


120Kt  A+B 


145Kt  A+B 


Vib 

Mag 


0.141422 


0.20833 


0.17198 


0.284234 


0.119453 


0.155949 


0.372482 


0.167834 


0.070084 


Vib 

Phase 

Pred  Pred 

Mag  Phase 

Yellow 

Wt 

After 

Mag 

After 

Phase 

Yellow  After 

P/L  Mag 

After 

Phase 

Yellow 

Tab 

After 

Mag 

After 

Phase 

0.146776 


0.174379 


0.222382 


0.189913 


0.299123 


0.137225 


0.156142 


0.3721 


0.16723 


0.072039 


0.094671 


0.238706 


0.150418 


0.246274 


0.122686 


0.168902 


0.356858 


0.153804 


0.143597 


0.128042 


0.212487 


0.161925 


0.268491 


0.120867 


0.161464 


0.365229 


0.168196 


0.14812 


B-3 


Flight# 


11 


FpgIOO  A-B 


Hover  A-B 


80Kt  A-B 


120Kt  A-B 


145Kt  A-B 


Hover  A+B 


80Kt  A+B 


120Kt  A+B 


145Kt  A+B 


Vib 

Mag 


0.094397 


0.096156 


0.189777 


0.15312 


0.144485 


0.223468 


0.276624 


0.199556 


0.30721 


Pred  Yellow 
Phase  Wt 


0.028446 


0.124352 


0.139536 


0.093446 


0.091786 


0.228018 


0.230493 


0.035045 


0.156248 


After 

Mag 


0.061013 


0.112297 


0.205902 


0.167124 


0.158739 


0.223659 


0.276358 


0.200924 


0.310082 


After  Yellow 
Phase  P/L 


After 

Mag 


0.128226 


0.126183 


0.161227 


0.11329 


0.166455 


0.236355 


0.266003 


0.244749 


0.39244 


After  Yellow 
Phase  Tab 


After 

Mag 


0.118311 


0.098419 


0.178636 


0.13684 


0.153684 


0.228745 


0.273132 


0.246428 


0.393314 


Flight# 


12 


FpgIOO  A-B 


Hover  A-B 


80Kt  A-B 


120Kt  A-B 


145Kt  A-B 


Hover  A+B 


80Kt  A+B 


120Kt  A+B 


145Kt  A+B 


Vib 

Mag 


0.088372 


0.138989 


0.074038 


0.104176 


0.189213 


0.023186 


0.246438 


0.278277 


0.322668 


Pred  Yellow 
Phase  Wt 


0.040087 


0.044676 


0.03868 


0.065874 


0.070583 


0.048864 


0.087278 


0.064167 


0.07209 


After 

Mag 


0.142768 


0.156993 


0.087231 


0.121074 


0.207116 


0.022969 


0.247341 


0.279056 


0.324228 


After  Yellow 
Phase  P/L 


After 

Mag 


0.025421 


0.148198 


0.076244 


0.07831 1 


0.166941 


0.009461 


0.27691 


0.303205 


0.378857 


After  Yellow 
Phase  Tab 


After 

Mag 


0.061214 


0.135446 


0.070957 


0.09079 


0.181655 


0.016388 


0.272864 


0.310735 


0.386504 


Flight# 


13 


FpgIOO  A-B 


Hover  A-B 


80Kt  A-B 


120Kt  A-B 


145Kt  A-B 


Hover  A+B 


80Kt  A+B 


120Kt  A+B 


145Kt  A+B 


Vib 

Mag 


0.457242 


0.147058 


0.171925 


0.214342 


0.13403 


0.160661 


0.176284 


0.20973 


0.10178 


Pred  Yellow 
Phase  Wt 


0.028405 


0.053746 


0.071035 


0.103124 


0.118282 


0.131253 


0.097421 


0.131168 


0.079618 


After 

Mag 


0.405699 


0.142737 


0.155019 


0.197622 


0.137117 


0.160441 


0.175006 


0.207565 


0.100729 


After  Yellow 
Phase  P/L 


After 

Mag 


0.489631 


0.105347 


0.189928 


0.225867 


0.08882 


0.147 


0.131059 


0.14083 


0.109462 


After  Yellow 
Phase  Tab 


After 

Mag 


0.482902 


0.135795 


0.18043 


0.220478 


0.115646 


0.152469 


0.142106 


0.149039 


0.118944 


Flight# 


14 


FpgIOO  A-B 


Hover  A-B 


80Kt  A-B 


120Kt  A-B 


145Kt  A-B 


Hover  A+B 


80Kt  A+B 


120Kt  A+B 


145Kt  A+B 


Vib 

Mag 


0.043871 


0.250402 


0.076905 


0.164148 


0.334668 


0.228697 


0.444824 


0.401959 


0.644317 


Pred 

Mag 


0.077443 


0.131275 


0.02997 


0.171274 


0.185128 


0.172602 


0.233458 


0.102478 


0.075764 


Pred  Yellow 
Phase  Wt 


After 

Mag 


0.103587 


0.263887 


0.085938 


0.181595 


0.339695 


0.228478 


0.443627 


0.399825 


0.641461 


After  Yellow 
Phase  P/L 


After 

Mag 


0.029493 


0.192677 


0.053517 


0.156085 


0.287052 


0.215184 


0.405011 


0.33457 


0.560175 


After  Yellow 
Phase  Tab 


After 

Mag 


0.012464 


0.232501 


0.069111 


0.156969 


0.316082 


0.220307 


0.412136 


0.341561 


0.560864 


Flight#  |  Vib  Vib 
Mag  Phase 


FpgIOO  A-B 


Hover  A-B 


80Kt  A-B 


120Kt  A-B 


145Kt  A-B 


Hover  A+B 


80Kt  A+B 


120Kt  A+B 


145Kt  A+B 


0.363445 


0.177773 


0.256493 


0.409914 


0.214618 


0.074173 


0.304634 


0.579386 


0.479439 


mzvawmmsMi 


0.128494  4’ 


EEB  HEBEI  ESB 


Pred  Yellow  After  After  Yellow  After  After  Yellow  After  After 

Phase  Wt  Mag  Phase  P/L  Mag  Phase  Tab  Mag  Phase 


0.18547 


0.152435  1 


0.137422 


0.372994 


0.189922 


0.26965 


0.42562 


0.217979 


0.074143 


0.304827 


0.579373 


0.47938 


0.328144 


0.216621 


0.253309 


0.412123 


0.25717 


0.071581 


0.309864 


0.574319 


0.489705 


0.358174 


0.184672 


0.252327 


0.40729 


0.231814 


0.077653 


0.313194 


0.58879 


0.499899 


B-4 


Flight# 


16 


FpgIOO  A-B 


Hover  A-B 


80Kt  A-B 


120Kt  A-B 


145Kt  A-B 


Hover  A+B 


80Kt  A+B 


120Kt  A+B 


145Kt  A+B 


Vib 

Mag 


0.374782 


0.161809 


0.222461 


0.418754 


0.04873 


0.163275 


0.155919 


0.416578 


0.624346 


Pred  Yellow 
Phase  Wt 


0.090561 


0.158173 


0.045509 


0.24027 


0.212789 


0.086611 


0.113375 


0.260885 


0.279851 


After 

Mag 


0.403091 


0.178094 


0.239399 


0.435896 


0.048966 


0.163106 


0.156862 


0.417116 


0.624957 


After  Yellow 
Phase  P/L 


After 

Mag 


0.3231 


0.116625 


0.207602 


0.390419 


0.027232 


0.151703 


0.19397 


0.44748 


0.642935 


After  Yellow 
Phase  Tab 


After 

Mag 


0.359229 


0.145586 


0.214348 


0.40572 


0.03286 


0.159739 


0.181803 


0.428989 


0.633334 


Flight# 


17 


FpgIOO  A-B 


Hover  A-B 


80Kt  A-B 


120Kt  A-B 


145Kt  A-B 


Hover  A+B 


80Kt  A+B 


120Kt  A+B 


145Kt  A+B 


Vib 

Mag 


0.031618 


0.117707 


0.073185 


0.076669 


0.224682 


0.033688 


0.232903 


0.399486 


0.574942 


Pred  Yellow 
Phase  Wt 


0.050796 


0.079119 


0.041995 


0.059912 


0.117649 


0.069537 


0.096333 


0.084833 


0.089234 


After 

Mag 


0.090624 


0.13068 


0.080352 


0.082855 


0.2376 


0.033768 


0.23267 


0.399681 


0.575259 


After  Yellow 
Phase  P/L 


After 

Mag 


0.033769 


0.157628 


0.08699 


0.112975 


0.2482 


0.042207 


0.233676 


0.420277 


0.585438 


After  Yellow 
Phase  Tab 


After 

Mag 


0.005643 


0.124324 


0.075076 


0.087716 


0.234957 


0.034119 


0.227422 


0.402528 


0.575537 


Flight# 


18 


FpgIOO  A-B 


Hover  A-B 


80Kt  A-B 


120Kt  A-B 


145Kt  A-B 


Hover  A+B 


80Kt  A+B 


120Kt  A+B 


145Kt  A+B 


Vib 

Mag 


0.10891 


0.097986 


0.058938 


0.065231 


0.119797 


0.115527 


0.144347 


0.450975 


0.709512 


Pred  Yellow 
Phase  Wt 


0.050044 


0.056474 


0.003668 


0.023473 


0.060929 


0.144143 


0.099493 


0.091338 


0.080946 


After 

Mag 


0.049047 


0.107456 


0.066904 


0.075608 


0.117981 


0.115383 


0.145216 


0.452588 


0.711889 


After  Yellow 
Phase  P/L 


After 

Mag 


0.166536 


0.148636 


0.072991 


0.096004 


0.167161 


0.108682 


0.180683 


0.511043 


0.778737 


After  Yellow 
Phase  Tab 


After 

Mag 


0.140336 


0.108782 


0.060728 


0.073046 


0.138326 


0.106987 


0.168551 


0.493166 


0.773768 


Flight# 


19 


FpgIOO  A-B 


Hover  A-B 


80Kt  A-B 


120Kt  A-B 


145Kt  A-B 


Hover  A+B 


80Kt  A+B 


120Kt  A+B 


145Kt  A+B 


Vib 

Mag 


0.065843 


0.198795 


0.095744 


0.226021 


0.321385 


0.164591 


0.588989 


0.947801 


1 .537488 


Pred  Yellow 
Phase  Wt 


0.148752 


0.175332 


0.077975 


0.123039 


0.255459 


0.056255 


0.199789 


0.314875 


0.547287 


After 

Mag 


0.111037 


0.180366 


0.07789 


0.209508 


0.303864 


0.164717 


0.590158 


0.949888 


1 .540229 


After  Yellow 
Phase  P/L 


After 

Mag 


0.025188 


0.220988 


0.119454 


0.260607 


0.353312 


0.171195 


0.632033 


1.01837 


1.617413 


After  Yellow 
Phase  Tab 


After 

Mag 


0.050843 


0.206389 


0.105878 


0.240341 


0.331272 


0.172808 


0.619995 


1 .004095 


1.614839 


Flight  #  |  Vib  Vib 

Mag  Phase 


Pred  Yellow  After  After  Yellow  After  After  Yellow  After  After 

Phase  Wt  Mag  Phase  P/L  Mag  Phase  Tab  Mag  Phase 


FpgIOO  A-B 


Hover  A-B 


80Kt  A-B 


120Kt  A-B 


145Kt  A-B 


Hover  A+B 


80Kt  A+B 


120Kt  A+B 


145Kt  A+B 


0.442209 


0.174997 


0.138393 


0.135376 


0.341969 


0.049651 


0.459304 


0.678619 


1.151947 


0.133645 


0.154564 


0.103606  1 


0.036912  4 


0.405373 


0.18026 


0.122097 


0.125465 


0.34101 


0.04972 


0.46043 


0.680173 


1.15381 


0.454472 


0.230026 


0.166499 


0.178895 


0.388051 


0.053708 


0.49719 


0.725206 


1.212255 


0.460267 


0.188225 


0.149506 


0.151953 


0.36026 


0.056565 


0.490632 


0.727371 


1.219836 


B-5 


Before  Before 


Pred  Yellow 


After  Yellow 


After  Yellow 


FpgIOO  A-B 


Hover  A-B 


80Kt  A-B 


120Kt  A-B 


145Kt  A-B 


Hover  A+B 


80  Kt  A+B 


120Kt  A+B 


145Kt  A+B 


Mag  Phase  Mag  Phase  Wt  Mag  Phase  P/L  Mag  Phase  Tab  Mag  |  Phase 


0.032707 


0.083903 


0.045817 


0.121328 


0.064676 


0.164067 


0.117083 


0.149976 


0.181986 


0.0197 


0.01359 


0.047175 


0.072691 


0.02042 


0.156546 


0.074428 


0.008687 


0.087178 


0.993553 


0.217932 


0.293496 


0.165645 


0.231657 


0.164172 


0.118343 


0.151708 


0.183913 


1.813004 


1.909229 


0.867441 


1.345842 


1 .437952 


0.565151 


1.222679 


1.904632 


2.252873 


0.285169 


0.232944 


0.127949 


0.252069 


0.208562 


0.244803 


0.272916 


0.434758 


0.634939 


Flight#  |  Vib  Vib 
Mag  Phase 


Pred  Yellow  After  After  Yellow  After  After  Yellow  After  After 

Phase  Wt  Mag  Phase  P/L  Mag  Phase  Tab  Mag  Phase 


FpgIOO  A-B 


Hover  A-B 


80Kt  A-B 


120Kt  A-B 


145Kt  A-B 


Hover  A+B 


80  Kt  A+B 


120Kt  A+B 


145Kt  A+B 


0.123726 


0.027149 


0.004418 


0.029976 


0.04869 


0.163819 


0.125195 


0.183337 


0.176624 


0.047083 


0.052165 


0.02832 


0.02966 


0.067846 


0.155879 


0.132697 


0.125869 


0.089133 


0.894825 


0.317759 


0.289636 


0.282459 


0.343268 


0.164222 


0.124965 


0.181511 


0.173296 


1.960873 


1.881627 


0.842267 


1.309011 


1.369656 


0.416319 


1.215092 


2.071132 


2.564042 


0.405038 


0.17746 


0.106583 


0.179457 


0.149647 


0.123108 


0.230021 


0.558862 


0.931555 


Flight#  |  Vib  Vib 
Mag  Phase 


Pred  Yellow  After  After  Yellow  After  After  Yellow  After  After 

Phase  Wt  Mag  Phase  P/L  Mag  Phase  Tab  Mag  Phase 


FpgIOO  A-B 


Hover  A-B 


80Kt  A-B 


120Kt  A-B 


145Kt  A-B 


Hover  A+B 


80  Kt  A+B 


120Kt  A+B 


145Kt  A+B 


0.03638 


0.084397 


0.141022 


0.188868 


0.104637 


0.081093 


0.078236 


0.067702 


0.139082 


0.023388 


0.127113 


0.121409 


0.171848 


0.093494 


0.087003 


0.044569 


0.078086 


0.117409 


1.03471 


0.225102 


0.416202 


0.467965 


0.215097 


0.080646 


0.07812 


0.068932 


0.135666 


1.813185 


1.888799 


0.807597 


1.272416 


1 .497649 


0.380723 


1.380062 


2.091228 


2.527629 


0.254791 


0.212508 


0.119095 


0.199781 


0.260065 


0.133662 


0.4014 


0.614206 


0.886971 


Flight#  |  Vib  Vib 
Mag  Phase 


Pred  Yellow  After  After  Yellow  After  After  Yellow  After  After 

Phase  Wt  Mag  Phase  P/L  Mag  Phase  Tab  Mag  Phase 


FpgIOO  A-B 


Hover  A-B 


80Kt  A-B 


120Kt  A-B 


145Kt  A-B 


Hover  A+B 


80  Kt  A+B 


120Kt  A+B 


145Kt  A+B 


0.065383 


0.030869 


0.042057 


0.088015 


0.024951 


0.098902 


0.046546 


0.152539 


0.087242 


0.035052 


0.081646 


0.025012 


0.069855 


0.028618 


0.103693 


0.027918 


0.101525 


0.014692 


1 .047449 


0.330151 


0.252749 


0.287559 


0.283846 


0.099073 


0.047845 


0.154667 


0.09065 


1.826668 


1.849551 


0.878326 


1.366173 


1.394808 


0.335922 


1.285933 


1.928613 


2.399551 


0.252719 


0.139282 


0.144195 


0.231705 


0.156444 


0.013274 


0.312348 


0.46573 


0.780793 


Flight#  |  Vib  Vib 
Mag  Phase 


Pred  Yellow  After  After  Yellow  After  After  Yellow  After  After 

Phase  Wt  Mag  Phase  P/L  Mag  Phase  Tab  Mag  Phase 


FpgIOO  A-B 


Hover  A-B 


80Kt  A-B 


120Kt  A-B 


145Kt  A-B 


Hover  A+B 


80  Kt  A+B 


120Kt  A+B 


145Kt  A+B 


0.164886 


0.110756 


0.138877 


0.165934 


0.15523 


0.048921 


0.195477 


0.073415 


0.249527 


0.052168 


0.0843 


0.040756 


0.057836 


0.016932 


0.035917 


0.126838 


0.104117 


0.062149 


1.135672 


0.371266 


0.362199 


0.405185 


0.356016 


0.048475 


0.194455 


0.07329 


0.249061 


1.781108 


1.755722 


0.732769 


1.112223 


1.25242 


0.38724 


1.296957 


2.002688 


2.254299 


0.198906 


0.059665 


0.112171 


0.02774 


0.031173 


0.11165 


0.331686 


0.507608 


0.621972 


B-6 


Pred  Yellow 


After  Yellow 


After  Yellow 


FpgIOO  A-B 


Hover  A-B 


80Kt  A-B 


120Kt  A-B 


145Kt  A-B 


Hover  A+B 


80  Kt  A+B 


120Kt  A+B 


145Kt  A+B 


Mag  Phase  Mag  Phase  Wt  Mag  Phase  P/L  Mag  Phase  Tab  Mag  |  Phase 


0.071631 


0.235183 


0.121482 


0.163799 


0.283022 


0.125287 


0.40637 


0.27158 


0.085521 


0.039252 


0.073489 


0.027723 


0.013991 


0.116078 


0.088157 


0.284325 


0.062613 


0.167117 


0.923969 


0.150736 


0.195814 


0.13535 


0.122263 


0.124841 


0.405112 


0.269737 


0.083311 


1.882477 


1.869642 


0.899455 


1.336653 


1.440219 


0.376096 


1.481608 


2.08241 


2.412573 


0.356794 


0.310739 


0.205243 


0.272268 


0.317302 


0.171045 


0.558549 


0.5765 


0.784034 


Pred  Yellow 


After  Yellow 


After  Yellow 


FpgIOO  A-B 


Hover  A-B 


80Kt  A-B 


120Kt  A-B 


145Kt  A-B 


Hover  A+B 


80  Kt  A+B 


120Kt  A+B 


145Kt  A+B 


Mag  |  Phase 


0.343202 


0.091605 


0.224173 


0.39421 


0.110259 


0.091364 


0.080271 


0.265887 


0.139524 


Phase  |  Wt  |  Mag  Phase  P/L  Mag  Phase  Tab  Mag  |  Phase 


0.046506 


0.161325 


0.078613 


0.219135 


0.110074 


0.092323 


0.084629 


0.15114 


0.088816 


1.154576 


0.388697 


0.487856 


0.676483 


0.382061 


0.090976 


0.081306 


0.268098 


0.141326 


2.129105 


1.934472 


0.928616 


1 .449609 


1.372509 


0.473521 


1.298521 


2.168018 


2.577311 


0.32885 


0.099537 


0.190742 


0.347157 


0.237481 


0.11232 


0.352912 


0.603242 


0.67606 


Flight#  |  Vib  Vib 
Mag  Phase 


Pred  Yellow  After  After  Yellow  After  After  Yellow  After  After 

Phase  Wt  Mag  Phase  P/L  Mag  Phase  Tab  Mag  Phase 


FpgIOO  A-B 


Hover  A-B 


80Kt  A-B 


120Kt  A-B 


145Kt  A-B 


Hover  A+B 


80  Kt  A+B 


120Kt  A+B 


145Kt  A+B 


0.090777 


0.110572 


0.212216 


0.32375 


0.018327 


0.111057 


0.361304 


0.461297 


0.333519 


0.023738 


0.134221 


0.13177 


0.176085 


0.16691 


0.138663 


0.189272 


0.166561 


0.161466 


0.92103 


0.26753 


0.423857 


0.497562 


0.295825 


0.111411 


0.362454 


0.463386 


0.335067 


1.930995 


1.812448 


0.67382 


0.970806 


1.387674 


0.499837 


1.011703 


1.649025 


2.112659 


0.375149 


0.162829 


0.163615 


0.195454 


0.152936 


0.165589 


0.236642 


0.349527 


0.483904 


Flight#  |  Vib  Vib 
Mag  Phase 


Pred  Yellow  After  After  Yellow  After  After  Yellow  After  After 

Phase  Wt  Mag  Phase  P/L  Mag  Phase  Tab  Mag  Phase 


FpgIOO  A-B 


Hover  A-B 


80Kt  A-B 


120Kt  A-B 


145Kt  A-B 


Hover  A+B 


80  Kt  A+B 


120Kt  A+B 


145Kt  A+B 


0.247293 


0.029988 


0.154695 


0.075135 


0.073574 


0.172685 


0.118899 


0.078203 


0.311925 


0.053574 


0.114688 


0.10163 


0.035242 


0.017912 


0.209098 


0.180323 


0.072629 


0.084439 


0.806763 


0.274036 


0.37334 


0.321469 


0.300273 


0.173007 


0.120161 


0.075867 


0.3086 


2.088473 


1.890617 


0.725277 


1 .208787 


1.342468 


0.561767 


1 .22529 


2.117212 


2.568264 


0.523627 


0.193434 


0.123841 


0.089399 


0.106645 


0.227578 


0.272465 


0.615087 


0.935508 


Flight#  |  Vib  Vib 
Mag  Phase 


Pred  Yellow  After  After  Yellow  After  After  Yellow  After  After 

Phase  Wt  Mag  Phase  P/L  Mag  Phase  Tab  Mag  Phase 


FpgIOO  A-B 


Hover  A-B 


80Kt  A-B 


120Kt  A-B 


145Kt  A-B 


Hover  A+B 


80  Kt  A+B 


120Kt  A+B 


145Kt  A+B 


0.141422 


0.20833 


0.17198 


0.284234 


0.119453 


0.155949 


0.372482 


0.167834 


0.070084 


0.092054 


0.191414 


0.092927 


0.17768 


0.058056 


0.176863 


0.272648 


0.01736 


0.146776 


1 .060332 


0.468413 


0.466076 


0.552079 


0.411222 


0.156229 


0.373768 


0.170091 


0.070879 


1.722304 


1 .94552 


0.712456 


1.037024 


1.391404 


0.550626 


1.20532 


1.997786 


2.495251 


0.247485 


0.284686 


0.086278 


0.143805 


0.206292 


0.218743 


0.428946 


0.568251 


0.8742 


B-7 


Pred  Yellow 


After  Yellow 


After  Yellow 


FpgIOO  A-B 


Hover  A-B 


80Kt  A-B 


120Kt  A-B 


145Kt  A-B 


Hover  A+B 


80  Kt  A+B 


120Kt  A+B 


145Kt  A+B 


Mag  Phase  Mag  Phase  Wt  Mag  Phase  P/L  Mag  Phase  Tab  Mag  |  Phase 


0.094397 


0.096156 


0.189777 


0.15312 


0.144485 


0.223468 


0.276624 


0.199556 


0.30721 


0.028446 


0.124352 


0.139536 


0.093446 


0.091786 


0.228018 


0.230493 


0.035045 


0.156248 


0.928254 


0.383688 


0.473525 


0.42181 


0.411954 


0.223759 


0.277892 


0.200341 


0.305386 


1.874974 


1.897601 


0.656232 


1.144496 


1.458929 


0.617465 


1.251084 


2.178004 


2.758963 


0.36207 


0.192219 


0.087989 


0.026426 


0.266749 


0.282001 


0.392955 


0.729814 


1.125045 


Pred  Yellow 


After  Yellow 


After  Yellow 


FpgIOO  A-B 


Hover  A-B 


80Kt  A-B 


120Kt  A-B 


145Kt  A-B 


Hover  A+B 


80  Kt  A+B 


120Kt  A+B 


145Kt  A+B 


Mag  |  Phase 


0.088372 


0.138989 


0.074038 


0.104176 


0.189213 


0.023186 


0.246438 


0.278277 


0.322668 


Phase  |  Wt  |  Mag  Phase  P/L  Mag  Phase  Tab  Mag  |  Phase 


0.040087 


0.044676 


0.03868 


0.065874 


0.070583 


0.048864 


0.087278 


0.064167 


0.07209 


1.069339 


0.43392 


0.346404 


0.388874 


0.483466 


0.02305 


0.246874 


0.279688 


0.323924 


1.747431 


1.857802 


0.831048 


1.203984 


1.299999 


0.378737 


1 .467644 


2.137508 


2.639592 


0.212009 


0.182418 


0.101003 


0.08389 


0.187353 


0.071079 


0.53305 


0.742155 


1 .062665 


Flight#  |  Vib  Vib 
Mag  Phase 


Pred  Yellow  After  After  Yellow  After  After  Yellow  After  After 

Phase  Wt  Mag  Phase  P/L  Mag  Phase  Tab  Mag  Phase 


FpgIOO  A-B 


Hover  A-B 


80Kt  A-B 


120Kt  A-B 


145Kt  A-B 


Hover  A+B 


80  Kt  A+B 


120Kt  A+B 


145Kt  A+B 


0.457242 


0.147058 


0.171925 


0.214342 


0.13403 


0.160661 


0.176284 


0.20973 


0.10178 


0.028405 


0.053746 


0.071035 


0.103124 


0.118282 


0.131253 


0.097421 


0.131168 


0.079618 


0.626516 


0.306239 


0.14673 


0.120104 


0.345596 


0.160571 


0.17695 


0.210795 


0.105101 


2.075305 


1.762414 


0.952014 


1.33749 


1 .278389 


0.244898 


1.131844 


1.859153 


2.416867 


0.69725 


0.135084 


0.255815 


0.300237 


0.037634 


0.094956 


0.137755 


0.335527 


0.765277 


Flight#  |  Vib  Vib 
Mag  Phase 


Pred  Yellow  After  After  Yellow  After  After  Yellow  After  After 

Phase  Wt  Mag  Phase  P/L  Mag  Phase  Tab  Mag  Phase 


FpgIOO  A-B 


Hover  A-B 


80Kt  A-B 


120Kt  A-B 


145Kt  A-B 


Hover  A+B 


80  Kt  A+B 


120Kt  A+B 


145Kt  A+B 


0.043871 


0.250402 


0.076905 


0.164148 


0.334668 


0.228697 


0.444824 


0.401959 


0.644317 


0.077443 


0.131275 


0.02997 


0.171274 


0.185128 


0.172602 


0.233458 


0.102478 


0.075764 


1.03779 


0.516871 


0.335205 


0.448978 


0.511539 


0.228629 


0.444938 


0.402841 


0.646459 


1.797667 


1.636209 


0.772848 


1 .236283 


1.070727 


0.181253 


0.928988 


1.679193 


1.799874 


0.246057 


0.089559 


0.072465 


0.158519 


0.165439 


0.153421 


0.166951 


0.171379 


0.244248 


Flight#  |  Vib  Vib 
Mag  Phase 


Pred  Yellow  After  After  Yellow  After  After  Yellow  After  After 

Phase  Wt  Mag  Phase  P/L  Mag  Phase  Tab  Mag  Phase 


FpgIOO  A-B 


Hover  A-B 


80Kt  A-B 


120Kt  A-B 


145Kt  A-B 


Hover  A+B 


80  Kt  A+B 


120Kt  A+B 


145Kt  A+B 


0.363445 


0.177773 


0.256493 


0.409914 


0.214618 


0.074173 


0.304634 


0.579386 


0.479439 


0.044457 


0.128494 


0.109086 


0.195748 


0.18547 


0.152435 


0.037687 


0.144156 


0.137422 


1 .083352 


0.427051 


0.506661 


0.672653 


0.378148 


0.073728 


0.305708 


0.581379 


0.482247 


1.630054 


1.956539 


0.831888 


1.338673 


1.586923 


0.387663 


1.349217 


2.060193 


2.50061 


0.411868 


0.276902 


0.23179 


0.406464 


0.375719 


0.133231 


0.487299 


0.857859 


1.010004 


B-8 


Pred  Yellow 


After  Yellow 


After  Yellow 


FpgIOO  A-B 


Hover  A-B 


80Kt  A-B 


120Kt  A-B 


145Kt  A-B 


Hover  A+B 


80  Kt  A+B 


120Kt  A+B 


145Kt  A+B 


Mag  Phase  Mag  Phase  Wt  Mag  Phase  P/L  Mag  Phase  Tab  Mag  |  Phase 


0.374782 


0.161809 


0.222461 


0.418754 


0.04873 


0.163275 


0.155919 


0.416578 


0.624346 


0.090561 


0.158173 


0.045509 


0.24027 


0.212789 


0.086611 


0.113375 


0.260885 


0.279851 


1.195644 


0.451571 


0.506889 


0.702885 


0.294016 


0.162933 


0.154722 


0.414355 


0.621235 


1.526207 


1.735992 


0.732074 


1.037558 


1.358156 


0.275524 


1.429665 


2.243364 


2.606872 


0.328622 


0.067048 


0.157888 


0.307763 


0.126949 


0.151998 


0.447117 


0.748348 


1.039969 


Pred  Yellow 


After  Yellow 


After  Yellow 


FpgIOO  A-B 


Hover  A-B 


80Kt  A-B 


120Kt  A-B 


145Kt  A-B 


Hover  A+B 


80  Kt  A+B 


120Kt  A+B 


145Kt  A+B 


Mag  |  Phase 


0.031618 


0.117707 


0.073185 


0.076669 


0.224682 


0.033688 


0.232903 


0.399486 


0.574942 


Phase  |  Wt  |  Mag  Phase  P/L  Mag  Phase  Tab  Mag  |  Phase 


0.050796 


0.079119 


0.041995 


0.059912 


0.117649 


0.069537 


0.096333 


0.084833 


0.089234 


1 .024037 


0.383469 


0.320983 


0.311485 


0.47158 


0.034126 


0.231857 


0.397408 


0.571986 


1.805875 


1.921829 


0.864868 


1.335737 


1.503437 


0.419734 


1.307397 


2.178215 


2.53441 


0.258641 


0.227202 


0.129634 


0.206848 


0.343935 


0.091715 


0.352028 


0.685962 


0.961197 


Flight#  |  Vib  Vib 
Mag  Phase 


Pred  Yellow  After  After  Yellow  After  After  Yellow  After  After 

Phase  Wt  Mag  Phase  P/L  Mag  Phase  Tab  Mag  Phase 


FpgIOO  A-B 


Hover  A-B 


80Kt  A-B 


120Kt  A-B 


145Kt  A-B 


Hover  A+B 


80  Kt  A+B 


120Kt  A+B 


145Kt  A+B 


0.10891 


0.097986 


0.058938 


0.065231 


0.119797 


0.115527 


0.144347 


0.450975 


0.709512 


0.050044 


0.056474 


0.003668 


0.023473 


0.060929 


0.144143 


0.099493 


0.091338 


0.080946 


0.887043 


0.349833 


0.315864 


0.322076 


0.290945 


0.115804 


0.14311 


0.448725 


0.706326 


1.933096 


1.932445 


0.857654 


1.314334 


1.517209 


0.352808 


1.415628 


2.442037 


3.034649 


0.398909 


0.231435 


0.120082 


0.185505 


0.28878 


0.046329 


0.43109 


0.931503 


1.411913 


Flight#  |  Vib  Vib 
Mag  Phase 


Pred  Yellow  After  After  Yellow  After  After  Yellow  After  After 

Phase  Wt  Mag  Phase  P/L  Mag  Phase  Tab  Mag  Phase 


FpgIOO  A-B 


Hover  A-B 


80Kt  A-B 


120Kt  A-B 


145Kt  A-B 


Hover  A+B 


80  Kt  A+B 


120Kt  A+B 


145Kt  A+B 


0.065843 


0.198795 


0.095744 


0.226021 


0.321385 


0.164591 


IIifjfcMiTikll 


0.947801 


1 .537488 


0.148752 


0.175332 


0.077975 


0.123039 


0.255459 


0.056255 


0.199789 


0.314875 


0.547287 


1.031178 


0.119829 


0.19798 


0.097675 


0.068969 


0.164279 


0.58801 


0.946059 


1.534801 


1.773744 


1.917394 


0.914788 


1.459156 


1.621236 


0.497889 


1.872481 


3.000183 


3.918124 


0.25029 


0.308531 


0.194438 


0.366807 


0.434272 


0.24661 


0.890396 


1.491208 


2.292662 


Flight#  |  Vib  Vib 
Mag  Phase 


Pred  Yellow  After  After  Yellow  After  After  Yellow  After  After 

Phase  Wt  Mag  Phase  P/L  Mag  Phase  Tab  Mag  Phase 


FpgIOO  A-B 


Hover  A-B 


80Kt  A-B 


120Kt  A-B 


145Kt  A-B 


Hover  A+B 


80  Kt  A+B 


120Kt  A+B 


145Kt  A+B 


0.442209 


0.174997 


0.138393 


0.135376 


0.341969 


0.049651 


0.459304 


0.678619 


1.151947 


0.064965 


0.077556 


0.10312 


0.133645 


0.154564 


0.115483 


0.103606 


0.036912 


0.134108 


0.777295 


0.372328 


0.171421 


0.224545 


0.417944 


0.04926 


0.459367 


0.679114 


1.152754 


1.934348 


2.010072 


0.973705 


1.413227 


1.727227 


0.413426 


1.70118 


2.533648 


3.341251 


0.646767 


0.312077 


0.24117 


0.287996 


0.509111 


0.128947 


0.765271 


1.183963 


1 .858797 


B-9 


Appendix  C:  AH-64A  Database 


Note:  All  data  in  this  appendix  was  manually  typed  into  Microsoft  Excel  spreadsheets. 
The  data  was  originally  presented  with  a  graphic  user  interface  in  the  PC-GBS  program. 


Flight  Log 

Flight  # 

Type 

BUNO 

Date 

Time 

Vib  Category 

1 

AH-64A 

86-08964 

6/5/2003 

153548 

good 

2 

AH-64A 

86-08962 

11/14/2003 

132815 

good 

3 

AH-64A 

86-08964 

6/5/2003 

145701 

good 

4 

AH-64A 

86-08963 

8/11/2003 

94214 

good 

5 

AH-64A 

91-00113 

1/28/2005 

35232 
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Appendix  D:  AH-64D  Database 


Note:  All  data  in  this  appendix  was  manually  typed  into  Microsoft  Excel  spreadsheets. 
The  data  was  originally  presented  with  a  graphic  user  interface  in  the  PC-GBS  program. 


Flight  Log 


Flight  # 

Type 

BUNO 

Date 

Time 

Vib  Category 

1 

AH-64D 

00-05184 

6/8/2005 

25649 

good 

2 

AH-64D 

01-05276 

9/29/2005 

50641 

good 

3 

AH-64D 

01-05283 

10/14/2005 

81022 

good 

4 

AH-64D 

96-05024 

9/15/2005 

44633 

good 

5 

AH-64D 

99-05134 

10/18/2005 

95536 

good 

6 

AH-64D 

96-05024 

9/15/2005 

40530 

above 

7 

AH-64D 

99-05134 

10/18/2005 

40724 

above 

8 

AH-64D 

00-05184 

8/20/2005 

102558 

above 

9 

AH-64D 

01-05283 

10/13/2005 

101630 

above 

10 

AH-64D 

96-05024 

9/7/2005 

72518 

above 

11 

AH-64D 

00-05184 

8/20/2005 

105108 

caution 

12 

AH-64D 

01-05283 

10/13/2005 

95431 

caution 

13 

AH-64D 

01-05276 

9/29/2005 

43748 

caution 

14 

AH-64D 

01-05277 

10/18/2005 

55723 

caution 

15 

AH-64D 

01-05277 

11/22/2005 

25833 

caution 

16 

AH-64D 

01-05283 

10/13/2005 

92020 

exceed 

17 

AH-64D 

96-05134 

6/9/2005 

23649 

exceed 

18 

AH-64D 

01-05276 

9/29/2005 

21844 

exceed 

19 

AH-64D 

01-05277 

10/18/2005 

31631 

exceed 

20 

AH-64D 

01-05277 

10/18/2005 

43503 

exceed 
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CO 

CD 

O 

CM 

CD 

CD 

CO 

CD 

h- 
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CD 
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CO 
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CD 
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LO 
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LO 
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CO 
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0.937375 

CD 

h- 

h- 

CM 

CD 

CO 

CD 

| 

0.94531 

#1  Blade 
Wt 

Pred 

Phase 

'=t 

CD 

LO 

3- 

co 

CO 

CM 

110 

CD 

CD 

LO 

Pred 

Mag 

0.105433 
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LO 

CD 
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00 
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CD 
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in 
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I 
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CO 

CO 
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n 

o 

d 

CD 

CD 

O 

CD 

O 
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in 

CD 

CD 

CD 

p 
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After 

Phase 

CD 

in 

00 

CO 

CD 
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In 
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CD 

in 
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in 
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CD 

in 

CN 

After 
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CO 

CO 

s 

o 

d 

CD 

CD 

o 

CD 

O 

d 

2.307009 

h- 

CD 
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#1  Blade 
Tab  6 

in 
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Phase 

CD 

in 

co 

in 
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CN] 

o 

in 

CNJ 

CD 

CN] 

CD 

Nj- 

CN 

After 

Mag 

CO 

CO 

Nj- 

n 

o 

d 

CD 

CD 

O 

Nj- 

CD 

O 

d 

r-- 

CD 

s 

1 .693232 

1 .730929 

o 

CD 

CO 

CN 

CD 

CD 

p 

CN 

#1  Blade 

Tab  8 

in 
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Phase 

CO 

CNJ 

CD 

in 

o 

CD 

CNJ 

CNJ 

in 

CNJ 

in 

CNJ 

in 

CN] 
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CN 

After 

Mag 

0.444271 

in 

CD 

CD 

cq 

1 .068535 

1 .554797 

2.136197 

CO 

r-- 

co 

in 

p 

CNj 

3.799389 

0 

TD 
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5 

CNJ 
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Phase 

1 

I 

After 

Mag 

0.604825 

r- 

CNJ 

CD 

CD 

cq 

d 

r-- 

CD 

r-- 

in 

1.31363 

1.272377 

1.256579 

1.348435 

#1  Blade 
Wt 

r- 

o 

Pred 

Phase 

trgt 
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CD 

CNJ 

o 

Is- 

1781- 

CD 

in 

Pred 
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0.105433 

0.064069 

0.151942 

0.130403 

0.027601 

CD 

CD 

O 

CD 

CD 

O 

d 

0.170182 
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Phase 

CD 

in 

CO 

CO 

^r 

CNJ 

CO 
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CN] 

CD 

CO 

CN 

CN 

II 
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d 

0.064069 
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d 
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CD 

d 
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d 
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CD 
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CD 

d 
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Appendix  E:  UH-60  Matlab  Code 


Note:  The  Matlab  codes  contained  in  this  appendix  were  each  used  to  perform  one  or 
more  specific  functions  in  the  course  of  this  research.  There  are  several  instances  where 
parts  of  the  code  have  been  commented  out  with  the  %  symbol.  In  order  to  reproduce  all 
of  the  analysis  of  this  thesis,  some  lines  of  code  may  need  to  be  un-commented  and  other 
lines  commented  out. 


%  Capt  Nathan  A  Miller 

%  This  program  creates  MAT  data  files  from  imported  Excel  database 
%  Ensure  that  the  two  databases  are  imported  into  matlab  prior  to 
running 

function  Excel  to  Mat 
global  data 
global  data2 
clc 


n  =  length (data2  (:, 1 ))  ; 
for  j  =1 : n 

%  Read  Vibration  Magnitude  and  Phase  as  well  as  prediction  magnitude 
%  and  phase  from  Excel  File  and  assign  to  data  matrix 


UH60VIBES ( j , 1) 

=  data (2+ ( j -1 ) 

*13,2)  ; 

UH60VIBES ( j , 2) 

=  data (2+ ( j -1 ) 

*13,3)  ; 

UH60VIBES ( j , 3) 

=  data (2+ ( j -1 ) 

*13,4)  ; 

UH60VIBES ( j , 4) 

=  data (2+ ( j -1 ) 

*13,5)  ; 

UH60VIBES ( j , 5) 

=  data (3+ ( j -1 ) 

*13,2)  ; 

UH60VIBES ( j , 6) 

=  data (3+ ( j -1 ) 

*13,3)  ; 

UH60VIBES ( j , 7) 

=  data (3+ ( j -1 ) 

*13,4)  ; 

UH60VIBES ( j , 8) 

=  data (3+ ( j -1 ) 

*13,5)  ; 

UH60VIBES ( j , 9) 

=  data (4+ ( j -1 ) 

*13,2)  ; 

UH60VIBES ( j , 10) 

=  data (4+ ( j -1 

)  *13,3) ; 

UH60VIBES ( j , 11) 

=  data (4+ ( j -1 

)  *13,4) ; 

UH60VIBES ( j , 12) 

=  data (4+ ( j -1 

)  *13,5) ; 

UH60VIBES ( j , 13) 

=  data (8+ ( j -1 

) *13,2)  ; 

UH60VIBES ( j , 14) 

=  data (8+ ( j -1 

) *13,3)  ; 

UH60VIBES ( j , 15) 

=  data (8+ ( j -1 

)  *13,4) ; 

UH60VIBES ( j , 16) 

=  data (8+ ( j -1 

)  *13,5) ; 

UH60VIBES ( j , 17) 

=  data (5+ ( j -1 

) *13,2)  ; 

UH60VIBES ( j , 18) 

=  data (5+ ( j -1 

) *13,3)  ; 

UH60VIBES ( j , 19) 

=  data (5+ ( j -1 

)  *13,4) ; 

UH60VIBES ( j , 20) 

=  data (5+ ( j -1 

)  *13,5) ; 

UH60VIBES ( j , 21) 

=  data (9+ ( j -1 

) *13,2)  ; 

UH60VIBES ( j , 22) 

=  data (9+ ( j -1 

) *13,3)  ; 

UH60VIBES ( j , 23) 

=  data (9+ ( j -1 

)  *13,4) ; 

UH60VIBES ( j , 24) 

=  data (9+ ( j -1 

)  *13,5) ; 

UH60VIBES ( j , 25) 

=  data (6+ ( j -1 

) *13,2)  ; 

UH60VIBES ( j , 26) 

=  data (6+ ( j -1 

) *13,3)  ; 

UH60VIBES ( j , 27) 

=  data (6+ ( j -1 

)  *13,4) ; 

UH60VIBES ( j , 28) 

=  data (6+ ( j -1 

)  *13,5) ; 
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UH60VIBES ( j , 29)  =  data (10+ (j-1) *13,2)  ; 

UH60VIBES ( j , 30)  =  data (10+ (j-1) *13, 3) ; 

UH60VIBES ( j , 31)  =  data (10+ (j-1) *13, 4)  ; 

UH 60 VIBES ( j , 32 )  =  data (10+ (j-1) *13, 5)  ; 

UH60VIBES ( j , 33)  =  data (7+ (j-1) *13,2)  ; 

UH 60 VIBES ( j , 34 )  =  data  (7+ (j-1) *13, 3) ; 

UH60VIBES ( j , 35)  =  data (7+ ( j -1) *13, 4 )  ; 

UH60VIBES ( j , 36)  =  data (7+ ( j -1) *13, 5)  ; 
end 

save  UH60VIBES .mat  UH60VIBES 

UH60  VMEP  Adjustments  =  zeros  (n, 12); 
for  j  =1 : n 

for  k=l : 12 

UH60_VMEP_Adj ustments ( j , k)  =  data2 ( j , k+1 ) ; 
end 

end 

save  UH60  VMEP  Adjustments .mat  UH60  VMEP  Adjustments 

UH60_Coef f_Data  =  zeros  (n, 57); 
for  j  =1 : n 

UH60_Coef f_Data ( j , 1 ) =data (2  + (j-1)  *13, 6); 

UH60_Coef f_Data ( j , 20 ) =data (2+ ( j -1 ) *  13 , 9) ;  %  Records  manual 

adj  ustment 

UH60_Coef f_Data ( j , 39) =data (2+ ( j -1 )  *  13 , 12 )  ;  %  that  produces  delta 

vibes 

for  k=l : 3 

UH60_Coeff_Data ( j ,  [2  3]  +  (k-l)*19)  =  data ( 2+ ( j -1 ) *  13 ,  [ 7  8]  +  (k-l)*3);  %  FPG100 (A-B) 

UH60_Coeff_Data ( j , [4  5]  +  (k-l)*19)  =  data ( 3+ ( j -1 ) *  13 ,  [ 7  8]  +  (k-l)*3);  %  0  Kts (A-B) 

UH60_Coeff_Data ( j , [6  7]  +  (k-l)*19)  =  data ( 4+ ( j -1 ) *  13 ,  [ 7  8]  +  (k-l)*3);  %  80  Kts (A-B) 

UH60_Coeff_Data ( j , [8  9]+(k-l)*19)  =  data (8+ ( j-1) *13, [7  8]+(k-l)*3);  %  80  Kts (A+B) 

UH60_Coeff_Data ( j , [10  11] + (k-1) *19)  =  data (5+ (j-1) *13, [7  8]+(k-l)*3);  %  120  Kts (A-B) 

UH60_Coeff_Data ( j , [12  13]  +  (k-1) *19)  =  data ( 9+ ( j -1 ) *  13 ,  [ 7  8]  +  (k-l)*3);  %  120  Kts (A+B) 

UH60_Coeff_Data ( j ,  [14  15]  +  (k-1) *19)  =  data ( 6+ ( j -1 ) *  13 ,  [ 7  8]  +  (k-1) *3) ;  %  145  Kts (A-B) 

UH60_Coeff_Data ( j ,  [16  17]  +  (k-1) *19)  =  data ( 10+ ( j -1 ) *  13 ,  [ 7  8]  +  (k-1) *3) ;  %  145  Kts (A+B) 

UH60_Coeff_Data ( j ,  [18  19]  +  (k-1) *19)  =  data ( 7+ ( j -1 ) *  13 ,  [ 7  8 ]  +  (k-1 ) *3 ) ;  %  0  Kts (A+B)  Non 
AVA  quantity 
end 

UH60_Coeff_Data ( j ,  [58  59])  =  data (2+ ( j -1 ) *  13 ,  [ 15  16]); 

UH60_Coef f_Data ( j ,  [ 60  61])  =  data ( 3+ ( j -1 ) *  13 ,  [ 15  16]); 

UH60_Coef f_Data ( j ,  [ 62  63])  =  data ( 4+ ( j -1 ) *  13 ,  [ 15  16]); 

UH60_Coef f_Data ( j ,  [ 64  65])  =  data ( 8+ ( j -1 ) *  13 ,  [ 15  16]); 

UH60_Coef f_Data ( j ,  [ 66  67])  =  data ( 5+ ( j -1 ) *  13 ,  [ 15  16]); 

UH60_Coef f_Data ( j ,  [ 68  69])  =  data ( 9+ ( j -1 ) *  13 ,  [ 15  16]); 

UH60_Coef f_Data ( j ,  [ 70  71])  =  data ( 6+ ( j -1 ) *  13 ,  [ 15  16]); 

UH60_Coeff_Data ( j ,  [72  73])  =  data ( 10+ ( j -1 ) *  13 ,  [ 15  16]); 

UH60_Coef f_Data ( j ,  [ 74  75])  =  data ( 7+ ( j -1 ) *  13 ,  [ 15  16]); 

end 

save  UH60  Coeff  Data. mat  UH60  Coeff  Data 


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
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%  Capt  Nathan  A  Miller 
%  Thesis  work 

%  This  program  loads  vibe  data  and  calculates  an  AVA  adjustment  matrix 

function  Vibs2AVA  ADJ 

load  UH60VIBES .mat 
n  =  length (UH60VIBES (:,  1) )  ; 

for  j  =1 : n 

mag  =  UH60VIBES ( j , [ 1  5  9  13  17  21  25  29])'; 
phase  =  UH60VIBES ( j , [2  6  10  14  18  22  26  30])'; 

vibsCmplx  =  mag.* (cos (  phase*pi/180  )  +  i*sin(  phase*pi/180  )  ); 

pmag  =  UH60VIBES ( j , [ 3  7  11  15  19  23  27  31])'; 
pphase  =  UH60VIBES ( j , [4  8  12  16  20  24  28  32])'; 

pvibsCmplx  =  pmag.* (cos (  pphase*pi/180  )  +  i*sin(  pphase*pi/180  )  ); 

delvib  =  -vibsCmplx  +  pvibsCmplx; 

adj(j)  =  LinearUH60 (  delvib  ); 
end 

for  k=l : n 

magnitude ( k, : )  =  adj(k).mag'; 

Phase (k, :)  =  adj ( k) . phase ' ; 

end 

AVAadj .mag  =  magnitude; 

AVAadj . phase  =  Phase; 
save  AVAadj. mat  AVAadj 


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

%  Capt  Miller 

%  This  function  converts  a  complex  vibset  into  an  adjustment  set  based 
%on  the  AVA  sensitivity  coefficients. 

function  adj  =  LinearUH60 (  vibsCmplx  ) 


%  Matrix  of  coefficient  magnitudes 


mag  = 

[0.00993,  0.0, 

o 

o 

O, 

o 

FPG100 

(A 

- 

B 

0.0,  0.0460, 

0.0162; 

o, 

o 

0 

kts 

(A 

- 

B 

0.0,  0.0289, 

0.0105; 

o, 

o 

80 

kts 

(A 

- 

B 

0.0,  0.0410, 

0.0296; 

o, 

o 

80 

kts 

(A 

+ 

B 

0.0,  0.0369, 

0.0113; 

g, 

o 

120 

kts 

(A 

- 

B 

0.0,  0.0516, 

0.0413; 

o, 

o 

120 

kts 

(A 

+ 

B 

0.0,  0.0460, 

0.0180; 

o, 

o 

145 

kts 

(A 

- 

B 

0.0,  0.0660,  0.0530]; 

%  Matrix  of  coefficient  phases 

%  145 
(deg) 

kts 

(A 

+ 

B 

phase 

=  [338.0,  0.0, 

0.0; 

g, 

o 

FPG100 

(A 

- 

B 

0.0,  147.3, 

146.7; 

g, 

o 

0 

kts 

(A 

- 

B 

0.0,  126.7, 

122.4; 

g, 

o 

80 

kts 

(A 

- 

B 

0.0,  196.2, 

196.2; 

g, 

o 

80 

kts 

(A 

+ 

B 

0.0,  138.0, 

136.7; 

g, 

o 

120 

kts 

(A 

- 

B 
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0.0,  189.6,  191.7;  %  120  kts  (A  +  B) 

0.0,  147.0,  139.0;  %  145  kts  (A  -  B) 

0.0,  192.7,  192.0];  %  145  kts  (A  +  B) 

%  Convert  the  coefficients  to  complex  numbers 
A  =  mag.* (cos (  phase*pi/180  )  +  i*sin(  phase*pi/180  )  ); 

%  Calculate  the  adjustments 
adjCmplx  =  A\vibsCmplx; 

%  Convert  the  adjustments  to  magnitude  and  phase 

adj .mag  =  abs (  adjCmplx  ) ; 

ad j .phase  =  angle (  adjCmplx  ) *180/pi; 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

%  Capt  Nathan  A  Miller 

%  Converts  VMEP  adjustments  into  magnitude  and  phase  of  adjustments. 

function  VMEP  2_MagPhase 

load  UH60  VMEP  Adj ustments . mat 
n  =  length (UH60  VMEP  Adjustments (:, 1) ) ; 

a=UH60  VMEP  Adjustments; 
clear  UH60  VMEP  Adjustments 

for  j  =1 : n 

%  Yellow  Blue  Red  Blck 

adjcmplx  ( j  ,  1)  =  +a (j , 1) -i*a ( j , 4) -a (j , 7) +i*a  (j , 10) ; 
adjcmplx  ( j  ,  2)  =  +a ( j ,2) -i*a (j , 5) -a (j , 8) +i*a  (j , 11) ; 
adjcmplx (j , 3)  =  +a  (j , 3) -i*a  (j , 6) -a  (j , 9) +i*a  (j , 12) ; 

end 

VMEPadj .mag=abs (adjcmplx) ; 

VMEPadj . phase  =  angle (adjcmplx) *180/pi; 

save  VMEPadj .mat  VMEPadj 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

%  Capt  Miller 

%  This  program  converts  the  coefficient  matrix  data  into  actual 
%  coefficient  matrices.  The  program  loads  the  raw  vibe  data  from 
%  UH60Vibes .mat  and  computes  a  delta  vibe  from  the  predicted  vibe 
values 

%  stored  in  UH60  Coeff  Data. mat.  The  coefficient  matrix  is  computed 
%  column  by  column.  Once  computed,  the  coefficient  matrix  operates  on 
a 

%  vibe  set  to  produce  a  linear  adjustment  set. 

function  VMEP^Coef f icient  Matrix 

load  UH60VIBES .mat 
load  UH60  Coeff  Data. mat 

n=length (UH60VIBES ( :  ,  1) ) 
for  j  =1 : n 


%  Weight  Adj 
%  P/L  Adj 
%  Tab  Adj 
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wt  =  UH60_Coef f_Data ( j , 1 ) ; 
pi  =  UH60_Coeff_Data ( j  ,  20)  ; 
tab  =  UH60_Coef f_Data ( j  ,  39)  ; 

for  k=l : 9  %  loops  over  the  8  flight  conditons,  FPG100 (A-B) , 
Hover (A-B) , 

%  80Kts (A-B) (A+B) ,  120Kts (A-B) (A+B) , 145Kts (A-B) (A+B)  and 

finally 

%  the  VMEP  specific  Hover (A+B) condition  as  #9 
vibmag=UH 60 VIBES ( j , 1+4* (k-1) ) ; 
vibphase=UH60VIBES ( j  ,  2  +  4* (k-1)  )  ; 
vibcmplx=vibmag* (cos (  vibphase*pi/180  )  +  ... 

i*sin(  vibphase*pi/180  )  ); 

pvibwtmag=UH60  Coeff  Data (j  ,  2+2* (k-1)  )  ; 
pvibwtphase=UH60_Coef f_Data ( j  ,  3+2* (k-1) )  ; 
pvibwtcmplx=pvibwtmag* (cos (  pvibwtphase*pi/180  )  +  ... 

i*sin(  pvibwtphase*pi/180  )  ); 

pvibplmag=UH60  Coeff  Data (j , 21+2* (k-1) ) ; 
pvibplphase=UH60  Coeff  Data (j , 22+2* (k-1) ) ; 
pvibplcmplx=pvibplmag* (cos (  pvibplphase*pi/180  )  +  ... 

i*sin(  pvibplphase*pi/180  )  ); 

pvibtabmag=UH60  Coeff  Data  (j , 40+2* (k-1) ) ; 
pvibtabphase=UH60  Coeff  Data (j , 41+2* (k-1) )  ; 
pvibtabcmplx=pvibtabmag* (cos (  pvibtabphase*pi/180  )  +  ... 

i*sin(  pvibtabphase*pi/180  )  ); 

dvibwt=pvibwtcmplx-vibcmplx; 
dvibpl=pvibplcmplx-vibcmplx; 
dvibtab=pvibtabcmplx-vibcmplx; 

A ( k, 1 ) =dvibwt/wt ; 

A ( k, 2 ) =dvibpl/pl ; 

A ( k, 3 ) =dvibtab/ tab; 

end 

Coeff { j } =A; 

end 


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 


%  Capt  Miller 

%  This  program  creates  polar  plots  of  the  ad  hoc  coefficient 
%  adjustments  as  well  as  the  VMEP  and  AVA  adjustments.  It  also 
%  calculates  the  maximum  difference  between  the  VMEP  adjustments  and 
%  the  ad  hoc  adjustments. 

clear ; clc; close  all; 

load  Coeff. mat 
load  UH60VIBES .mat 
load  VMEPadj.mat 

%  Runs  the  Least  Squares  Coefficient  function  to  generate  the  LSC 
%  matrix  from  the  Coeff  cell. 

[LSC,  standard] =Least_Squares_Coeff (Coeff )  ; 
r=0 
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for  z=l:5 
r=r+l ; 
vib  num=z; 

b=l  ;%  Sets  beginning  Coefficient  Matrix  to  consider 
e=length (Coef f ) ;  %  Sets  ending  Coefficient  Matrix  to  consider 

%(set  to  length (Coef f)  to  do  all) 

%  Convert  Pre-Adjustment  vibe  data  into  complex  array 
for  j =1 : length (Coef f) 

mag  =  UH60VIBES ( j , [ 1  5  9  13  17  21  25  29  33])'; 

phase  =  UH60VIBES ( j , [2  6  10  14  18  22  26  30  34])'; 

vibsCmplx  =  mag.* (cos (  phase*pi/180  )  +  i*sin(  phase*pi/180  )  ); 

pmag  =  UH60VIBES ( j , [ 3  7  11  15  19  23  27  31  35])'; 

pphase  =  UH60VIBES ( j , [4  8  12  16  20  24  28  32  36])'; 

pvibsCmplx  =  pmag.* (cos (  pphase*pi/180  )  +  i*sin(  pphase*pi/180  ) 

)  ; 

vibset ( : , j )  =  -vibsCmplx+pvibsCmplx; 

end 

%  Select  discrete  vibe  set  to  use  for  comparison 
discrete  vib=vibset ( : , vib  num); 

for  j =1 : length (Coef f) 

adj  ( : , j ) =Coef f { j } \vibset (  : ,  j  )  ; 

%  adj ( : , j ) =Coef f { j } \discrete_vib; 

Cmag ( : , j )  =  abs (adj  ( : , j ) ) ;  %  Calculate  real-number  Mag  and 

Cphase ( : , j )  =  (angle (adj (:, j ))) *180/pi;  %  Phase  values  for  later 

use 

end 

VMEP2AVA.mag  =  Cmag';  %  Stores  Mag  and  Phase  data  to  disk 

VMEP2AVA. phase  =  Cphase'; 
save  VMEP2AVA.mat  VMEP2AVA 

%  Calculate  the  Least  Squares  Method  adjustment 
A  =  LSC.mag. * (cos (  LSC.phase*pi/180  )  +  i*sin(  LSC.phase*pi/180  )  ); 

LS_adj  =A\discrete_vib; 

AVA  adj =LinearUH60 (  discrete  vib (1:8)-  pvibsCmplx ( 1 : 8 )) ; 

%AVA  adj =LinearUH60 (  discrete  vib (1:8)); 

VMEPmag=VMEPadj .mag (vib  num,  : )  ; 

VMEPphase=VMEPadj .phase (vib  num, : ) ; 


%  Determine  the  max  magnitude  so  we  can  normalize  our  adjustments  for 
%  use  in  plotting  on  a  single  graph 
for  k=l:3 

M=max (VMEPmag (k) ,  abs (LS  adj (k))); 

Max ( k) =max (M, AVA  adj.mag(k)); 

end 

subplot (2, 3, r) 
polar (0,1, ' . w ' ) 
hold; 

%  Plot  PCR  adjustments 
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polar (VMEPphase (1) *pi/l 80 , VMEPmag ( 1 ) /Max ( 1 ) , ' *r ' ) 
polar (angle (adj  (1,  z)  )  ,  abs (adj  (1,  z) ) /Max (1) ,  '  ro  '  ) 

%  polar (angle (LS_adj  ( 1 )  )  ,  abs (LS_adj  ( 1 ) ) /Max ( 1 ) ,  ' sr ' ) 

%polar (AVA  adj .phase (1) *pi/180,AVA  adj . mag ( 1 ) /Max ( 1 ) ,  'or') 

%  Plot  Tab  adjustments 

polar (VMEPphase (2) *pi/l 80 , VMEPmag (2 ) /Max (2 )  ,  ' *k ' ) 
polar (angle (adj (2, z) ) , abs (adj (2, z) ) /Max (2) , ' ok' ) 

%  polar (angle (LS_adj (2 ) ) , abs (LS_adj (2 ) ) /Max (2 ) , ' sk ' ) 

%  polar (AVA  adj .phase (2) *pi/180,  AVA  adj .mag (2) /Max (2) ,  'ok') 

%  Plot  Tab  adjustments 

polar (VMEPphase (3) *pi/ 1 80 , VMEPmag ( 3 ) /Max ( 3 ) , ' *b ' ) 
polar (angle (adj  (3,  z)  )  ,  abs (adj  (3,  z) ) /Max (3) ,  '  ob '  ) 

%  polar (angle (LS_adj  (3)  )  ,  abs (LS_adj  (3) ) /Max (3) ,  ' sb ' ) 

%  polar (AVA_adj .phase (3) *pi/ 180,  AVA^adj .mag (3) /Max (3) ,  ' ob ' ) 

title (' Red-wt,  Black-PCR,  Blue-Tab') 

%  figure (2 ) 

%  subplot  (1,3, 1 ) 

%  polar (angle (adj  ( 1 , b : e) ) , abs (adj  ( 1 , b : e) ) ,  ' xr  '  ) 

%  title (' Weight ' ) 

%  subplot (1,3,2) 

%  polar (angle (adj (2 , b : e) ) , abs (adj (2 , b : e) ) , ' xr ' ) 

%  title ('PCR') 

%  subplot  (1,3,3) 

%%  polar (angle (adj  ( 3 , b : e) ) , abs (adj  ( 3 , b : e) ) ,  '  xr  '  ) 

%  title ( 'Trim  Tab' ) 

end 

VMEPcplx  =  VMEPadj .mag. * (cos (VMEPadj .phase*pi/180)  +  ... 

i*sin (VMEPadj .phase*pi/180)  )  ; 
for  j=l : 20 

diff ( j ,  : ) =abs (VMEPcplx ( j  ,  :  )  -  (A\vibset ( :  ,  j ) )  .  ' ) 

end 

max (diff) 
figure (2 ) 

polar (AVA  adj .phase (2) *pi/180, AVA  adj .mag (2 ) , ' or ' ) 
hold  on 

polar (VMEPphase (2) *pi/l 80 , VMEPmag (2 )  ,  ' *r ' ) 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 


%  Capt  Miller 

%  This  function  creates  a  single  coefficient  matrix  from  N  coefficient 
%  matrices  by  summing  the  squares  of  the  N  coefficient  matrix 
components 

%  and  then  taking  the  square  root.  It  also  calculates  the  standard 
%  deviation  of  all  ad  hoc  coefficients  from  their  mean  value. 

function  [LSC,  standard] =Least_Squares_Coef f (Coeff) 

numcoef f=length (Coef f ) ; 
numrows=length (Coef f  { 1 }  ( :  ,  1 )  )  ; 
numcollumns=length (Coef f  { 1 }  ( 1 ,  : )  )  ; 
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C=zeros (numrows , numcollumns ) ; 


for  j=l:numcoeff 

A.mag=abs (Coef f {  j  } )  ; 

A.phase=angle (Coeff  {  j  } )  ; 

C=C+ (A. mag. A2)  . * (cos (  A. phase  )  +  i*sin(  A. phase  )  ); 

end 

LSC.mag=sqrt (abs (C) /numcoeff)  ; 

LSC . phase=angle (C) *180/pi; 
for  j  =1 : numrows 

for  k=l : numcollumns 

if  LSC . phase ( j , k)  <  0 

LSC . phase ( j , k)  =  LSC. phase (j , k) +360; 

end 

end 

end 

LSCcmplx  =  LSC. mag. * (cos ( LSC. phase *pi/ 180) +i*sin ( LSC. phase *pi/ 180) ) ; 
for  j=l:numcoeff 

dif f { j } =abs (Coeff { j } -LSCcmplx) ; 

end 

for  j=l: numrows 

for  k=l : numcollumns 
for  1=1: numcoeff 

temp ( 1 ) =dif f { 1 } ( j , k) ; 

end 

standard ( j , k)  =  std(temp); 

end 

end 


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

%  Capt  Miller 

%  This  program  will  maesure  the  difference  between  the  VMEP  adjustment 
and 

%  the  ad  hoc  adjustment  for  each  flight.  Intent  is  to  determine 
%  if  the  minus  zero  produces  any  differences  in  adjustments. 

clear; clc; close  all; 

load  Coeff. mat 
load  UH60VIBES .mat 
load  UH60  Coeff  Data. mat 
load  VMEPadj.mat 


n=2  0 

[LSC,  standard] =Least_Squares_Coeff (Coeff) 
for  j  =1 : n 

mag  =  UH60VIBES ( j , [ 1  5  9  13  17  21  25  29  33])'; 
phase  =  UH60VIBES ( j , [2  6  10  14  18  22  26  30  34])'; 
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vibsCmplx  =  mag.* (cos (  phase*pi/180  )  +  i*sin(  phase*pi/180  )  ); 

pmag  =  UH60VIBES ( j , [ 3  7  11  15  19  23  27  31  35])'; 

pphase  =  UH60VIBES ( j , [4  8  12  16  20  24  28  32  36])'; 

pvibsCmplx  =  pmag.* (cos (  pphase*pi/180  )  +  i*sin(  pphase*pi/180  ) 


vibset ( : , j )  =  -vibsCmplx+pvibsCmplx; 

end 

VMEPcplx  =  VMEPadj .mag. * (cos (VMEPadj .phase*pi/180)  +  ... 
i*sin (VMEPadj .phase*pi/180)  )  ; 

A  =  LSC.mag. * (cos (  LSC.phase*pi/180  )  +  i*sin(  LSC.phase*pi/180  )  ); 

for  j=l : 10 

tmp { 2* j -1 , : } =Coef f { 2* j } ; %  tmp(odd)  =  A (even) 
tmp { 2 * j ,  : } =Coef f { 2* j -1 } ; %  tmp (even)  =  A (odd) 

end 

%clear  Coeff;  %  Do  not  rem  out  if  you  wish  to  swap  around 
%Coeff=tmp;  %  ad  hoc  matrices 

for  j  =1 : n 

%adj ( j , : ) =A\vibset ( : , j ) ;  %  This  one  for  plotting  RMS  ad  hoc 
adj ( j , : ) =Coef f { j } \vibset ( : , j ) ; %This  one  for  plotting  individual  ad 
hoes 
end 

for  type=l:3 

%f igure (type) 
subplot (1,3, type) 

polar (angle (adj  ( 1 : n, type) ) , abs (adj  ( 1 : n, type) ) ,  '  ko  '  ) 
hold  on 

polar (angle (VMEPcplx ( 1 : n, type) ) , abs (VMEPcplx ( 1 : n, type) ) , ' k+ ' ) 
end 

for  j=l : 20 

diff ( j , : ) =abs (VMEPcplx ( j , : )  -  adj ( j , : ) ) ; 

end 

diff 

max (diff) 

%f igure ( 1 ) 
subplot (1,3,1) 
title ( ' Weight ' ) 

%legend('Ad  Hoc  Adj ustment ',' Real  Adjustment') 

%f igure (2 ) 
subplot  (1,3,2) 
title (' Pitch  Link') 

%legend('Ad  Hoc  Adj ustment ',' Real  Adjustment') 

%f igure ( 3 ) 
subplot (1,3,3) 
title ( 'Trim  Tab' ) 

%legend('Ad  Hoc  Adj ustment ',' Real  Adjustment') 


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 


E-9 


%  Capt  Nathan  A  Miller 

%  Determining  difference  in  AVA  and  VMEP  methods  using  Histograms. 

clear;  clc; 
load  AVAadj.mat 
load  VMEPadj.mat 

AVAcomplex  =  AVAadj .mag . * (cos (  AVAadj .phase*pi/180  )  +... 

i*sin(  AVAadj .phase*pi/180  )  ); 

VMEPcomplex  =  VMEPadj .mag. * (cos (  VMEPadj .phase*pi/180  )  +... 

i*sin(  VMEPadj .phase*pi/180  )  ); 

Magnitude  diff  =  abs (AVAcomplex  -  VMEPcomplex) ; 

Rel  mag  diff  =  Magnitude  diff . /AVAadj .mag; 

for  i=l : 20 

if  abs (AVAcomplex ( i , 2 ) )  <  1/2 

%  if  abs (VMEPcomplex ( i , 2 ) )  <  113/2 
Rel  mag  diff(i,2)=0; 

abs (AVAcomplex ( i , 2 ) ) 
abs (VMEPcomplex ( i , 2 ) ) 

%end 

end 

if  abs (AVAcomplex ( i , 3 ) )  <  2/2 

%  if  abs (VMEPcomplex ( i , 3 ) )  <  0.5/2 
Rel_mag_dif f ( i , 3 ) =0 ; 

abs (AVAcomplex ( i , 3 ) ) 
abs (VMEPcomplex (i, 3) ) 

%end 

end 

end 

X  =  ( [012345678  9]+0.5)/10; 
for  i=l:2 

[n ( i , : ) , xout ( i , : ) ] =hist (Rel_mag_dif f (  : , i  +  1 ) , X) ; 
subplot ( 1 , 2 , i ) ; 

bar(xout(i,  : ) *  1 0  0 , n ( i ,  : ) / 0 . 2 ,  ' k '  )  ; 

AXIS ( [0  100  0  50] ) 

ylabel (' Percentage  of  Total  Flights  Sampled') 

end 

subplot  (1,2, 1 ) 
title (' Pitch  Link') 

subplot  (1,2,2) 
title ( 'Trim  Tab' ) 

text (-60, -3, 'Percent  Difference  of  PC-GBS  Adjustments  from  AVA 
Algorithm  Adjustments') 


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
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%  Capt  Miller  Thesis 

%  This  is  a  test  program  to  determine  if  any  mistakes  were  made 
%  during  data  entry.  Any  errors  in  vib  data  entry  will  appear  as  gross 
%  differences  in  the  plotted  ad  hoc  coefficients. 

clear  all;  elf;  clc;close  all 
load  Coeff.mat 

[LSC,  standard] =Least_Squares_Coef f (Coeff) ; 

[LSCgood,  standardgood] =Least_Squares_Coef f (Coef f (1:5) ) ; 

[LSCexd,  standardexd] =Least_Squares_Coef f (Coef f (15:20) ) ; 
b=l  ; 

e=length (Coeff )  ; 
n=l; 

FLG=-1  %positive  for  single  fig 

state=6; 

for  h=l : 9 

for  j=b:e 

value (h, j  +l-b,  : ) =Coef f { j  }  (h,  : )  ; 

end 

end 

latvalue  =  value ([1  2  3  5  7] ,  : ,  :  )  ; 
vertvalue  =  value ( [9  4  6  8]  ,  : ,  : )  ; 

for  h=l : 9 

maxi (h) =max (abs (value (h,  : ,  1 )  )  )  ; 
max2 (h) =max (abs (value  (h,  : ,  2 )  )  )  ; 
max3 (h) =max (abs (value  (h,  : , 3)  )  )  ; 

end 

for  h=l:5 

latmaxl (h) =max (abs (latvalue (h,  : , 1)  )  )  ; 
latmax2 (h) =max (abs (latvalue (h,  : , 2)  )  )  ; 
latmax3 (h) =max (abs (latvalue (h, :  , 3)  )  )  ; 

end 

for  h=l:4 

vertmaxl (h) =max (abs (vertvalue (h, : ,  1 )  )  )  ; 
vertmax2 (h) =max (abs (vertvalue (h, : ,  2 )  )  )  ; 
vertmax3 (h) =max (abs (vertvalue (h, : ,  3)  )  )  ; 

end 

for  h=l : 9 
%figure (h) 
subplot ( 3 , 3 , h) 

%polar (angle (value (h, : , 1 ) ) , abs (value (h, : , 1 ) ) /maxi (h) , ' bx ' ) 

%hold  on 

polar (angle (value (h,  : ,  n)  )  ,  abs (value (h,  : , n) ) ,  '  rx  '  ) 

%f igure (2 ) 

%subplot ( 3 , 3 , h) 

%polar (angle (value (h,  : , 2 ) ) , abs (value (h,  : , 2 ) ) /max2 (h) ,  '  rx '  ) 

%f igure ( 3 ) 

%subplot ( 3 , 3 , h) 

%polar (angle (value (h, : , 3) ) , abs (value (h, : , 3) ) /max3 (h) , ' gx ' ) 
end 

subplot (3,3,1) 
title ( ' FPG1 00  A-B' ) 
subplot  (3,3,2) 
title ('Hover  A-B') 
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subplot (3,3,3) 
title (' 80  Kts  A-B') 
subplot (3,3,4) 
title (' 80  Kts  A+B ' ) 
subplot  (3,3,5) 
title ( ' 120  Kts  A-B' ) 
subplot (3,3,6) 
title ( ' 120  Kts  A+B' ) 
subplot (3,3,7) 
title ('145  Kts  A-B') 
subplot (3,3,8) 
title ('145  Kts  A+B') 
subplot (3,3,9) 
title ('Hover  A+B') 
if  FLG>0 
figure  ( 1 ) 
elf ; 

polar (angle (value (state, : , n) ) , abs (value (state, : , n) ) , ' rx ' ) 
end 

figure  (2 ) 
for  i=l:5 

temp  =  max (latmaxl (i) *10, latmax2 (i) ) ; 
maxval  =  max (temp, latmax3 (i) ) ; 

subplot (2, 3, i)  %  Each  subplot  is  for  a  different  lateral  regime 
polar ( 0 , maxval , ' w . ' ) 
hold  on 

polar (angle ( latvalue ( i ,  : , 1 ) ) , abs ( latvalue (i,  : ,  1 ) ) *  10 ,  ' kx  '  ) 
polar (angle (latvalue (i, : , 2) ) , abs (latvalue ( i , : , 2 ) ) , ' kd ' ) 
polar (angle (latvalue (i,  : , 3) ) , abs (latvalue (i,  : , 3) ) ,  ' kA ' ) 
end 

figure ( 3 ) 
for  i=l : 4 

temp  =  max (vertmaxl ( i ) , vertmax2 ( i ) ) ; 
maxval  =  max (temp, vertmax3 (i) )  ; 

subplot (2, 2, i)  %  Each  subplot  is  for  a  different  lateral  regime 
polar ( 0 , maxval , ' w . ' ) 
hold  on 

polar (angle (vertvalue ( i , : , 1 ) ) , abs (vertvalue (i, : , 1 ) ) *10 , 'kx ')%%%%%%% 

S-S-S-S-S-S-S-S-S-S-S-S-S-S-S-S-S-S-S- 

ooooooooooooooooooo 

polar (angle (vertvalue (i,  : , 2) ) , abs (vertvalue ( i ,  : ,  2 ) )  ,  '  kd '  ) 
polar (angle (vertvalue (i, : , 3) ) , abs (vertvalue ( i , : , 3 ) ) , ' kA ' ) 
end 

abs (LSCexd.mag-LSCgood.mag) . /LSCgood.mag*100 

figure  ( 4 ) 
i=5 

polar  (0,  . 06,  ' w. ' ) 
hold  on 

polar (angle (latvalue ( i ,  : , 1 ) ) , abs (latvalue (i,  : , 1 ) ) *  10 ,  'kx') 
polar (angle (latvalue (i,  : , 2) ) , abs (latvalue (i,  : , 2) )  ,  ' kd' ) 
polar (angle (latvalue (i,  : , 3) ) , abs (latvalue  (i,  : , 3) ) ,  '  kA  '  ) 
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Appendix  F:  AH-64A  Matlab  Code 


Note:  The  Matlab  codes  contained  in  this  appendix  were  each  used  to  perform  one  or 
more  specific  functions  in  the  course  of  this  research.  There  are  several  instances  where 
parts  of  the  code  have  been  commented  out  with  the  %  symbol.  In  order  to  reproduce  all 
of  the  analysis  of  this  thesis,  some  lines  of  code  may  need  to  be  un-commented  and  other 
lines  commented  out. 

%  Capt  Nathan  A  Miller 

%  This  program  creates  MAT  data  files  from  imported  Excel  database 
%  Ensure  that  the  two  databases  are  imported  into  matlab  prior  to 
running 

function  Excel  to  MAT 
global  data  data2 
n  =  length (data2 ( : , 2 ) ) -2 

%  Read  Vibration  Magnitude  and  Phase  as  well  as  prediction  magnitude 
%  and  phase  from  Excel  File  and  assign  to  data  matrix 
AH64A  Vibes  =  zeros (n, 28); 
for  j  =1 : n 

for  k=l : 7 


AH64A_Vibes ( j  , 1+4* (k-1) ) 
AH64A_Vibes ( j  , 2+4* (k-1) ) 
AH64A_Vibes ( j  ,  3  +  4* (k-1)  ) 
AH64A_Vibes ( j  ,  4  +  4* (k-1)  ) 

end 


=  data ( l  +  k+ ( j -1 ) *11 , 2 )  ; 
=  data ( l+k+ ( j -1 ) *11 , 3 ) ; 
=  data (l  +  k+ ( j-1) *11,4)  ; 
=  data ( l+k+ ( j -1 ) *11,5); 


end 

save  AH64A  Vibes. mat  AH64A  Vibes 
AH64A  VMEP  Adjustments  =  zeros (n, 20); 
for  j  =1 : n 

for  k=l : 20 


AH64A  VMEP  Ad j ustments ( j , k)  =  data2 ( j +2 , k+1 ) ; 
end 


end 

save  AH64A  VMEP  Adjustments .mat  AH64A  VMEP  Adjustments 
AH64A  Coeff  Data  =  zeros (n, 75); 
for  j  =1 : n 

AH64A_Coef f_Data ( j , 1 ) =data (2+ (j— 1)  *11,  6)  ; 

AH64A_Coeff_Data ( j , 16) =data (2+ ( j-1) *11, 9) ;  %  Records  manual  adjustment 

AH64A_Coef f_Data ( j , 31 ) =data (2+ ( j -1 ) *11 , 12 ) ;  %  that  produces  delta  vibes 
AH64A_Coef f_Data ( j , 4  6) =data (2+(j-l)*ll,15); 

AH64A_Coef f_Data ( j , 61 ) =data (2+(j— 1)  *11,18)  ; 
for  k=l : 7 


for  z=l : 5 

AH64A_Coef f_Data ( j , 2*k+15* (z-1) )  =  data (l+k+ ( j-1) *11, 7  +  3*  (z-1) ) ; 
AH64A_Coef f_Data ( j , 2*k+l+15* (z-1) ) =data (l+k+ ( j-1) *11, 8+3*  (z-1) ) ; 
end 


end 


end 


save  AH64A  Coeff  Data. mat  AH64A  Coeff  Data 


oooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooo 
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%  Capt  Nathan  A  Miller 
%  Thesis  work 

%  This  program  loads  vibe  data  and  calculates  an  AVA  adjustment  matrix 

function  Vibs2AVA  ADJ 

load  AH64A  Vibes. mat 
n  =  length (AH64A  Vibes (:,  1)  ); 

for  j  =1 : n 

mag  =  AH64A  Vibes (j,[l  5  9  13  17  21  25])'; 
phase  =  AH64A  Vibes ( j ,  [2  6  10  14  18  22  26]  )  '  ; 

vibsCmplx  =  mag.* (cos (  phase*pi/180  )  +  i*sin(  phase*pi/180  )  ); 

pmag  =  AH64A  Vibes ( j ,  [3  7  11  15  19  23  27] )  '  ; 
pphase  =  AH64A  Vibes ( j ,  [4  8  12  16  20  24  28] )  '  ; 

pvibsCmplx  =  pmag.* (cos (  pphase*pi/180  )  +  i*sin(  pphase*pi/180  )  ); 

delvib  =  -vibsCmplx  +  pvibsCmplx; 
adj(j)  =  LinearAH64A(  delvib  ); 

end 

for  k=l : n 

magnitude ( k, : )  =  adj(k).mag'; 

Phase (k, :)  =  adj ( k) .phase'; 


end 

AVAadj.mag  =  magnitude; 
AVAadj. phase  =  Phase; 
save  AVAadj. mat  AVAadj 


S-S-S-S-S-S-S-S-S-S-S-S-S-S-S-S-S-S-S-S-S-S-S-S-S-S-S-S-S-S-S-S-S-S-S-S-S-S-S-S-S-S-S-S-S-S-S-S-S-S-S-S-S-S-S-S-S-S-S-S-S-S-S-S-S-S-S- 

ooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooo 


%  Capt  Nathan  A  Miller 

%  Converts  VMEP  adjustments  into  magnitude  and  phase  of  adjustments. 

function  VMEP  2  MagPhase 

load  AH64A  VMEP  Adjustments .mat 
n  =  length (AH64A  VMEP  Adjustments (:, 1) ) ; 

a=AH64A  VMEP  Adjustments; 
clear  AH64A  VMEP  Adjustments 

for  j  =1 : n 


blade  #  1 


adj  cmplx ( j , 1 ) 

=  +a ( j , 1) -i 

adj  cmplx ( j , 2 ) 

=  +a (j , 2) -i 

adj  cmplx ( j , 3 ) 

=  +a ( j , 3) -i 

adj  cmplx ( j , 4 ) 

=  +a ( j , 4) -i 

adj  cmplx ( j , 5 ) 

=  +a ( j , 5) -i 

2  3  4 

*a(j,  6)  -a(j,  11)  +i*a  (j  ,  16)  ; 
*a(j,7)-a(j,12)+i*a(j,17) ; 
*a(j, 8) -a(j, 13) +i*a ( j , 18) ; 
*a(j, 9) -a(j, 14) +i*a (j  ,  19)  ; 
*a(j, 10) -a(j, 15) +i*a(j,20) 


%  Weight  Adj 
%  P/L  Adj 
%  Tab  Adj  8 
%  Tab  Adj  6 
%  Tab  Adj  4 


end 
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VMEPadj .mag=abs (adjcmplx) ; 

VMEPadj . phase  =  angle (adjcmplx) *180/pi; 

save  VMEPadj. mat  VMEPadj 


ooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooo 


%  Capt  Miller 

%  This  program  generates  plots  of  the  ad  hoc  adjustments 

clear; clc; close  all; 

load  Coeff.mat 
load  AH64A  Vibes. mat 

for  j =1 : length (Coeff) 

mag  =  AH64A  Vibes ( j ,  [1  5  9  13  17  21  25] )  '  ; 
phase  =  AH64A  Vibes ( j ,  [2  6  10  14  18  22  26]  )  1  ; 

vibsCmplx  =  mag.* (cos (  phase*pi/180  )  +  i*sin(  phase*pi/180  )  ); 

pmag  =  AH64A  Vibes ( j ,  [3  7  11  15  19  23  27] )  '  ; 
pphase  =  AH64A  Vibes ( j ,  [4  8  12  16  20  24  28] ) ' ; 

pvibsCmplx  =  pmag.* (cos (  pphase*pi/180  )  +  i*sin(  pphase*pi/180  ) 

)  ; 


vib ( : , j )  =  vibsCmplx; %  -  pvibsCmplx; 


end 

v=vib ( : , 1 ) 

for  j =1 : length (Coeff ) 

adj ( : , j ) =Coef f { j } \v; 

magnitude ( : , j )  =  abs (adj  ( : ,  j ) )  ; 

Faze ( : , j )  =  (angle (adj ( : , j ) ) ) *180/pi; 


end 


ava=LinearAH64A (  v  ); 

avaadj =ava . mag . * (cos (  ava.phase*pi/180  )  +  i*sin(  ava.phase*pi/180 


VMEP2AVA.mag  =  magnitude'; 

VMEP2AVA. phase  =  Faze'; 
save  VMEP2AVA.mat  VMEP2AVA 

figure ( 1 ) 

if  max (abs (ava .mag ( 1 ) ) )  >  max (abs (adj ( 1 , : ) ) ) 
polar (ava . phase (1) *pi/180, ava . mag ( 1 ) , ' *r ' ) 
hold 

polar (angle (adj  ( 1 ,  : ) ) , abs (adj  ( 1 ,  : ) ) ) 
else 

polar (angle (adj ( 1 , : ) ) , abs (adj ( 1 , : ) ) ) 
hold 

polar (ava . phase (1) *pi/l 80 , ava .mag ( 1 ) , ' *r ' ) 
end 
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ava .mag 
ava . phase 
figure  (2 ) 

if  max (abs (ava .mag (2 ) ) )  >  max (abs (adj (2 , : ) ) ) 
polar (ava . phase (2) *pi/l 80 , ava . mag (2 ) , ' *r ' ) 
hold 

polar (angle (adj  (2 , : ) ) , abs (adj  (2 ,  : )  )  ) 
else 

polar (angle (adj  (2 ,  : ) ) , abs (adj  (2 ,  : )  )  ) 
hold 

polar (ava . phase (2) *pi/l 80 , ava . mag (2 )  ,  ' *r '  ) 
end 

figure  ( 3 ) 

if  max  (abs  (ava. mag  (3)  )  )  >  max  (abs  (adj  (3,  : )  )  ) 
polar (ava . phase (3) *pi/l 80 , ava . mag ( 3 )  ,  ' *r ' ) 
hold 

polar (angle (adj  (3,  : ) ) , abs (adj  (3,  : ) ) ) 
else 

polar (angle (adj  (3,  : ) ) , abs (adj  (3,  : ) ) ) 
hold 

polar (ava . phase (3) *pi/l 80 , ava . mag ( 3 )  ,  ' *r '  ) 
end 

figure  ( 4 ) 

if  max  (abs  (ava. mag  (4)  )  )  >  max  (abs  (adj  (4,  : )  )  ) 
polar (ava . phase (4) *pi/l 80 , ava . mag ( 4 )  ,  ' *r ' ) 
hold 

polar (angle (adj (4, : ) ) , abs (adj (4, : ) ) ) 
else 

polar (angle (adj (4, : ) ) , abs (adj (4, : ) ) ) 
hold 

polar (ava . phase (4) *pi/l 80 , ava . mag ( 4 )  ,  '  *r  '  ) 
end 

figure  ( 5 ) 

if  max  (abs  (ava. mag  (5)  )  )  >  max  (abs  (adj  (5,  : )  )  ) 
polar (ava . phase (5) *pi/l 80 , ava . mag ( 5 )  ,  ' *r ' ) 
hold 

polar (angle (adj (5, : ) ) , abs (adj (5, : ) ) ) 
else 

polar (angle (adj (5, : ) ) , abs (adj (5, : ) ) ) 
hold 

polar (ava . phase (5) *pi/l 80 ,  ava . mag ( 5 )  ,  ' *r  '  ) 
end 


%  Capt  Miller 

%  This  program  converts  the  coefficient  matrix  data  into  actual 
%  coefficient  matrices.  The  program  loads  the  raw  vibe  data  from 
%  AH64A  Vibes. mat  and  computes  a  delta  vibe  from  the  predicted  vibe 
%  values  stored  in  AH64A  Coeff  Data. mat.  The  coefficient  matrix  is 
%  computed  column  by  column.  Once  computed,  the  coefficient  matrix 
%  operates  on  a  vibe  set  to  produce  a  quasi-linear  adjustment  set. 

%  This  process  is  repeated  for  each  coefficient  data  set.  The  same 
%  vibe  set  is  always  used  for  comparison. 

function  VMEP  Coefficient  Matrix 

load  AH64A  Vibes. mat 
load  AH64A  Coeff  Data. mat 
n=length (AH64A  Vibes ( : , 1 ) ) ; 
for  j  =1 : n 

wt  =  AH64A  Coeff  Data(j,l); 
pi  =  AH64A_Coeff_Data ( j  ,  16)  ; 
tab8  =  AH64A_Coeff_Data ( j  ,  31) ; 
tab6  =  AH64A_Coeff_Data ( j  ,  46)  ; 
tab4  =  AH64A_Coef f_Data ( j  ,  61 )  ; 

for  k=l : 7  %  loops  over  the  7  flight  conditons,  FPG100  through  140  Knots 
vibmag=AH64A  Vibes ( j , 1+4* (k-1) ) ; 
vibphase=AH64A  Vibes ( j , 2  +  4* (k-1) )  ; 
vibcmplx=vibmag* (cos (  vibphase*pi/180  )  +  ... 
i*sin(  vibphase*pi/180  )  ); 

pvibwtmag=AH64A  Coeff  Data (j  ,  2+2* (k-1) )  ; 
pvibwtphase=AH64A  Coeff  Data (j  ,  3+2* (k-1)  )  ; 
pvibwtcmplx=pvibwtmag* (cos (  pvibwtphase*pi/180  )  +  ... 
i*sin(  pvibwtphase*pi/180  )  ); 

pvibplmag=AH64A  Coeff  Data (j  ,  17+2* (k-1)  )  ; 
pvibplphase=AH64A  Coeff  Data (j  ,  18+2* (k-1)  )  ; 
pvibplcmplx=pvibplmag* (cos (  pvibplphase*pi/180  )  +  ... 
i*sin(  pvibplphase*pi/180  )  ); 

pvibtab8mag=AH64A  Coeff  Data (j  ,  32+2* (k-1)  )  ; 
pvibtab8phase=AH64A  Coeff  Data (j , 33+2* (k-1) )  ; 
pvibtab8cmplx=pvibtab8mag* (cos (  pvibtab8phase*pi/180  )  +  ... 
i*sin(  pvibtab8phase*pi/180  )  ); 

pvibtab6mag=AH64A  Coeff  Data (j , 47+2* (k-1) ) ; 
pvibtab6phase=AH64A  Coeff  Data (j  ,  48+2* (k-1)  )  ; 
pvibtab6cmplx=pvibtab6mag* (cos (  pvibtab6phase*pi/180  )  +  ... 
i*sin(  pvibtab6phase*pi/180  )  ); 

pvibtab4mag=AH64A  Coeff  Data (j , 62+2* (k-1) ) ; 
pvibtab4phase=AH64A  Coeff  Data (j  ,  63+2* (k-1)  )  ; 
pvibtab4cmplx=pvibtab4mag* (cos (  pvibtab4phase*pi/180  )  +  ... 
i*sin(  pvibtab4phase*pi/180  )  ); 

dvibwt=pvibwtcmplx-vibcmplx; 
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dvibpl=pvibplcmplx-vibcmplx; 

dvibtab8=pvibtab8cmplx-vibcmplx; 

dvibtab6=pvibtab6cmplx-vibcmplx; 

dvibtab4=pvibtab4cmplx-vibcmplx; 

A ( k, 1 ) =dvibwt/wt; 

A ( k, 2 ) =dvibpl/pl ; 

A ( k, 3 ) =dvibtab8 / tab8 ; 

A ( k, 4 ) =dvibtab6/tab6; 

A ( k, 5 ) =dvibtab4 / tab4 ; 

end 

Coef f { j } =A; 

end 

save  Coeff.mat  Coeff 

oooooooooooooooooooooooooooooooooooooooooooooooooooooooooooooo 

%  Capt  Miller  Thesis 

%  This  is  a  test  program  to  determine  if  any  mistakes  were  made 

%  during  data  entry. 

clear;  elf ;  clc; 

load  Coeff.mat 

b=l 

e=length (Coeff) 
n=l 

for  h=l : 7 

for  j=b:e 

value (h, j  +  l-b,  : ) =Coef f { j  }  (h,  : )  ; 

end 

end 

for  h=l : 7 

maxi (h) =max (abs (value (h,  : , 1) ) )  ; 
max2 (h) =max (abs (value  (h,  : , 2 ) ) )  ; 
max3 (h) =max (abs (value  (h,  : , 3) ) )  ; 
max4 (h) =max (abs (value (h,  : , 4 ) ) )  ; 
max5 (h) =max (abs (value  (h,  : , 5)  )  )  ; 

end 

for  h=l : 7 
%figure (h) 
subplot ( 3 , 3 , h) 

%polar (angle (value (h, : , 1 ) ) , abs (value (h, : , 1 ) ) /maxi (h) , ' bx ' ) 
polar (angle (value (h,  : ,  n)  )  ,  abs (value (h,  : , n) ) ,  '  rx  1  ) 

%hold  on 
%f igure  (2 ) 

%subplot ( 3 , 3 , h) 

%polar (angle (value (h, : , 2 ) ) , abs (value (h, : , 2 ) ) /max2 (h) , ' rx ' ) 

%f igure  ( 3 ) 

%subplot ( 3 , 3 , h) 

%polar (angle (value (h, : , 3) ) , abs (value (h, : , 3) ) /max3 (h) , ' gx ' ) 
%polar (angle (value (h, : , 4 ) ) , abs (value (h, : , 4 ) ) /max4 (h) , ' mx ' ) 
%polar (angle (value (h, : , 5) ) , abs (value (h, : , 5) ) /max5 (h) , ' cx ' ) 
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end 

subplot (3,3, 1 ) 
title (' FPG100  Lat') 
subplot (3,3,2) 
title ('Hover  Lat') 
subplot (3,3,3) 
title (' 60  Kts  Vert') 
subplot (3,3,4) 
title (' 80  Kts  Vert') 
subplot (3,3,5) 
title ('100  Kts  Vert') 
subplot (3,3,6) 
title ('120  Kts  Vert') 
subplot (3,3,7) 
title ('140  Kts  Vert') 

%f igure  ( 1 ) 

%clf  ; 

%polar (angle (value ( 4 ,  : , 1 ) ) , abs (value ( 4 ,  : ,  1 ) )  ,  ' bx ' ) 
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%  Capt  Miller 

%  Rough  plots  of  coeff  matrix  results 

clear; clc; elf ; 

load  Coeff. mat 
load  AH64A  Vibes. mat 
load  VMEPadj.mat 

[LSC,  standard] =Least_Squares_Coeff (Coeff ) ; 

%  Runs  the  Least  Squares  Coefficient  function  to 
%  generate  the  LSC  matrix  from  the  Coeff  cell, 
vib  num=20 
b=l_ 

e=length (Coeff) 

for  j =1 : length (Coeff ) 

mag  =  AH64A  Vibes ( j , [1  5  9  13  17  21  25] ) ' ; 
phase  =  AH64A  Vibes ( j , [2  6  10  14  18  22  26] ) ' ; 

vibsCmplx  =  mag.* (cos (  phase*pi/180  )  +  i*sin(  phase*pi/180  )  ); 

pmag  =  AH64A  Vibes (j,[3  7  11  15  19  23  27])'; 
pphase  =  AH64A  Vibes ( j ,  [4  8  12  16  20  24  28] )  '  ; 

pvibsCmplx  =  pmag.* (cos (  pphase*pi/180  )  +  i*sin(  pphase*pi/180  ) 

)  ; 

vibset ( : , j )  =  -vibsCmplx  +  pvibsCmplx; 

end 

discrete  vib=vibset ( : , vib  num) ;  %  Select  discrete  vibe  set  to  use  for 
comparison 

for  j =1 : length (Coeff ) 

adj ( : , j ) =Coef f { j } \discrete_vib; 

Cmag ( : , j )  =  abs ( adj  ( : ,  j ) )  ; 

Cphase ( : , j )  =  (angle (adj ( : , j ) ) ) *180/pi; 

end 

VMEP2AVA.mag  =  Cmag'; 

VMEP2AVA. phase  =  Cphase'; 
save  VMEP2AVA.mat  VMEP2AVA 
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%  Calculate  the  Least  Squares  Method  adjustment 
A  =  LSC.mag. * (cos (  LSC.phase*pi/180  )  +  i*sin(  LSC.phase*pi/180  )  ); 

LS_adj  =A\discrete_vib; 

AVA  adj =LinearAH64A (  discrete_vib  ); 

VMEPmag=VMEPadj .mag (vib  num,  : ) 

VMEPphase=VMEPadj .phase (vib  num, : ) 

%  Determine  the  max  magnitude  so  we  can  normalize  our  adjustments  for 
%  use  in  plotting  on  a  single  graph 
for  k=l : 5 

Maxl=max (max (Cmag  ( k, b :  e) )  ,  AVA  adj. mag (k) ) ; 

%Maxl=max (abs (AVA  adj .mag (k) ) ) ; 

Max2=max (VMEPmag (k) ,  abs (LS  adj (k) ) ) ; 

Max ( k) =max (Maxi , Max2 ) ; 

%Max ( k ) =max ( Cmag ( k , b : e ) ) ; 

end 

%f igure ( 1 ) 
subplot (2,3,1) 
polar  (0,1,  ' . w ' ) 
hold 

%  Plot  weight  adjustments 

%  polar (Cphase ( 1 , b : e) *pi/l 80 , Cmag ( 1 , b : e) /Max ( 1 )  ,  ' xr  '  ) 

polar (angle (LS_adj ( 1 ) ) , abs (LS_adj ( 1 ) ) /Max ( 1 ) , ' sr ' ) 
polar (AVA  adj .phase (1) *pi/ 180, AVA  adj . mag ( 1 ) /Max ( 1 ) , ' or ' ) 
polar (VMEPphase (1) *pi/l 80 , VMEPmag ( 1 ) /Max ( 1 ) , ' *r ' ) 

%  Plot  PCR  adjustments 

%  polar (Cphase (2 , b : e) *pi/l 80 , Cmag (2 , b : e) /Max (2 ) , ' xk ' ) 

polar (angle (LS_adj (2 ) ) , abs (LS_adj (2 ) ) /Max (2 ) , ' sk ' ) 
polar (AVA  adj .phase (2) *pi/180, AVA  adj .mag (2 ) /Max (2 ) , ' ok ' ) 
polar (VMEPphase (2) *pi/l 80 , VMEPmag (2 ) /Max (2 ) , ' *k ' ) 

%  Plot  tab  8  adjustments 

%  polar (Cphase ( 3 , b : e) *pi/180, Cmag ( 3 , b : e) /Max ( 3 ) , ' xb ' ) 

polar (angle (LS_adj  (3) ) , abs (LS_adj  (3) ) /Max (3) ,  '  sb '  ) 
polar (AVA_adj . phase (3) *pi/l 80 , AVA_adj . mag ( 3 ) /Max ( 3 ) , ' ob ' ) 
polar (VMEPphase (3) *pi/l 80 , VMEPmag ( 3 ) /Max ( 3 ) , ' *b ' ) 

%  Plot  tab  6  adjustments 

%  polar (Cphase ( 4 , b : e) *pi/l 80 , Cmag ( 4 , b : e) /Max ( 4 ) , ' xg ' ) 

polar (angle (LS_adj  ( 4 ) ) , abs (LS_adj  ( 4 ) ) /Max ( 4 ) ,  ' sg  '  ) 
polar (AVA  adj .phase (4) *pi/180, AVA  adj . mag ( 4 ) /Max ( 4 ) , ' og ' ) 
polar (VMEPphase (4) *pi/l 80 , VMEPmag ( 4 ) /Max ( 4 )  ,  ' *g ' ) 

%  Plot  tab  4  adjustments 

%  polar (Cphase ( 5 , b : e) *pi/l 80 , Cmag ( 5 , b : e) /Max ( 5 ) , ' xm ' ) 

polar (angle (LS_adj  (5) ) , abs (LS_adj  (5) ) /Max (5) ,  ' sm'  ) 
polar (AVA_adj .phase (5) *pi/180, AVA_adj .mag (5) /Max (5) , ' om' ) 
polar (VMEPphase (5) *pi/l 80 , VMEPmag ( 5 ) /Max ( 5 ) , ' *m ' ) 

title (' Red-Wt,  Black-Pcr,  Blue-Tab8,  Green-Tab6,  Magenta-Tab4 ' ) 
subplot  (2,3,2) 

polar (Cphase ( 1 , b : e) *pi/180, Cmag ( 1 , b : e) , ' xr ' ) 
title ( ' Weight ' ) 
subplot (2,3,3) 

polar (Cphase (2 , b : e) *pi/l 80 , Cmag (2 , b : e) , 'xk' ) 
title (' Pitch  Link') 


F-8 


subplot  (2,3,4) 

polar (Cphase ( 3 , b : e) *pi/180, Cmag ( 3 , b : e) , ' xb ' ) 
title ( ' Tab  8-10  '  ) 
subplot (2,3,5) 

polar (Cphase ( 4 , b : e) *pi/180, Cmag ( 4 , b : e) , ' xg ' ) 
title ( ' Tab  6-10  '  ) 
subplot (2,3,6) 

polar (Cphase ( 5 , b : e) *pi/180, Cmag ( 5 , b : e) , ' xm ' ) 
title ('Tab  4-10') 

VMEPcplx  =  VMEPadj .mag. * (cos (VMEPadj .phase*pi/180)  +  ... 
i*sin (VMEPadj .phase*pi/180)  )  ; 
for  j=l : 20 

diff ( j ,  : ) =abs (VMEPcplx ( j  ,  :  )  - 

(A\vibset ( : , j ) )  .  ' ) ; % (Coef f { j } \vibset ( : , j ) )  . ' )  ;  % 

end 

max (diff) 
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%  Capt  Miller 

%  This  function  calculates  AVA  adjustments  for  the  AH-64A. 
function  adj  =  LinearAH64A(  vibsCmplx  ) 

O, 

o 

%  This  function  uses  the  AH-64A  AVA  coefficients  to  calculate 


o, 

o 

adjustments  based  on 

the  values  of 

vibsCmplx 

o, 

o 

Matrix  of  coefficient 

magnitudes 

mag  =  [0.000330,  0.044, 

0.000,  0.000, 

0.000; 

%  FPG100 

(Lat) 

0.000370,  0.143, 

0.000,  0.000, 

0.000; 

%  Hover 

(Lat) 

0.000331,  0.054, 

0.141,  0.255, 

0.286; 

%  60 

kts 

(Vert) 

0.000269,  0.062, 

0.227,  0.329, 

0.363; 

%  80 

kts 

(Vert) 

0.000390,  0.106, 

0.264,  0.434, 

0.485; 

%  100 

kts 

(Vert) 

0.000369,  0.156, 

0.405,  0.613, 

0.630; 

%  120 

kts 

(Vert) 

0.000287,  0.224, 

0.436,  0.664, 

0.689]  ; 

%  140 

kts 

(Vert) 

o, 

o 

Matrix  of  coefficient 

phases  (deg) 

phase  =  [165,  21,  0, 

0,  0; 

%  FPG100 

(Lat) 

169,  51,  0, 

o 

o 

%  Hover 

(Lat) 

231,  274,  255, 

246,  273; 

%  60  kts 

(Vert) 

234,  283,  267, 

260,  270; 

%  80  kts 

(Vert) 

232,  265,  269, 

260,  273; 

%  100  kts 

(Vert) 

242,  240,  262, 

253,  266; 

%  120  kts 

(Vert) 

250,  239,  249, 

245,  267]; 

%  140  kts 

(Vert) 

o, 

o 

Convert  the  coefficients  to  complex 

numbers 

A 

=  mag.* (cos (  phase*pi/180  )  +  i*sin(  phase*pi/180  ) 

) ; 

%  Calculate  the  adjustments 
adjCmplx  =  A\vibsCmplx; 

%  Convert  the  adjustments  to  magnitude  and  phase 

adj .mag  =  abs (  adjCmplx  ) ; 

adj .phase  =  angle (  adjCmplx  ) *180/pi; 
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%  Capt  Miller 

%  This  function  creates  a  single  coefficient  matrix  from  N  coefficient 
%  matrices  by  summing  the  squares  of  the  N  coefficient  matrix 
%  components  and  then  taking  the  square  root. 

function  [LSC,  standard] =Least_Squares_Coeff (Coef f) 

numcoef f=length (Coef f )  ; 
numrows=length (Coef f  { 1 }  (  :  ,  1 )  )  ; 
numcollumns=length (Coef f  { 1 }  ( 1 ,  :)); 

C=zeros (numrows , numcollumns )  ; 

for  j=l: numcoef f 

A.mag=abs (Coeff {  j  } )  ; 

A.phase=angle (Coeff  {  j } )  ; 

C=C+ (A. mag . A2 ) . * (cos (  A. phase  )  +  i*sin(  A. phase  )  ); 

end 

LSC.mag=sqrt (abs (C) /numcoeff)  ; 

LSC . phase=angle (C) *180/pi; 
for  j  =1 : numrows 

for  k=l : numcollumns 

if  LSC . phase ( j , k)  <  0 

LSC . phase ( j , k)  =  LSC .phase (j , k) +360; 

end 

end 

end 

LSCcmplx  =  LSC. mag. * (cos ( LSC. phase *pi/ 180) +i*sin (LSC. phase *pi/l 80) ) ; 
for  j=l:numcoeff 

dif f { j } =abs (Coeff { j } -LSCcmplx) ; 


end 

for  j  =1 : numrows 

for  k=l : numcollumns 
for  1=1: numcoeff 

temp ( 1 ) =dif f { 1 }  ( j  ,  k) ; 

end 

standard ( j , k)  =  std(temp); 

end 

end 


ooooooooooooooooooooooooooooooooooooooooooooooooooooooooo 
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Capt  Miller 


%  This  program  will  maesure  the  difference  between  the  VMEP  adjustment 

%  and  the  ad  hoc  adjustment  for  each  flight.  This  is  done  for  normal 

%  ad  hoc  coeffcients  and  ad  hoc  minus  zero  coefficients.  Intent  is  to 

%  determine  if  the  minus  zero  produces  any  differences  in  adjustments. 

clear ; clc; close  all; 

load  Coeff.mat 
load  AH64A  Vibes. mat 
load  AH64A  Coeff  Data. mat 
load  VMEPadj .mat 

[LSC,  standard] =Least_Squares_Coef f (Coeff ) ; 
n=2  0; 

for  j =1 : length (Coeff ) 

mag  =  AH64A  Vibes ( j ,  [1  5  9  13  17  21  25] )  '  ; 
phase  =  AH64A  Vibes ( j ,  [2  6  10  14  18  22  26] )  '  ; 

vibsCmplx  =  mag.* (cos (  phase*pi/180  )  +  i*sin(  phase*pi/180  )  ); 

pmag  =  AH64A  Vibes (j , [3  7  11  15  19  23  27])'; 
pphase  =  AH64A  Vibes ( j ,  [4  8  12  16  20  24  28] )  '  ; 

pvibsCmplx  =  pmag.* (cos (  pphase*pi/180  )  +  i*sin(  pphase*pi/180  ) 

)  ; 

vibset ( : , j )  =  -vibsCmplx  +  pvibsCmplx; 

end 

for  j=l : 10 

tmp { 2 * j -1 ,  : } =Coef f { 2 * j } ; %  tmp(odd)  =  A (even) 
tmp { 2* j , : } =Coef f { 2* j -1 } ; %  tmp (even)  =  A (odd) 

end 

%clear  Coeff; 

%Coef f=tmp; 

VMEPcplx  =  VMEPadj .mag. * (cos (VMEPadj .phase*pi/180)  +  ... 
i*sin (VMEPadj .phase*pi/180) ) 

A  =  LSC .mag .* (cos (  LSC.phase*pi/180  )  +  i*sin(  LSC.phase*pi/180  )  ); 

for  j=l : 20 

%adj ( j , : ) =A\vibset ( : , j ) ;  %  This  one  for  plotting  RMS  ad  hoc 
adj ( j , : ) =Coef f { j } \vibset ( : , j ) ; %This  one  for  plotting  individual  ad 
hoes 
end 

for  type=l:5 

%  figure (type) 
subplot (2,3, type) 

polar (angle (adj  ( 1 : n, type) ) , abs (adj  ( 1 : n, type) ) ,  '  ko  '  ) 
hold  on 

polar (angle (VMEPcplx ( 1 ; n, type) ) , abs (VMEPcplx ( 1 : n, type) ) , ' k+ ' ) 
end 

for  j=l : 20 

diff  ( j ,  : ) =abs (VMEPcplx ( j ,  : )  -  adj  ( j ,  : ) ) ; 

end 

diff 

max (diff) 
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%f igure  ( 1 ) 
subplot  (2,3, 1) 
title ( ' Weight ' ) 

%legend('Ad  Hoc  Ad j ustment Actual  Adjustment') 
%f igure (2 ) 
subplot (2,3,2) 
title (' Pitch  Link') 

%legend('Ad  Hoc  Ad j ustment ',' Actual  Adjustment') 

%f igure ( 3 ) 

subplot (2,3,3) 

title ('Trim  Tab  8-10') 

%legend('Ad  Hoc  Ad j ustment ',' Actual  Adjustment') 

subplot (2,3,4) 

title ('Trim  Tab  6-10') 

subplot (2,3,5) 

title ('Trim  Tab  4-10') 
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%  Capt  Nathan  A  Miller 

%  This  program  determines  the  difference  between  the  AVA  algorithm  and 
the 

%  PC-GBS  and  creates  a  histogram  of  the  results  for  each  adjustment 
type. 

clear;  clc; 
load  AVAadj .mat 
load  VMEPadj .mat 

AVAcomplex  =  AVAadj .mag. * (cos (  AVAadj .phase*pi/180  )  +... 

i*sin(  AVAadj .phase*pi/180  )  ); 

VMEPcomplex  =  VMEPadj .mag. * (cos (  VMEPadj .phase*pi/180  )  +... 

i*sin(  VMEPadj .phase*pi/180  )  ); 

Magnitude  diff  =  abs (AVAcomplex  -  VMEPcomplex) ; 

Rel  mag  diff  =  Magnitude  diff . /max (VMEPadj .mag,  AVAadj .mag) ; 

for  i=l : 20 

if  abs (AVAcomplex ( i , 1 ) )  <  113/2 

if  abs (VMEPcomplex ( i , 1 ) )  <  113/2 
Rel  mag  diff(i,l)=0; 

abs (AVAcomplex ( i , 1 ) ) ; 
abs (VMEPcomplex ( i ,  1 )  )  ; 
end 

end 

if  abs (AVAcomplex ( i , 2 ) )  <  0.5/2 

if  abs (VMEPcomplex ( i , 2 ) )  <  0.5/2 
Rel  mag  diff(i,2)=0; 

abs (AVAcomplex ( i , 2 ) ) ; 
abs (VMEPcomplex ( i ,  2  )  )  ; 
end 

end 
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if  abs (AVAcomplex ( i , 3 ) )  <  0.5/2 

if  abs (VMEPcomplex ( i , 3 ) )  <  0.5/2 
Rel_mag_dif f ( i , 3 ) =0 ; 
end 

end 

if  abs (AVAcomplex ( i , 4 ) )  <  0.5/2 

if  abs (VMEPcomplex ( i , 4 ) )  <  0.5/2 
Rel  mag  diff(i,4)=0; 
end 

end 

if  abs (AVAcomplex ( i , 5 ) )  <  0.5/2 

if  abs (VMEPcomplex ( i , 5 ) )  <  0.5/2 

Rel_mag_dif f ( i , 5 ) =0 ; 

end 

end 

end 

X  =  ( [012345678  9]+0.5)/10; 
for  i=l:5 

[n  ( i ,  : ) , xout ( i ,  : ) ] =hist (Rel_mag_dif f (  :  ,  i )  ,  X) ; 
subplot ( 1 , 5 , i ) ; 

bar (xout (i,  : ) *  100 , n  ( i ,  : ) / 0 . 2 ,  ' k ' ) ; 

AXIS ( [0  100  0  100] ) 

end 

subplot (1,5,1) 

title ( 'Weight ' , ' Font Size ',16) 

ylabel (' Percentage  of  Total  Flights  Sampled' ,' FontSize ', 16) 
subplot (1,5,2) 

title (' Pitch  Link', ' FontSize ', 16) 
subplot (1,5,3) 

title ('Tab  8-10 ',' FontSize ', 16) 

xlabel (' Percent  Difference  of  PC-GBS  Adjustments  from  AVA  Algorithm 
Adjustments ' , ' FontSize ' , 16) 
subplot  (1,5,4) 

title ('Tab  6-10 ',' FontSize ', 16) 
subplot (1,5,5) 

title ('Tab  4-10 ',' FontSize ', 16) 
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%  Capt  Miller  Thesis 

%  This  is  a  test  program  to  determine  if  any  mistakes  were  made 
%  during  data  entry, 
clear  all;  close  all;  clc; 

load  Coeff.mat 

[LSC,  standard] =Least_Squares_Coef f (Coeff) ; 

[LSCgood,  standardgood] =Least_Squares_Coef f (Coef f (1:5) ) ; 

[LSCexd,  standardexd] =Least_Squares_Coef f (Coef f (15:20) ) ; 
b=l  ; 

e=length (Coef f )  ; 
n=l  ; 

for  h=l : 7 
for  j=b:e 

value (h, j  +  l-b,  : ) =Coef f { j  }  (h,  : )  ; 
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end 


end 

vertvalue  =  value  (3 : 7, : ,  : )  ; 
latvalue  =  value  ( 1 : 2 ,  : ,  : )  ; 
for  h=l : 7 

maxi (h) =max (abs (value (h,  : , 1)  )  )  ; 
max2 (h) =max (abs (value  (h,  : ,  2 )  )  )  ; 
max3 (h) =max (abs (value (h,  : , 3)  )  )  ; 
max4 (h) =max (abs (value  (h,  : , 4 ) ) )  ; 
max5 (h) =max (abs (value (h,  : , 5)  )  )  ; 

end 

for  h=l : 5 

vertmaxl (h) =max (abs (vertvalue (h, : ,  1)  )  )  ; 
vertmax2 (h) =max (abs (vertvalue (h,  : ,  2 )  )  )  ; 
vertmax3 (h) =max (abs (vertvalue (h,  : ,3) ) ) ; 
vertmax4 (h) =max (abs (vertvalue (h, : , 4 ) ) ) ; 
vertmax5 (h) =max (abs (vertvalue (h,  : ,  5)  )  )  ; 

end 

for  h=l:2 

latmaxl (h) =max (abs (latvalue (h,  : , 1)  )  )  ; 
latmax2 (h) =max (abs (latvalue (h,  : , 2)  )  )  ; 
latmax3 (h) =max (abs (latvalue (h, : , 3)  )  )  ; 
latmax4 (h) =max (abs (latvalue (h,  : , 4) ) )  ; 
latmax5 (h) =max (abs (latvalue (h, :  , 5)  )  )  ; 

end 

for  h=l : 7 
subplot ( 3 , 3 ,  h) 

%polar (angle (value (h,  : ,  1 )  )  ,  abs (value (h,  : , 1 ) ) /maxi (h) ,  ' bx '  ) 

polar (angle (value (h,  : , n) ) , abs (value (h,  : , n) ) ,  ' rx  '  ) 

end 

subplot  (3,3,1) 
title (' FPG100  Lat') 
subplot (3,3,2) 
title ('Hover  Lat') 
subplot (3,3,3) 
title (' 60  Kts  Vert') 
subplot  (3,3,4) 
title (' 80  Kts  Vert') 
subplot (3,3,5) 
title ('100  Kts  Vert') 
subplot (3,3,6) 
title ('120  Kts  Vert') 
subplot (3,3,7) 
title  ('140  Kts  Vert') 

figure  (2 ) 
for  i=l:5 

tempi  =  max (vertmaxl ( i ), vertmax2 ( i )) ; 
temp2  =  max (vertmax3 ( i ) , vertmax4 ( i ) ) ; 
temp3  =  max (tempi ,temp2); 


F-14 


maxval  =  max (temp3, vertmax5 (i) ) ; 

subplot (2, 3, i)  %  Each  subplot  is  for  a  different  lateral  regime 
polar ( 0 , maxval , '  w . ' ) 
hold  on 

polar (angle (vertvalue ( i ,  : , 1 ) ) , abs (vertvalue ( i ,  : ,  1 ) )  ,  ' kx '  ) 
polar (angle (vertvalue (i, : , 2) ) , abs (vertvalue ( i , : , 2 ) ) , ' kd ' ) 
polar (angle (vertvalue (i,  : , 3) ) , abs (vertvalue ( i ,  : , 3 ) )  ,  ' kA ' ) 
polar (angle (vertvalue (i, : , 4) ) , abs (vertvalue ( i , : , 4 ) ) , ' kv ' ) 
polar (angle (vertvalue (i, : , 5) ) , abs (vertvalue ( i , : , 5 ) ) , ' kp ' ) 
end 

figure ( 3 ) 
for  i=l:2 

tempi  =  max (latmaxl (i) , latmax2 (i) ) ; 
temp2  =  max (latmax3 (i) , latmax4 (i) ) ; 
temp3  =  max (tempi , temp2 ) ; 
maxval  =  max (temp3, latmax5 (i) ) ; 

subplot ( 1 , 2 , i )  %  Each  subplot  is  for  a  different  lateral  regime 

polar ( 0 , maxval , '  w . ' ) 
hold  on 

polar (angle ( latvalue ( i ,  : , 1 ) ) , abs ( latvalue (i,  :  ,  1 ) *  100 )  ,  'kx') 
polar (angle (latvalue (i,  : , 2) ) , abs (latvalue (i,  : , 2) ) ,  ' kd' ) 
polar (angle (latvalue (i,  : , 3) ) , abs (latvalue (i,  : , 3) ) ,  ' kA ' ) 
polar (angle (latvalue (i, : , 4) ) , abs (latvalue (i, : , 4) ) , ' kv' ) 
polar (angle (latvalue (i,  : , 5) ) , abs (latvalue (i,  : , 5) ) ,  ' kp ' ) 
end 

figure  ( 4 ) 

polar ( 0 , maxval , ' w . ' ) 
hold  on 

polar (angle (latvalue (i, : , 1) ) , abs (latvalue (i, : , 1) *100) , ' ko ' ) 
polar (angle (latvalue (i,  : , 2) ) , abs (latvalue (i,  : , 2) )  ,  ' kd' ) 
polar (angle (latvalue (i,  : , 3) ) , abs (latvalue (i,  : , 3) ) ,  ' kA ' ) 
polar (angle (latvalue (i, : , 4) ) , abs (latvalue (i, : , 4) ) , ' kv' ) 
polar (angle (latvalue (i,  : , 5) ) , abs (latvalue (i,  : , 5) ) ,  ' kp ' ) 
title ('Ad  Hoc  Coefficients  for  Hover (Lat)  for  AH-64A') 

abs (LSCexd.mag-LSCgood.mag) . /LSCgood.mag*100 
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Appendix  G:  AH-64D  Matlab  Code 

Note:  The  Matlab  codes  contained  in  this  appendix  were  each  used  to  perform  one  or 
more  specific  functions  in  the  course  of  this  research.  There  are  several  instances  where 
parts  of  the  code  have  been  commented  out  with  the  %  symbol.  In  order  to  reproduce  all 
of  the  analysis  of  this  thesis,  some  lines  of  code  may  need  to  be  un-commented  and  other 
lines  commented  out. 

%  Capt  Nathan  A  Miller 

%  This  program  creates  MAT  data  files  from  imported  Excel  database 
%  Ensure  that  the  two  databases  are  imported  into  matlab  prior  to 
running 

function  Excel  to  MAT 
global  data  data2 
n  =  length (data2 ( : , 2 ) ) -2 

%  Read  Vibration  Magnitude  and  Phase  as  well  as  prediction  magnitude 
%  and  phase  from  Excel  File  and  assign  to  data  matrix 
Longbow  Vibes  =  zeros (n, 28); 
for  j  =1 : n 

for  k=l : 7 

Longbow_Vibes ( j , 1  +  4*  (k-1) )  =  data ( l  +  k+ ( j -1 ) *11 , 2 ) ; 

Longbow_Vibes ( j , 2+4* ( k— 1 ) )  =  data ( l+k+ ( j -1 ) *11 , 3) ; 

Longbow_Vibes ( j , 3+4* (k-1) )  =  data ( l+k+ ( j -1 ) *11 , 4 ) ; 

Longbow_Vibes ( j , 4  +  4* (k-1) )  =  data ( l  +  k+ ( j -1 ) *11 , 5)  ; 

end 

end 

save  Longbow  Vibes. mat  Longbow  Vibes 
Longbow  VMEP  Adjustments  =  zeros (n, 20); 
for  j  =1 : n 

for  k=l : 20 

Longbow  VMEP  Adj ustments ( j , k)  =  data2 ( j +2 , k+1 ) ; 
end 

end 

save  Longbow  VMEP  Adj ustments . mat  Longbow  VMEP  Adjustments 
Longbow  Coeff  Data  =  zeros (n, 75); 
for  j  =1 : n 

Longbow_Coef f_Data ( j  ,  1 ) =data (2+ (j— 1)  *11,  6)  ; 

Longbow_Coef f_Data ( j , 16) =data (2+ ( j-1) *11, 9) ;  %  Records  manual  adjustment 
Longbow_Coef f_Data ( j , 31 ) =data (2+ ( j -1 ) *11 , 12 ) ;  %  that  produces  delta  vibes 

Longbow_Coef f_Data ( j , 4  6) =data (2+ (j— 1) *11, 15) ; 

Longbow_Coef f_Data ( j , 61 ) =data (2  +  (j— 1) *11, 18) ; 
for  k=l : 7 
for  z=l:5 

Longbow_Coeff _Data ( j, 2*k+15* (z-1) )  =  data (l+k+ (j-1) *11, 7  +  3*  (z-1) ) ; 
Longbow_Coeff_Data ( j, 2*k+l  +  15* (z-1) ) =data (l+k+ (j-1) *11,  8  +  3* (z-1) ) ; 
end 

Longbow  Coeff  Data ( j , 74+2 *k) =data ( l  +  k+ ( j -1 ) *  1 1 , 2 1 )  ; 

Longbow  Coeff  Data ( j , 75+2 *k) =data ( l  +  k+ ( j -1  )* 1 1 , 22 )  ; 

end 

end 

save  Longbow_Coef f  Data. mat  Longbow  Coeff  Data 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
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%  Capt  Miller 

%  This  program  converts  the  coefficient  matrix  data  into  actual 
%  coefficient  matrices.  The  program  loads  the  raw  vibe  data  from 
%  AH64A  Vibes. mat  and  computes  a  delta  vibe  from  the  predicted  vibe 
%  values  stored  in  AH64A  Coeff  Data. mat.  The  coefficient  matrix  is 
%  computed  column  by  column.  Once  computed,  the  coefficient  matrix 
%  operates  on  a  vibe  set  to  produce  a  quasi-linear  adjustment  set. 

%  This  process  is  repeated  for  each  coefficient  data  set.  The  same 
%  vibe  set  is  always  used  for  comparison. 

function  VMEP  Coefficient  Matrix 

FLG=-1; 

delta=0 ; 

load  Longbow  Vibes. mat 
load  Longbow_Coef f  Data. mat 

n=length (Longbow  Vibes (:, 1) ); 
for  j  =1 : n 

wt  =  Longbow_Coef f_Data ( j , 1 ) ; 
pi  =  Longbow_Coeff_Data ( j  ,  16)  ; 
tab8  =  Longbow_Coeff_Data ( j , 31) ; 
tab6  =  Longbow_Coeff_Data ( j  ,  46)  ; 

tab4  =  Longbow_Coef f_Data ( j  ,  61 )  ; 

for  k=l  :  7  %  loops  over  the  7  flight  conditons,  FPG100  through  140  Knots 
vibmag=Longbow  Vibes ( j  ,  1  +  4* (k-1)  )  ; 
vibphase=Longbow  Vibes ( j , 2  +  4* (k-1)  )  ; 
vibcmplx=vibmag* (cos (  vibphase*pi/180  )  +  ... 
i*sin(  vibphase*pi/180  )  ); 


pvibwtmag=Longbow  Coeff  Data (j , 2+2* (k-1) ) ; 
pvibwtphase=Longbow  Coeff  Data (j , 3+2* (k-1) ) ; 
pvibwtcmplx=pvibwtmag* (cos (  pvibwtphase*pi/180  )  +  ... 
i*sin(  pvibwtphase*pi/180  )  ); 

pvibplmag=Longbow  Coeff  Data (j , 17+2* (k-1) ) ; 
pvibplphase=Longbow  Coeff  Data (j , 18+2* (k-1) ) ; 
pvibplcmplx=pvibplmag* (cos (  pvibplphase*pi/180  )  +  ... 
i*sin(  pvibplphase*pi/180  )  ); 

pvibtab8mag=Longbow  Coeff  Data (j , 32+2* (k-1) ) ; 
pvibtab8phase=Longbow  Coeff  Data (j , 33+2* (k-1) ) ; 
pvibtab8cmplx=pvibtab8mag* (cos (  pvibtab8phase*pi/180  )  +  ... 
i*sin(  pvibtab8phase*pi/180  )  ); 

pvibtab6mag=Longbow  Coeff  Data (j , 47+2* (k-1) ) ; 
pvibtab6phase=Longbow  Coeff  Data (j , 48+2* (k-1) ) ; 
pvibtab6cmplx=pvibtab6mag* (cos (  pvibtab6phase*pi/180  )  +  ... 
i*sin(  pvibtab6phase*pi/180  )  ); 

pvibtab4mag=Longbow  Coeff  Data (j , 62+2* (k-1) ) ; 
pvibtab4phase=Longbow  Coeff  Data (j , 63+2* (k-1) ) ; 
pvibtab4cmplx=pvibtab4mag* (cos (  pvibtab4phase*pi/180  )  +  ... 
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i*sin(  pvibtab4phase*pi/180  )  ); 

if  FLG  >  0 

zeromag=Longbow  Coeff  Data (j , 76+2* ( k— 1 ) ) ; 
zerophase=Longbow  Coeff  Data (j , 77+2* ( k— 1 ) )  ; 
zerocmplx=zeromag* (cos (  zerophase*pi/180  )  +  ... 

i*sin(  zerophase*pi/180  )  ); 

delta=vibcmplx-zerocmplx; 

end 

dvibwt=pvibwtcmplx-vibcmplx-delta; 

dvibpl=pvibplcmplx-vibcmplx-delta; 

dvibtab8=pvibtab8cmplx-vibcmplx-delta; 

dvibtab6=pvibtab6cmplx-vibcmplx-delta; 

dvibtab4=pvibtab4cmplx-vibcmplx-delta; 

A ( k, 1 ) =dvibwt/wt ; 

A ( k, 2 ) =dvibpl/pl ; 

A ( k, 3 ) =dvibtab8 / tab8 ; 

A ( k, 4 ) =dvibtab6/tab6; 

A ( k, 5 ) =dvibtab4 /tab4 ; 


end 

Coeff { j } =A; 


end 

save  Coeff. mat  Coeff 


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 


%  Capt  Nathan  A  Miller 

%  Converts  VMEP  adjustments  into  magnitude  and  phase  of  adjustments, 
function  VMEP  2  MagPhase 

load  Longbow  VMEP  Adjustments .mat 


n  =  length (Longbow  VMEP  Adjustments (:, 1) ) ; 

a=Longbow  VMEP  Adjustments; 
clear  Longbow  VMEP  Adjustments 

for  j  =1 : n 


O, 

o 

adj  cmplx ( j , 1 ) 
adj  cmplx ( j , 2 ) 
adj  cmplx ( j , 3 ) 
adj  cmplx ( j , 4 ) 
adj  cmplx ( j , 5 ) 


12  3  4 

+a  (j , 1) -i*a(j, 6) -a  (j , 11) +i*a  ( j , 16) ; 
+a ( j , 2) -i*a (j , 7) -a (j  ,  12) +i*a (j  ,  17)  ; 
+a (j , 3) -i*a (j , 8) -a (j , 13) +i*a  (j , 18) ; 
+a  (j  ,  4)  -  i*a(j,  9)  -a  (j  ,  14)  +i*a  (j  ,  19)  ; 
+a ( j , 5) -i*a ( j , 10) -a ( j , 15) +i*a(j,20) 


%  Weight  Adj 
%  P/L  Adj 
%  Tab  Adj  8 
%  Tab  Adj  6 
%  Tab  Adj  4 


end 

VMEPadj .mag=abs (adj cmplx) ; 

VMEPadj . phase  =  angle (adjcmplx) *180/pi; 

save  VMEPadj. mat  VMEPadj 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
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%  Capt  Nathan  A  Miller 
%  Thesis  work 

%  This  program  loads  vibe  data  and  calculates  an  AVA  adjustment  matrix 

function  Vibs2AVA  ADJ 

load  Longbow  Vibes. mat 
n  =  length (Longbow  Vibes (:, 1) ); 

for  j  =1 : n 

mag  =  Longbow  Vibes (j,[l  5  9  13  17  21  25])'; 
phase  =  Longbow  Vibes (j, [2  6  10  14  18  22  26])'; 

vibsCmplx  =  mag.* (cos (  phase*pi/180  )  +  i*sin(  phase*pi/180  )  ); 

pmag  =  Longbow  Vibes (j, [3  7  11  15  19  23  27])'; 
pphase  =  Longbow  Vibes (j, [4  8  12  16  20  24  28])'; 

pvibsCmplx  =  pmag.* (cos (  pphase*pi/180  )  +  i*sin(  pphase*pi/180  )  ); 

delvib  =  -vibsCmplx  +  pvibsCmplx; 

adj(j)  =  LinearLongbow (  delvib  ); 

end 

for  k=l : n 

magnitude ( k, : )  =  adj(k).mag'; 

Phase (k, :)  =  ad j ( k) . phase ' ; 

end 

AVAadj.mag  =  magnitude; 

AVAadj. phase  =  Phase; 
save  AVAadj .mat  AVAadj 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

%  Capt  Miller 

%  This  file  computes  ad  hoc  coefficients  using  all  flights,  only  good 
%  flights,  and  only  exceedence  flights.  These  coefficients  are  then 
%  compared  by  producing  adjustment  sets  for  the  20  sampled  flights. 

%  The  differences  in  adjustments  are  then  compared  based  on  the  "Good" 
%  flight  ad  hoc  matrix  and  the  "Exceed"  flight  ad  hoc  matrix.  If  the 
two 

%  matrices  produce  adjustments  that  are  within  the  tolerance  of 
mechanical 

%  limitations  of  the  aircraft,  then  the  adjustment  types  are 
identical . 

clear,  close  all,  clc; 

load  Coeff.mat 

load  Longbow  Vibes. mat 

[LSC,  standard] =Least_Squares_Coeff (Coef f) ; 

[LSCgood,  standardgood] =Least_Squares_Coef f (Coef f (1:5) ) ; 

[LSCexd,  standardexd] =Least_Squares_Coef f (Coef f (15:20) ) ; 

LSCgoodcplx  =  LSCgood. mag. * (cos (LSCgood. phase*pi/180)  +  ... 
i*sin (LSCgood. phase*pi/180) ) ; 
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LSCexdcplx  =  LSCexd.mag. * (cos (LSCexd.phase*pi/180)  +  ... 
i*sin (LSCexd.phase*pi/180) ) ; 

for  j =1 : length (Coeff) 

mag  =  Longbow  Vibes ( j , [ 1  5  9  13  17  21  25])'; 
phase  =  Longbow  Vibes (j, [2  6  10  14  18  22  26])'; 

vibsCmplx  =  mag.* (cos (  phase*pi/180  )  +  i*sin(  phase*pi/180  )  ); 

pmag  =  Longbow  Vibes (j , [3  7  11  15  19  23  27])'; 
pphase  =  Longbow  Vibes (j, [4  8  12  16  20  24  28])'; 

pvibsCmplx  =  pmag.* (cos (  pphase*pi/180  )  +  i*sin(  pphase*pi/180  ) 

)  ; 

vibset ( : , j )  =  -vibsCmplx+pvibsCmplx; 

end 

for  j =1 : length (Coeff ) 

discrete_vib=vibset  (  : , j ) ; 

adjgood ( j , : ) =LSCgoodcplx\discrete_vib; 

adj  exceed ( j ,  : ) =LSCexdcplx\discrete_vib; 

end 

diff  =  abs (adjgood-adjexceed) 
maxdiff(l:5)  =  max (dif f ( : ,  1 : 5 )  ) 

PctDiffCoeff  =  abs (LSCgoodcplx-LSCexdcplx) . /abs (LSCgoodcplx) *100 


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

%  Capt  Miller 

%  This  program  plots  AVA,  VMEP,  and  ad  hoc  adjustments  together  for 
%  comparison.  It  also  plots  RMS  of  ad  hoc. 

clear; clc; elf ; 

load  Coeff. mat 

load  Longbow  Vibes. mat 

load  VMEPadj.mat 

%  Runs  the  Least  Squares  Coefficient  function  to 

%  generate  the  LSC  matrix  from  the  Coeff  cell. 
[LSC,  standard] =Least_Squares_Coeff (Coeff) 
vib  num=16 
b=l 

e=length (Coeff) 

for  j =1 : length (Coeff ) 

mag  =  Longbow  Vibes (j,[l  5  9  13  17  21  25])'; 
phase  =  Longbow  Vibes (j, [2  6  10  14  18  22  26])'; 

vibsCmplx  =  mag.* (cos (  phase*pi/180  )  +  i*sin(  phase*pi/180  )  ); 

pmag  =  Longbow  Vibes (j, [3  7  11  15  19  23  27])'; 
pphase  =  Longbow  Vibes (j, [4  8  12  16  20  24  28])'; 

pvibsCmplx  =  pmag.* (cos (  pphase*pi/180  )  +  i*sin(  pphase*pi/180  ) 

)  ; 

vibset (:, j )  =  -vibsCmplx+pvibsCmplx; 

end 

discrete  vib=vibset ( : , vib  num) ;  %  Select  discrete  vibe  set  to  use  for 
comparison 
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for  j =1 : length (Coeff) 

adj ( : , j ) =Coef f { j } \discrete_vib; 

Cmag ( : , j )  =  abs ( adj ( : ,  j  )  )  ; 

Cphase ( : , j )  =  (angle (adj ( : , j ) ) ) *180/pi; 

end 

VMEP2AVA.mag  =  Cmag'; 

VMEP2AVA. phase  =  Cphase'; 
save  VMEP2AVA.mat  VMEP2AVA 

%  Calculate  the  Least  Squares  Method  adjustment 
A  =  LSC.mag. * (cos (  LSC.phase*pi/180  )  +  i*sin(  LSC.phase*pi/180  )  ); 

LS_adj  =A\discrete_vib; 

AVA  adj =LinearLongbow (  discrete  vib  ); 

VMEPmag=VMEPadj .mag (vib  num,  : ) 

VMEPphase=VMEPadj .phase (vib  num,  : ) 

%  Determine  the  max  magnitude  so  we  can  normalize  our  adjustments  for 
%  use  in  plotting  on  a  single  graph 
for  k=l : 5 

Maxl=max (max (Cmag  ( k, b :  e) )  ,  AVA  adj. mag (k) ) ; 

%Maxl=max (abs (AVA  adj .mag (k) ) ) ; 

Max2=max (VMEPmag (k) ,  abs (LS  adj (k) ) ) ; 

Max ( k) =max (Maxi , Max2 ) ; 

%Max ( k ) =max ( Cmag ( k , b : e ) ) ; 

end 

figure ( 1 ) 
subplot (2,3,1) 
polar (0,1, ' . w ' ) 
hold 

%  Plot  weight  adjustments 

%  polar (Cphase ( 1 , b : e) *pi/l 80 , Cmag ( 1 , b : e) /Max ( 1 ) , ' xr ' ) 

polar (angle (LS_adj  ( 1 ) ) , abs (LS_adj  ( 1 ) ) /Max ( 1 ) ,  ' sr '  ) 
polar (AVA  adj .phase (1) *pi/180, AVA  adj .mag ( 1 ) /Max ( 1 ) , ' or ' ) 
polar (VMEPphase (1) *pi/180, VMEPmag ( 1 ) /Max ( 1 ) , ' *r ' ) 

%  Plot  PCR  adjustments 

%  polar (Cphase (2 , b : e) *pi/l 80 , Cmag (2 , b : e) /Max (2 ) , ' xk ' ) 

polar (angle (LS_adj  (2 ) ) , abs (LS_adj  (2 ) ) /Max (2 )  ,  ' sk '  ) 
polar (AVA  adj .phase (2) *pi/180, AVA  adj .mag (2) /Max (2) ,' ok' ) 
polar (VMEPphase (2) *pi/l 80 , VMEPmag (2 ) /Max (2 ) , ' *k ' ) 

%  Plot  tab  8  adjustments 

%  polar (Cphase ( 3 , b : e) *pi/l 80 , Cmag ( 3 , b : e) /Max ( 3 ) , ' xb ' ) 

polar (angle (LS_adj  (3) )  ,  abs (LS_adj  (3) ) /Max (3) ,  ' sb ' ) 
polar (AVA_adj . phase (3) *pi/l 80 , AVA_adj . mag ( 3 ) /Max ( 3 ) , ' ob ' ) 
polar (VMEPphase (3) *pi/l 80 , VMEPmag ( 3 ) /Max ( 3 ) , ' *b ' ) 

%  Plot  tab  6  adjustments 

%  polar (Cphase ( 4 , b : e) *pi/l 80 , Cmag ( 4 , b : e) /Max ( 4 )  ,  ' xg ' ) 
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polar (angle (LS_adj  ( 4 ) ) , abs (LS_adj  ( 4 ) ) /Max ( 4 ) ,  ' sg '  ) 
polar (AVA  adj .phase (4) *pi/180, AVA  adj . mag ( 4 ) /Max ( 4 ) , ' og ' ) 
polar (VMEPphase (4) *pi/l 80 , VMEPmag ( 4 ) /Max ( 4 ) , ' *g ' ) 

%  Plot  tab  4  adjustments 

%  polar (Cphase ( 5 , b : e) *pi/l 80 , Cmag ( 5 , b : e) /Max ( 5 ) , ' xm ' ) 

polar (angle (LS_adj  (5) ) , abs (LS_adj  (5) ) /Max (5) ,  ' sm'  ) 
polar (AVA_adj .phase (5) *pi/180, AVA_adj .mag (5) /Max (5) , ' om' ) 
polar (VMEPphase (5) *pi/l 80 , VMEPmag ( 5 ) /Max ( 5 ) , ' *m ' ) 

title (' Red-Wt,  Black-Pcr,  Blue-Tab8,  Green-Tab6,  Magenta-Tab4 ' ) 
subplot (2,3,2) 

polar (Cphase ( 1 , b : e) *pi/180, Cmag ( 1 ,  b : e)  ,  ' xr  '  ) 
subplot (2,3,3) 

polar (Cphase (2 , b : e) *pi/180, Cmag (2 , b : e) , ' xk ' ) 
subplot  (2,3,4) 

polar (Cphase ( 3 , b : e) *pi/180, Cmag ( 3 , b : e) , ' xb ' ) 
subplot (2,3,5) 

polar (Cphase ( 4 , b : e) *pi/l 80 , Cmag ( 4 , b : e) , ' xg ' ) 
subplot (2,3,6) 

polar (Cphase ( 5 , b : e) *pi/180, Cmag ( 5 , b : e) , ' xm ' ) 

VMEPcplx  =  VMEPadj .mag. * (cos (VMEPadj .phase*pi/180)  +  ... 

i*sin (VMEPadj .phase*pi/180) ) ; 
for  j=l : 20 

diff ( j , : ) =abs (VMEPcplx ( j , : )  -  (A\vibset ( : , j ) ) . ' ) 


end 

max (diff) 
figure  (2 ) 

polar (0,1, ' . w ' ) 
hold 

%  Plot  weight  adjustments 

polar (AVA  adj .phase (1) *pi/180, AVA  adj .mag ( 1 ) /Max ( 1 ) , ' or ' ) 
polar (VMEPphase (1) *pi/180, VMEPmag ( 1 ) /Max ( 1 )  ,  ' *r ' ) 

%  Plot  PCR  adjustments 

polar (AVA  adj .phase (2) *pi/180, AVA  adj .mag (2 ) /Max (2 ) , ' ok ' ) 
polar (VMEPphase (2) *pi/180, VMEPmag (2 ) /Max (2 ) , ' *k ' ) 

%  Plot  tab  8  adjustments 

polar (AVA_adj . phase (3) *pi/l 80 , AVA_adj .mag ( 3 ) /Max ( 3 ) , ' ob ' ) 
polar (VMEPphase (3) *pi/l 80 , VMEPmag ( 3 ) /Max ( 3 ) , ' *b ' ) 

%  Plot  tab  6  adjustments 

polar (AVA  adj .phase (4) *pi/180, AVA  adj .mag ( 4 ) /Max ( 4 ) , ' og ' ) 
polar (VMEPphase (4) *pi/l 80 , VMEPmag ( 4 ) /Max ( 4 ) , ' *g ' ) 

%  Plot  tab  4  adjustments 

polar (AVA_adj .phase (5) *pi/180, AVA_adj .mag (5) /Max (5) , ' om' ) 
polar (VMEPphase (5) *pi/180, VMEPmag ( 5 ) /Max ( 5 ) , ' *m ' ) 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
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function  adj  =  LinearLongbow (  vibsCmplx  ) 

%  This  function  calculates  adjustments  according  to  AVA  coefficients. 
%  Input: 

%  vibsCmplx  =  matrix  of  complex  vibration  data 
%  Output: 

%  adj  =  adjustment  vector  [hub  weight (oz),  pitch  link (notch),  3  X 
trim  tab (oz) ] 

O, 

o 


%  Matrix  of  AVA  coefficient  magnitudes 


0.0004936, 

0.04448, 

0, 

0, 

0; 

o, 

o 

FPG100 

(Lat) 

0.0004580, 

0.15600, 

0, 

0, 

0; 

o, 

o 

Hover 

(Lat) 

0.0004696, 

0.03924, 

0.1605, 

0.3385, 

0.6507; 

o, 

o 

60 

kts 

(Vert) 

0.0004546, 

0.06359, 

0.1714, 

0.2886, 

0.7391; 

o, 

o 

80 

kts 

(Vert) 

0.0004825, 

0.11370, 

0.1888, 

0.3168, 

0.6596; 

o, 

o 

100 

kts 

(Vert) 

0.0004547, 

0.18190, 

0.2161, 

0.3747, 

0.6992; 

o, 

o 

120 

kts 

(Vert) 

0.0004519, 

0.24340, 

0.3071, 

0.4467, 

0.9085] ; 

o, 

o 

140 

kts 

(Vert) 

%  Matrix  of  AVA  coefficient  phases  (deg) 


163.0, 

15.5, 

0, 

0, 

0; 

o, 

o 

FPG100 

(Lat) 

171.2, 

LO 

l> 

LO 

0, 

0, 

0; 

o, 

o 

Hover 

(Lat) 

211 . 9, 

286.6, 

263.5, 

270.6, 

256.1; 

o, 

o 

60  kts 

(Vert) 

204.7, 

273.2, 

261 . 6, 

261.3, 

259.2; 

g, 

o 

80  kts 

(Vert) 

215.9, 

262.0, 

258.8, 

268.4, 

256.7; 

o, 

o 

100  kts 

(Vert) 

219.6, 

256.1, 

255.8, 

258.6, 

255.0; 

g, 

o 

120  kts 

(Vert) 

235.6, 

247.2, 

250.3, 

260.5, 

250.5] ; 

g, 

o 

140  kts 

(Vert) 

%  Convert  the  coefficients  to  complex  numbers 
A  =  mag.* (cos (  phase*pi/180  )  +  i*sin(  phase*pi/180  )  ); 


%  Calculate  the  adjustments 
adjCmplx  =  A\vibsCmplx; 

%  Convert  the  adjustments  to  magnitude  and  phase 

adj .mag  =  abs (  adjCmplx  ) ; 

adj .phase  =  angle (  adjCmplx  ) *  180/pi; 


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 


%  Capt  Miller 

%  This  function  creates  a  single  coefficient  matrix  from  N  coefficient 
%  matrices  by  summing  the  squares  of  the  N  coefficient  matrix 
components 

%  and  then  taking  the  square  root. 

function  [LSC,  standard] =Least_Squares_Coef f (Coeff) 

numcoef f=length (Coef f )  ; 
numrows=length (Coeff  { 1 }  (  :  ,  1 )  )  ; 
numcollumns=length (Coeff { 1 }  ( 1 ,  :)); 

C=zeros (numrows , numcollumns )  ; 

for  j=l: numcoef f 

A.mag=abs (Coeff {  j  } )  ; 

A.phase=angle (Coeff  {  j  } )  ; 

C=C+ (A. mag. A2) . * (cos (  A. phase  )  +  i*sin(  A. phase  )  ); 

end 
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LSC.mag=sqrt (abs (C) /numcoeff)  ; 

LSC . phase=angle (C) *180/pi; 
for  j=l:numrows 

for  k=l : numcollumns 

if  LSC . phase ( j , k)  <  0 

LSC . phase ( j , k)  =  LSC. phase (j , k) +360; 

end 

end 

end 

LSCcmplx  =  LSC. mag. * (cos (LSC. phase *pi/ 180) +i*sin (LSC. phase *pi/ 180) ) ; 
for  j=l:numcoeff 

dif f { j } =abs (Coeff { j } -LSCcmplx) ; 

end 

for  j=l:numrows 

for  k=l : numcollumns 
for  1=1: numcoeff 

temp ( 1 ) =dif f { 1 }  ( j  ,  k)  ; 

end 

standard ( j , k)  =  std(temp); 

end 

end 


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

%  Capt  Miller 

%  This  program  will  maesure  the  difference  between  the  VMEP  adjustment 
and 

%  the  ad  hoc  adjustment  for  each  flight.  Plots  are  also  generated  of 
ad 

%  hoc  and  VMEP  adjustments. 


clear; clc; close  all; 

load  Coeff. mat 
load  Longbow  Vibes. mat 
load  Longbow_Coef f  Data. mat 
load  VMEPadj.mat 
n=2  0; 

[LSC,  standard] =Least_Squares_Coeff (Coeff) 
for  j =1 : length (Coeff ) 

mag  =  Longbow  Vibes (j , [1  5  9  13  17  21  25])'; 
phase  =  Longbow  Vibes (j, [2  6  10  14  18  22  26])'; 

vibsCmplx  =  mag.* (cos (  phase*pi/180  )  +  i*sin(  phase*pi/180  )  ); 

pmag  =  Longbow  Vibes (j, [3  7  11  15  19  23  27])'; 
pphase  =  Longbow  Vibes (j, [4  8  12  16  20  24  28])'; 

pvibsCmplx  =  pmag.* (cos (  pphase*pi/180  )  +  i*sin(  pphase*pi/180  ) 

)  ; 

vibset ( : , j )  =  -vibsCmplx+pvibsCmplx; 

end 
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VMEPcplx  =  VMEPadj .mag. * (cos (VMEPadj .phase*pi/180)  +  ... 
i*sin (VMEPadj .phase*pi/180) ) 

A  =  LSC.mag. * (cos (  LSC.phase*pi/180  )  +  i*sin(  LSC.phase*pi/180  )  ); 

for  j  =1 : 1 0 

tmp { 2* j -1 , : } =Coef f { 2* j } ; %  tmp(odd)  =  A (even) 
tmp { 2 * j ,  : } =Coef f { 2* j -1 } ; %  tmp (even)  =  A (odd) 

end 

clear  Coef f ; 

Coef f=tmp; 

for  j  =1 : n 

%adj ( j , : ) =A\vibset ( : , j ) ;  %  This  one  for  plotting  RMS  ad  hoc 

adj ( j , : ) =Coef f { j } \vibset ( : , j ) ; %This  one  for  plotting  individual  ad 

hoes 

end 

for  type=l:5 

subplot (2,3, type) 

polar (angle (adj ( 1 : n, type) ) , abs (adj ( 1 : n, type) ) , ' ko ' ) 
hold  on 

polar (angle (VMEPcplx ( 1 : n, type) ) , abs (VMEPcplx ( 1 : n, type) ) , ' k+ ' ) 
end 

for  j=l : 20 

diff  ( j ,  : ) =abs (VMEPcplx ( j ,  : )  -  adj  ( j ,  : ) ) ; 

end 

diff 

max (diff) 

%f igure ( 1 ) 
subplot (2,3,1) 
title ( ' Weight ' ) 

%legend('Ad  Hoc  Adj ustment Actual  Adjustment') 

%f igure (2 ) 
subplot  (2,3,2) 
title (' Pitch  Link') 

%legend('Ad  Hoc  Adj ustment ',' Actual  Adjustment') 

%f igure ( 3 ) 

subplot (2,3,3) 

title ( ' Trim  Tab  8-10') 

%legend('Ad  Hoc  Adj ustment ',' Actual  Adjustment') 

subplot  (2,3,4) 

title ('Trim  Tab  6-10') 

subplot (2,3,5) 

title ( ' Trim  Tab  4-10') 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

%  Capt  Nathan  A  Mille 

%  This  program  generates  the  histograms  showing  difference  between  AH- 
64D 

%  AVA  adjustments  from  PC-GBS  adjustments. 

clear;  clc; 
load  AVAadj .mat 
load  VMEPadj .mat 
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AVAcomplex  =  AVAadj .mag. * (cos (  AVAadj .phase*pi/180  )  +... 

i*sin(  AVAadj .phase*pi/180  )  ); 

VMEPcomplex  =  VMEPadj .mag . * (cos (  VMEPadj .phase*pi/180  )  +... 

i*sin(  VMEPadj .phase*pi/180  )  ); 

Magnitude  diff  =  abs (AVAcomplex  -  VMEPcomplex) ; 

Rel  mag  diff  =  Magnitude  diff . /AVAadj .mag; %max (VMEPadj .mag, 
AVAadj .mag) ; 

for  i=l : 20 

if  abs (AVAcomplex ( i , 1 ) )  <  113/2 

if  abs (VMEPcomplex ( i , 1 ) )  <  113/2 
Rel  mag  diff(i,l)=0; 

abs (AVAcomplex ( i , 1 ) ) ; 
abs (VMEPcomplex ( i ,  1 )  )  ; 
end 

end 

if  abs (AVAcomplex ( i , 2 ) )  <  0.5/2 

if  abs (VMEPcomplex ( i , 2 ) )  <  0.5/2 
Rel  mag  diff(i,2)=0; 

abs (AVAcomplex ( i , 2 ) ) ; 
abs (VMEPcomplex ( i ,  2  )  )  ; 
end 

end 

if  abs (AVAcomplex ( i , 3 ) )  <  0.5/2 

if  abs (VMEPcomplex ( i , 3 ) )  <  0.5/2 
Rel_mag_dif f ( i , 3 ) =0 ; 
end 

end 

if  abs (AVAcomplex ( i , 4 ) )  <  0.5/2 

if  abs (VMEPcomplex ( i , 4 ) )  <  0.5/2 
Rel  mag  diff(i,4)=0; 
end 

end 

if  abs (AVAcomplex ( i , 5 ) )  <  0.5/2 

if  abs (VMEPcomplex ( i , 5 ) )  <  0.5/2 

Rel_mag_dif f ( i, 5) =0 ; 

end 

end 

end 

X  =  ( [012345678  9]+0.5)/10; 
for  i=l:5 

[n  ( i ,  : ) , xout ( i ,  : ) ] =hist (Rel_mag_dif f (  :  ,  i )  ,  X) ; 
subplot ( 1 , 5 , i ) ; 

bar (xout (i,  : ) *  100 , n  ( i ,  : ) / 0 . 2 ,  ' k ' ) ; 
axis ( [0  100  0  100] ) 

end 

subplot  (1,5, 1 ) 

title ( 'Weight ' , ' Font Size ' , 16) 

ylabel (' Percentage  of  Total  Flights  Sampled' ,' FontSize ', 16) 
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subplot (1,5,2) 

title (' Pitch  Link', ' FontSize ' , 16) 
subplot (1,5,3) 

title ('Tab  8-10 FontSize ', 16) 

xlabel (' Percent  Difference  of  PC-GBS  Adjustments  from  AVA  Algorithm 
Adjustments ' , ' FontSize ' , 16) 
subplot  (1,5,4) 

title ('Tab  6-10 FontSize ', 16) 
subplot (1,5,5) 

title ('Tab  4-10 ',' FontSize ', 16) 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

%  Capt  Miller  Thesis 

%  This  is  a  test  program  to  determine  if  any  mistakes  were  made 
%  during  data  entry.  All  AH-64D  coefficients  are  plotted  together. 
Errors 

%  will  show  up  as  individual  coefficients  plotting  grossly  different 
from 

%  the  rest  of  the  coefficients, 
clear  all;  close  all;  clc; 

load  Coeff.mat 

[LSC,  standard] =Least_Squares_Coef f (Coeff) ; 

[LSCgood,  standardgood] =Least_Squares_Coef f (Coef f (1:5) ) ; 

[LSCexd,  standardexd] =Least_Squares_Coef f (Coef f (15:20) ) ; 
b=l  ; 

e=length (Coeff )  ; 
n=l ; 

for  h=l : 7 
for  j=b:e 

value (h, j  +  l-b,  : ) =Coef f { j  }  (h,  : )  ; 

end 

end 

vertvalue  =  value ( 3 : 7 ,  : ,  : )  ; 
latvalue  =  value (1:2,  :  ,  : )  ; 
for  h=l : 7 

maxi (h) =max (abs (value  (h,  : , 1) ) )  ; 
max2 (h) =max (abs (value (h,  : , 2 ) ) )  ; 
max3 (h) =max (abs (value  (h,  : , 3) ) )  ; 
max4 (h) =max (abs (value (h,  : , 4 ) ) )  ; 
max5 (h) =max (abs (value  (h,  : , 5) ) )  ; 

end 

for  h=l : 5 

vertmaxl (h) =max (abs (vertvalue (h, : , 1) ) ) ; 
vertmax2 (h) =max (abs (vertvalue (h, : , 2 ) ) ) ; 
vertmax3 (h) =max (abs (vertvalue (h, : , 3) ) ) ; 
vertmax4 (h) =max (abs (vertvalue (h,  : , 4 ) ) )  ; 
vertmax5 (h) =max (abs (vertvalue (h, : , 5) ) ) ; 

end 

for  h=l:2 

latmaxl (h) =max (abs (latvalue (h,  : ,  1) ) )  ; 
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Iatmax2 (h) =max (abs (latvalue (h,  : , 2)  )  )  ; 
latmax3 (h) =max (abs (latvalue (h, : , 3)  )  )  ; 
latmax4 (h) =max (abs (latvalue (h,  : , 4) ) )  ; 
latmax5 (h) =max (abs (latvalue (h, : , 5)  )  )  ; 

end 

for  h=l : 7 
%figure (h) 
subplot ( 3 , 3 , h) 

%polar (angle (value (h, : , 1 ) ) , abs (value (h, : , 1 ) ) /maxi (h) , ' bx ' ) 
polar (angle (value (h,  : ,  n)  )  ,  abs (value (h,  : , n) ) ,  '  rx  1  ) 

%hold  on 
%f igure  (2 ) 

%subplot ( 3 , 3 , h) 

%polar (angle (value (h, : , 2 ) ) , abs (value (h, : , 2 ) ) /max2 (h) , ' rx ' ) 

%f igure  ( 3 ) 

%subplot ( 3 , 3 , h) 

%polar (angle (value (h, : , 3) ) , abs (value (h, : , 3) ) /max3 (h) , ' gx ' ) 

%polar (angle (value (h, : , 4 ) ) , abs (value (h, : , 4 ) ) /max4 (h) , ' mx ' ) 

%polar (angle (value (h, : , 5) ) , abs (value (h, : , 5) ) /max5 (h) , ' cx ' ) 
end 

subplot (3,3, 1 ) 
title (' FPG100  Lat') 
subplot  (3,3,2) 
title ('Hover  Lat') 
subplot (3,3,3) 
title (' 60  Kts  Vert') 
subplot  (3,3,4) 
title (' 80  Kts  Vert') 
subplot (3,3,5) 
title  ('100  Kts  Vert') 
subplot (3,3,6) 
title ('120  Kts  Vert') 
subplot (3,3,7) 
title ('140  Kts  Vert') 

figure  (2 ) 
for  i=l:5 

tempi  =  max (vertmaxl ( i ) , vertmax2 ( i ) ) ; 
temp2  =  max (vertmax3 ( i ) , vertmax4 ( i ) ) ; 
temp3  =  max (tempi , temp2 ) ; 
maxval  =  max (temp3, vertmax5 (i) ) ; 

subplot (2, 3, i)  %  Each  subplot  is  for  a  different  lateral  regime 
polar ( 0 , maxval , ' w . ' ) 
hold  on 

polar (angle (vertvalue ( i , : , 1 ) ) , abs (vertvalue ( i , : , 1 ) ) , ' kx ' ) 
polar (angle (vertvalue (i, : , 2) ) , abs (vertvalue ( i , : , 2 ) ) , ' kd ' ) 
polar (angle (vertvalue (i,  : , 3) ) , abs (vertvalue ( i ,  : , 3 ) )  ,  ' kA ' ) 
polar (angle (vertvalue (i, : , 4) ) , abs (vertvalue ( i , : , 4 ) ) , ' kv ' ) 
polar (angle (vertvalue (i, : , 5) ) , abs (vertvalue ( i , : , 5 ) ) , ' kp ' ) 
end 

figure ( 3 ) 
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for  i=l:2 

tempi  =  max (latmaxl (i) , latmax2 (i) ) ; 
temp2  =  max (latmax3 (i) , latmax4 (i) ) ; 
temp3  =  max (tempi , temp2 ) ; 
maxval  =  max (temp3, latmax5 (i) ) ; 

subplot ( 1 , 2 , i )  %  Each  subplot  is  for  a  different  lateral  regime 

polar ( 0 , maxval , '  w . ' ) 
hold  on 

polar (angle ( latvalue ( i , : , 1 ) ) , abs ( latvalue (i,:,l)*100),'kx') 
polar (angle (latvalue (i,  : , 2) ) , abs (latvalue (i,  : , 2) ) ,  ' kd' ) 
polar (angle (latvalue (i,  : , 3) ) , abs (latvalue (i,  : , 3) ) ,  ' kA ' ) 
polar (angle (latvalue (i, : , 4) ) , abs (latvalue (i, : , 4) ) , ' kv' ) 
polar (angle (latvalue (i, : , 5) ) , abs (latvalue (i, : , 5) ) , ' kp ' ) 
end 

figure  ( 4 ) 

polar ( 0 , maxval , '  w . ' ) 
hold  on 

polar (angle (latvalue (i, : , 1) ) , abs (latvalue (i, : , 1) *100) , ' ko ' ) 
polar (angle (latvalue (i,  : , 2) ) , abs (latvalue (i,  : , 2) )  ,  ' kd' ) 
polar (angle (latvalue (i, : , 3) ) , abs (latvalue (i, : , 3) ) , ' kA ' ) 
polar (angle (latvalue (i, : , 4) ) , abs (latvalue (i, : , 4) ) , ' kv' ) 
polar (angle (latvalue (i,  : , 5) ) , abs (latvalue (i,  : , 5) ) ,  ' kp ' ) 
title ('Ad  Hoc  Coefficients  for  Hover (Lat)  for  AH-64A') 

abs (LSCexd.mag-LSCgood.mag) . /LSCgood.mag*100 
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